function dy=odef(t,y);
dy(1,1)=y(2);
dy(2,1)=-9.8;