a0,an,bn,n,k =var('a0,an,bn,n,k')
# f(x) is even
a0=1/(pi)*(integral(k,x,-pi/2,pi/2))
print 'a0=',a0
an=integral(k*cos(n*x),x,-pi/2,pi/2)*2/pi
print 'an=',an

SageDays/Kreyszig-10.4-11 (last edited 2010-08-12 13:33:41 by Aashita)