1

代码说明：

说明：  在命令窗输入如下内容就得到图形了: [t,x]=ode45( godhua ,[0,10],[0,0,0,1]) subplot(131) plot(t,x(:,4)) xlabel( t ) ylabel( y ) subplot(132) plot(t,x(:,1)) xlabel( t ) ylabel( v ) subplot(133) plot(t,x(:,2)) xlabel( t ) ylabel( \theta ) (At the command window type the following content to be the graphics: [t, x] = ode45 (godhua, [0,10], [0,0,0,1]) subplot (131) plot (t, x (:, 4 )) xlabel (t) ylabel (y) subplot (132) plot (t, x (:, 1)) xlabel (t) ylabel (v) subplot (133) plot (t, x (:, 2)) xlabel ( t) ylabel (heta))