[yfl]

function ydot=bysj(t,y)
ydot=[-t*exp(y(1)-y(2));
t*exp(y(1)-y(2))-y(2);
t*exp(y(1)-y(2))];
[t,y]=ode23t('bysj',[0,4],[0;4.2;0])
这样就可以啦。
然后你用plot（y）画图即可，我画了图，但由于附件太大不能上传。