%CALCULATION OF LYAPUNOV EXPONENTS FOR THE MLC SYSTEM
clear all; p.L=0.025; p.C=0.025; p.alpha=0.825; p.beta=10; p.r0=1.18;
mlcle(200.0,1000.0,0.1,odeset('RelTol',1e-3),p,[6.374674155 3.16227766 1.169288347]');
%Results: 
%Step: 0.10000, t: 1000.00000, 1.66460, 0.00920, -27.13465

clear all; p.L=0.025; p.C=0.025; p.alpha=0.825; p.beta=10; p.r0=1.18;
mlcle(200.0,1000.0,0.2,odeset('RelTol',1e-3),p,[6.374674155 3.16227766 1.169288347]');
%Results: 
%Step: 0.20000, t: 1000.00000, 1.68016, 0.00191, -27.16149

clear all; p.L=0.025; p.C=0.025; p.alpha=0.825; p.beta=10; p.r0=1.18;
mlcle(200.0,1000.0,0.5,odeset('RelTol',1e-3),p,[6.374674155 3.16227766 1.169288347]');
%Results: 
%Step: 0.50000, t: 1000.00000, 1.64175, 0.02229, -27.12453

%CALCULATION OF LYAPUNOV EXPONENTS FOR THE GENERIC VARIATIONAL EQUATION FOR THE MLC SYSTEM

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,101,200.0,200.0,0.5,odeset('RelTol',1e-4),p,[6.374674155 3.16227766 1.169288347]','mlc118legen1.res');
% resutls in mlc118legen1.res

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,101,200.0,200.0,0.5,odeset('RelTol',1e-4),p,[-4.442057487 -3.16227766 -0.01408704155]','mlc118legen2.res');
% resutls in mlc118legen2.res

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,101,500.0,1000.0,0.5,odeset('RelTol',1e-4),p,[6.374674155 3.16227766 1.169288347]','mlc118legen1a.res');
% resutls in mlc118legen1a.res


clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,101,500.0,1000.0,0.5,odeset('RelTol',1e-4),p,[-4.442057487 -3.16227766 -0.01408704155]','mlc118legen2a.res');
% resutls in mlc118legen2a.res

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,11,200.0,200.0,0.5,odeset('RelTol',1e-4),p,[-1.451976737 -3.16227766 -1.850734763]','mlc118legen3.res');
% resutls in mlc118legen3.res

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,11,200.0,200.0,0.5,odeset('RelTol',1e-4),p,[1.376941029 3.16227766 -1.949889025]','mlc118legen4.res');
% resutls in mlc118legen4.res

clear all; p.L=0.025;  p.C=0.025;  p.alpha=0.825;  p.beta=10;  p.r0=1.18;  
mlclegen(-0.2,0.0,101,500.0,1000.0,0.5,odeset('RelTol',1e-4),p,[-1.451976737 -3.16227766 -1.850734763]','mlc118legen3a.res');
% resutls in mlc118legen3a.res
