diff(sin(x), x, 2); plot(sin(x), x, -2pi, +2pi); do start:=time(); for(p:=3; p<50000; p:=p+2) do y:=1; for(r:=3; r