#Solving Half Range Expansions
f(x) = pi - x
c(x) = 0
s(x) = 0
for n in range(1,5):
    an = integrate((2*f(x)*cos((n*pi*x)/pi))/pi, (x,0,pi))                          #Even Periodic Extension
    bn = integrate((2*f(x)*sin((n*pi*x)/pi))/pi, (x,0,pi))                               #Odd  Periodic Extension
    
    c(x) = c(x) + an*cos((n*pi*x)/pi)
    s(x) = s(x) + bn*sin((n*pi*x)/pi)
    
print 'The Half Range Cosine expansion'
print c(x),
print "+..."
print 'The Half Range Sine expansion'
print s(x),
print "+..."

* Solution by:

Kreyszig-10.4-23-U (last edited 2010-12-18 07:19:20 by ArnabBasu)