function dx=fun1(t,x)
dx = [x(2);-sin(x(1))];

