[4739] in Athena Bugs

home help back first fref pref prev next nref lref last post

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--

home help back first fref pref prev next nref lref last post