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