On Dec 4, 1:58 am, Robert Dodier <[email protected]> wrote: > OK, here's the input I tried: > > kk:1; > ll:10; > dz:1; > ghomo(zl,zg):= 2/ll*sin(kk*zl)*sin(kk*(ll-zg))/(kk*sin(kk*ll)) ; > gfourn(z,zd,n):=sin(npi*n*z/ll)*sin(npi*n*zd/ll)/((npi*n/ll)^2 - > kk^2) ; > uhomo(z):=integrate(ghomo(zl,z)*f1(zl),zl,0,z)+integrate(ghomo(z,zg)*f1 > (zg),zg,z,ll); > ufourn(z,n):=makelist(integrate(f1(zd)*gfourn(z,zd,nn),zd,0,ll),nn, > 1,n); > f1(z):=z*(ll-z); > npi:%pi; > > and with this I get: (in Maxima built from CVS recently) > > ufourn(23/10,8); > => [4000*sin(23*%pi/100)/(%pi^3*(%pi^2/100-1)),0, > 4000*sin(69*%pi/100)/(27*%pi^3*(9*%pi^2/100-1)),0, > 32*sin(23*%pi/20)/(%pi^3*(%pi^2/4-1)),0, > 4000*sin(161*%pi/100)/(343*%pi^3*(49*%pi^2/100-1)),0] > float (%); > => [-94.65540344672326,0.0,-35.36737678051993,0.0, > -.3192995340015049,0.0,-0.0922486943844007,0.0] >
This also works in maxima_console() in Sage 4.2.1, so the problem must lie in something else Sage has done to create additional contexts. I'll try to look into this more. - kcrisman -- To post to this group, send an email to [email protected] To unsubscribe from this group, send an email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org
