Another speed up: load(lrats) an replace subst with lratsubst. Since ratsubst knows a bit more about the meaning of expressions than
does subst, you'll need to decide it that's OK.
--Barton
________________________________________
Using CRE expressions reduces the run time for the big loop in your code by a factor of about 3. To do this, insert something
like
f : rat(f)
C[0] : rat(1/(2*pi^(3/2)*sqrt(rho[0, 0]*d)*m*sigma));
into your code
--Barton
________________________________________
F