henon 映射

% The computation of Lyapunov for two-dimension Henon maping
% reference:
% 岳毅宏，韩文秀，基于系统Lyapunov指数分析的倍周期分岔研究，
% 控制与决策，17，2002: 814-819
clc;clear;close all;
% rand（‘state‘0）;
% Author‘s email:zjliu2001@163.com
% \copyright:zjliu
M=1000;
N=1000;
D2=1;
D3=0.3;
D4=0;
L1=0;
L2=0;
q=1;
for k=1:M;
    x=zeros（1N）;
    y=zeros（1N）;
    x（1）=rand;
    y（1）=rand;
    for L=1:N-1;
        x（L+1）=1-1.4.*x（L）^2+y（L）;
        y（L+1）=0.3*x（L）;
    end