t=0:.005:10;
n=0;
N=1000;
y=4/((2*n-1)*pi)+1.275 %*sin((n*2-1)*t);
n=1;
while n<=N
y=y+(4/((2*n-1)*pi)*sin((n*2-1)*t));
n=n+1;
end
plot(t,y);
xlabel('Seconds');
ylabel('Amplitude');