计算物理解常微分方程编程.pdfVIP

例 1： function jyjs [T,Y]=ode45(yjs,[0:10],[2,1]); plot(T,Y(:,1),T,Y(:,2)) function ydot=yjs(t,y) ydot=[y(2); 4 ]; 例 2：研究受空气阻尼的抛体运动 function znxp global m d e % 用全局变量传递参数 m=1;b=[0,0.2,0.2];p=[0,0,1]; % 传递给方程的参数 px=[4.6;4.5;4.5]; % 说明文字的 x 坐标