a0,an,n,bn=var('a0,an,n,bn')
a0=1/(2*4)*(integral(4+x,x,-4,0)+integral(4-x,x,0,4))
print 'a0=',a0
an=1/4*integral((4+x)*cos(n*x*pi/4),x,-4,0)+integral((4-x)*cos(n*x),x,0,4)
print 'an=',an
bn=1/4*integral((4+x)*sin(n*x*pi/4),x,-4,0)+integral(x*sin(n*x*pi/4),x,0,4)
print 'bn=',bn

SageDays/Kreyszig-10-20 (last edited 2010-08-12 11:25:28 by Aashita)