g(x) = (x**2)/2
a0 = integrate(g(x)/pi , (x,-pi,pi))
G(x)=a0
for n in range(1,8):
        
        an = integrate(g(x)*cos(n*x)/pi, (x,-pi,pi))
        bn = integrate(g(x)*sin(n*x)/pi, (x,-pi,pi))    
        G(x) = G(x)+an*cos(n*pi/pi*x) + bn*sin(n*pi/pi*x)
print G(x)

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