octave:2> x=linspace(0,pi,100);
octave:3> y=sin(x).^50;
octave:4> plot(x,y)
octave:5> 
octave:5> y=sin(x).^50;
octave:6> plot(x,y)
octave:7> y=sin(x-0.5).^50;
octave:8> plot(x,y)
octave:9> y=sin(x-1.5).^50;
octave:10> plot(x,y)
octave:11> y=sin(x-2).^50;
octave:12> plot(x,y)
octave:13> eulerocentrato
octave:14> eulerocentrato
octave:15> laxwendroff
octave:16> laxwendroff
octave:17> laxwendroff
octave:18> quit