function dydt=de(t,y)
dydt=y*(1-t^2)+2;