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