求問壹道微分方程組數值解法的問題，如果有大神做出，matlab代碼發

clc

phi=0;F=7500;C=2940;u=1;

f=@(t,y)([y(2);F/y(5)*sin(phi)-u/y(1)^2+y(1)*y(4)^2;y(4);-(F/y(5)*cos(phi)+2*y(2)*y(4))/y(1);-F/C]);

[t,Y]=ode45(f,[0 7500],[1752 0 0 9.7e-6 2400]);

plot(t,Y(:,1),t,Y(:,2),t,Y(:,3),t,Y(:,4),t,Y(:,5));

legend('r','v','xita','omiga','m')