lyaplogistic.m

来自「非线性动力学中lyapunov、duffing等方程的求解」· M 代码 · 共 18 行

x=rand(1,1)
save xxx x
lambda=4;
sum=0;
N=5000;
for i=1:N
    f=lambda*x*(1-x);
    df=lambda*(1-2*x);
    sum=sum+log(abs(df));
    x=f;
end
h=sum/N