超混沌迭代代码

function F=rk3（hTXf）
%h是步长
%T仿真总时间
%X给迭代变量赋初值
%f是函数
N=T/h;         %仿真总次数
d=length（X）;    %状态变量的维数
x1=zeros（N1）;   %状态变量d维，每维N个迭代值
x2=zeros（N1）;
x3=zeros（N1）;
x4=zeros（N1）;
% dx1=zeros（N1）;   %状态变量d维，每维N个迭代值
% dx2=zeros（N1）;
% dx3=zeros（N1）;
% dx4=zeros（N1）;
g0=zeros（1d）; %存储每次迭代产生的g0值的数组，大小为d维
g1=zeros（1d）;   
g2=zeros（1d）;   
g3=zeros（1d）;   
g4=zeros（1d）;
g5=zeros（1d）; 
i=1;             %初次迭代时，数组元素序号为1
for t1=0:h:T     %循环迭代
    X=X+h.*g5;   %求