计算物理 混沌 单摆

%例1：单摆的不同周期的位移曲线function djdb[t1,w1]=ode45(@f,[0:3*pi/200:6*pi],[pi/7,0],[]);[t2,w2]=ode45(@f,[0:3*p