[4739] in Athena Bugs
more macsyma
daemon@ATHENA.MIT.EDU (Joe Harrington)
Wed Apr 11 23:44:34 1990
Date: Wed, 11 Apr 90 23:43:33 -0400
From: Joe Harrington <jh@ATHENA.MIT.EDU>
To: bugs@ATHENA.MIT.EDU
Cc: jh@ATHENA.MIT.EDU
Reply-To: jh@ATHENA.MIT.EDU
Same thing happens using the romberg function to numerically integrate
it. It seems that it's the interaction of arrays of functions and
variables in limits of integration. This works:
test(up):=romberg(q^2,q,0,rinsurf(up)/rinsc);
but this fails:
k[0](q):=q^2;
test[l](up):=romberg(k[l](q),q,0,rinsurf(up)/rinsc);
test(0);
Bad argument to float ((mtimes simp) 3.076923076923077 ((mplus simp) 0.7 ((mtim
es simp) -0.025 ((mplus simp) -1 ((mtimes simp) 3 ((mexpt simp) $up 2))))))
--jh--