求微分方程数值解、画曲线

求数值解、画出（y（t），y'（t））的曲线
条件
y''（t）= -sin（y（t）），t∈[0,10]
y（0）= 3
y'（0）= 0

function  Test1()

y0 = [3 0]';
[T,Y] = ode45(@Dfun,[0 10],y0);
figure
plot(T,Y)

function dy = Dfun(t,y)

dy = [y(2);  -sin( -y(1) )];

结果：