Thread Subject:

A problem in plotting the fourier series !!!!

syms an bn a0;
syms t n;
syms x;
y = sym(0);

T=2;
N=13;
x = heaviside(t+0.5)-heaviside(t-0.5); %Square pulse
f = @(t) x;

an = int(f(t).*cos(n*pi*t/T),t,-T,T)*(1/T);
bn = int(f(t).*sin(n*pi*t/T),t,-T,T)*(1/T);
a0 = int(f(t),t,-T,T)*(1/T);
y = a0;

pretty(a0);
pretty(an);
pretty(bn);


for i=-N:2:N;
     y = y + subs(an, i) + subs(bn, i);
end


ezplot(i,y), grid on % Problem point

I am not able to print 'y' which is the fourier series of the square wave !!!!
Please help !!!!