*d0; z3 = z0new+e3z/d3*d0;
%
% ****************************************************
% * 求李雅谱诺夫指数
% ****************************************************
%
if j>= timesum1
lsum1 = lsum1+log(d1/d0);
lsum2 = lsum2+log(d2/d0);
lsum3 = lsum3+log(d3/d0);
count1 = count1+1;
% lamda1 = lsum1/count1/h0;
% lamda2 = lsum2/count1/h0;
% lamda3 = lsum3/count1/h0;
end
end
lamda1 = lsum1/count1/h0;
lamda2 = lsum2/count1/h0;
lamda3 = lsum3/count1/h0;
% Lyapunov = [lamda1, lamda2, lamda3];
fprintf('LE1=%f,LE2=%f,LE3=%f\n',lamda1,lamda2,lamda3);
end
% *
% ****************************************************
% * 函数: 龙格库塔迭代方程
% * 描述: 龙格库塔迭代方程
% * 参数: