李雅普诺夫指数的计算

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; znorm=cum;
%% Initial values分配初值，把所有的初值都变成一列都给了y的第一行，当阶数为一的时候y（4）=1=y（8）=y（10）
y(1:n)=ystart(:);
for i=1:n1
    y((n1+1)*i)=1.0;
end;
t=tstart;%给t分配一个起始值
%% Main loop  把时间的次数分配个iITERLYAP
for ITERLYAP