B:=k->(400/Pi)/(2*k+1); g:=(x,n)->sum(B(k)*sin((2*k+1)*x),'k'=0..n); g(x,1); plot(g(x,1),x=-Pi..Pi);