[图],可行的是()A、建立函数文件ty.m function?dy=ty(...
,可行的是()
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),'+')