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

SageDays/Kreyszig-10-18 (last edited 2010-08-12 11:21:50 by Aashita)