danbai.m

用欧拉方法和龙格－库塔方法求微分方程数值解

function y=danbai(t,x)
g=9.8;l=25;
y=[x(2),-g/l*sin(x(1))]';