shiyan41.m

来自「用欧拉方法和龙格－库塔方法求微分方程数值解」· M 代码 · 共 9 行

ts=0:0.1:11;
%a=0.1745;
a=0.5236;
x0=[a,0];
[t,x]=ode23('danbai',ts,x0);
g=9.8;l=25;
z=a*cos(sqrt(g/l)*t);
[t,x(:,1),z]