var('k')
f(x) = 2*k
f
a0 = integrate(f(x)/pi , (x,0,pi))
an = 0
g(x)= 0
for n in range(1,8):
        
        
        bn = integrate(f(x)*sin(n*x)/pi, (x,0,pi))      
                
        g(x) =  an*cos(n*pi/pi*x) + bn*sin(n*pi/pi*x) + g(x) 
        
print a0
print g(x)+ a0

Kreyszig-10.4-1 (last edited 2010-12-18 07:14:14 by AmitTewari)