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

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