var('k')
g(x) = k
p(x) = 0
a0 = integrate(g(x)/pi , (x,-pi/2,pi/2))
p0 = integrate(p(x)/pi , (x,pi/2,3*pi/2))
G(x)=a0
P(x)=p0
for n in range(1,8):
        
        an1 = integrate(g(x)*cos(n*x)/pi, (x,-pi/2,pi/2))
        bn1 = integrate(g(x)*sin(n*x)/pi, (x,-pi/2,pi/2))
        an2 = integrate(p(x)*cos(n*x)/pi, (x,pi/2,3*pi/2))
        bn2 = integrate(p(x)*sin(n*x)/pi, (x,pi/2,3*pi/2))      
        G(x) = G(x)+an1*cos(n*pi/pi*x) + bn1*sin(n*pi/pi*x)
        P(x) = P(x)+an2*cos(n*pi/pi*x) + bn2*sin(n*pi/pi*x)
print G(x)+P(x)

Kreyszig-10.5-11-U (last edited 2010-12-18 07:12:12 by Shubham)