[yfchenth2]

x0=[0 ;10];
tspan=[-2*pii/10:2*pi];
[T,X]=ode45('wangshang',tspan,x0)

function dx=wangshang(t,x)
dx=[x(2);-cos(x(1))*sin(x(1))-cos(x(1))];