Fourier series- Square wave - Sawtooth wave

#Fourier Series-Square Wave
#f(t)=a/2+(2a)/pi[sin(wt)+sin(3wt)/3+sin(5wt)/5+........]

from pylab import*
t=linspace(0.0,2*pi,100)
a=2.0
w=2.0
y=a/2
for i in range(1,200,2):
    y=y+2*a*sin(i*w*t)/(pi*i)
plot(t,y)
show()

#Fourier Series-Saw-tooth wave

#f(t)=a/2+a/pi[sin(wt)+sin(2wt)/2+sin(3wt)/3+sin(4wt)/4+sin(5wt)/5+........]

from pylab import*

t=linspace(0.0,2*pi,100)

a=2.0

w=5.0

y=a/2

for i in range(1,100):

    y=y+a*sin(i*w*t)/(pi*i)

plot(t,y)

show()