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