function pendulum=f(t,y)
pendulum= [y(2); -9.8*sin(y(1))];
%pendulum=[0 1;-9.8 0]*y;
%pendulum= [y(2); -9.8*y(1)];