t1=[0,0.9];

y0=[0;0];

[t,y]=ode45(@ dantong,t1,y0)

plot(t,y)