[4738] in Athena Bugs
macsyma
daemon@ATHENA.MIT.EDU (Joe Harrington)
Wed Apr 11 22:52:05 1990
Date: Wed, 11 Apr 90 22:51:03 -0400
From: Joe Harrington <jh@ATHENA.MIT.EDU>
To: bugs@ATHENA.MIT.EDU
Cc: jh@ATHENA.MIT.EDU
Reply-To: jh@ATHENA.MIT.EDU
The following three macsyma batch files should produce the same
output. The first generates a floating exception, the other two work.
They are part of a larger program which needs the sophistication
of the first example. Since this is a bug in third-party software, I
expect you will need to send it to Symbolics.
------example 1 (fails)------
showtime:true;
numer:true;
p[-1](x):=0;
p[0](x):=1;
p[1](x):=x;
p[l](x):=ratsimp((2*l-1)/(l)*x*p[l-1](x)-(l-1)/(l)*p[l-2](x));
rinsurf(u):=7/10-p[2](u)/20;
rinsc:rinsurf(1)/2;
ingrand[l](q,up):=q^(l+2);
inint[l](up):=integrate(ingrand[l](q,up),q,0,rinsurf(up)/rinsc);
------example 2 (works)------
showtime:true;
numer:true;
ingrand[l](q,up):=q^(l+2);
inint[l](up):=integrate(ingrand[l](q,up),q,0,29/13);
------example 3 (works)------
showtime:true;
numer:true;
p[-1](x):=0;
p[0](x):=1;
p[1](x):=x;
p[l](x):=ratsimp((2*l-1)/(l)*x*p[l-1](x)-(l-1)/(l)*p[l-2](x));
rinsurf(u):=7/10-p[2](u)/20;
rinsc:rinsurf(1)/2;
ingrand(q):=q^(2);
inint(up):=integrate(ingrand(q),q,0,rinsurf(up)/rinsc);
--jh--
