李雅普诺夫指数的计算

function [Texp,Lexp]=lyapunov(n,rhs_ext_fcn,fcn_integrator,tstart,stept,tend,ystart,ioutp) %% 常微分方程系统的Lyapunov指数计算 %% n=number of nonlinear odes n是非线性的阶数 %% n2=n*(n+1)=total number of odes，这个不懂，应该是龙哥库塔方法的阶数 n1=n; n2=n1*(n1+1); %% Number of steps  round是四舍五入的算法，tend和tstart，stept是给了的东西 nit = round((tend-tstart)/stept); %% Memory allocation 内存分配，给y分配了多一点的地方，这个不懂，给cum分配和阶数一样的地方 %其他的都不懂是什么？ y=zeros(n2,1); cum=zeros(n1,1); y0=y;%这里表明y和y0其实差不多的，是比较高阶的东西，因为n2就代表这比较高阶的东西 gsc=cum; zno

李雅普诺夫指数的计算最先出现在Python成神之路。