%octave
x=linspace(0,4*pi,200);
y=x/2.*sin(x);
plot(x,y)