RK45(dydt, trange, yinit, tol)：RK45就像matlab的ode45。 它可以用于八度音程，因为它...

% 函数像 ode45 一样求解常微分方程。 % dydt = 一阶微分方程组% trange = 自变量的范围% yinit = 初始条件％tol =误差绝对公差，默认值= 1e-6; ％ 例子： % [t,y] = RK45(@(t,y)[y(2); (1-y(1)^2)*y(2)-y(1)],[0 20],[2 ; 0]); % 子图 (1,2,2); % plot(t,y,'-o') % title('van der Pol 方程 (\mu = 1) with RK45'); % xlabel('时间 t'); % ylabel('解决方案 y'); ％图例（'y_1'，'y_2'）; % ％[t2，y2] = ode45（@（t，y）[y（2）;（1-y（1）^ 2）* y（2）-y（1）]，[0 20]，[2 ; 0]）; % 子图 (1,2,1); % plot(t2,y2,