求常微分方程组的解。 A. 建立函数文件ty.m。 function dy=ty(t,y) dy=[y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; 调用函数文件： >> [t, y]=ode45(@ty, [0, 12], [0, 1, 1]); >> plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+') B. 建立函数文件ty.m。 function dy=ty(t, y) dy=[ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; 调用函数文件： >> clear >> h=@ty; >> [t, y]=ode45(h, [0, 12], [0, 1, 1]); >> plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+') C. ty=@(t, y) [ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)]; [t, y]=ode45(ty, [0, 12], [0, 1, 1]); plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+') D. [t, y]=ode45(@(t, y) [ y(2)*y(3); -y(1)*y(3);-0.5*y(1)*y(2)], [0, 12], [0, 1, 1]); plot(t,y(:,1),'-',t,y(:,2),'*',t,y(:,3),'+')