% Sine series for f(x)=(1-x) on [0,1]
n=input('Choose N');
clear c;
for j=1:n
    c(j)=2/(j*pi);
end;
xgrid=0:0.01:1;
hold off
fn=zeros(size(xgrid));
for j=1:n
    fn=fn+ c(j)*sin(j*pi*xgrid);
end
plot(xgrid, (1-xgrid), xgrid, fn);