求助，Matlab編程解壹階非線性常微分方程組，給高分。誰知道程序怎麽修改？

clc

Yd=inline('[(YY(2)-1)*2*pi*50;(1-(1.47*1/2.82)*sin(YY(1)))/11.28]','t','YY');

[t,YY]=ode45(Yd,0.3,[31.54*pi/180;1])

x=YY(:,1);

y=YY(:,2);

plot(t,x*180/pi);

xlabel('t/s'),ylabel('delta/deg')