代碼如下:
x0?=?[1;1];?%初始條件
%求出t在[0,50]時,對應x,y的曲線
[t,xx]?=?ode45(y,[0,50],x0);
plot(t,xx(:,1),t,xx(:,2))
legend('x','y')
圖形結果如下: