发布网友 发布时间:2022-04-28 20:12
共1个回答
热心网友 时间:2022-06-23 01:24
你把下面的程序放到第一个m文件中
t0=[0,2];
x0=[0,0,12,5,5,0];
options = odeset('RelTol',1e-8,'AbsTol',1e-8,'OutputFcn',@odeplot);
[t,y] = ode45(@odefunc,t0,x0,options);
figure(2)
plot3(y(:,1),y(:,2),y(:,3))
把下面的程序放在第二个m文件中
function dy=odefunc(t,y)
g=-9.8;
dy=[y(4)
y(5)
y(6)
0
0
g];