Thanks, with this improvements maple still wins, but now it is only 3
times faster than maxima. ( 66 seconds against 196 seconds for N= 6 in
my loop). Is there possibility to improve this result more?
> 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
>
>
>