x = 0:0.01:pi;

figure(figure)

plot(x,x.*sin(x),'.g');
hold on

y = pi*sin(x)/2;
plot(x,y,'-r');


for lcv = 2:2:10,
  pause
  y = y - (8*lcv/(pi*(lcv*lcv-1)^2)).*sin(lcv.*x);
  plot(x,y,'-r');
end