clear all
p.L=0.025;
p.C=0.025;
p.alpha=0.825;
p.beta=10;
p.r0=1.0;
p.G=0.1;
p.n=3;
xstart=zeros(3*p.n,1);
xstart(1:3)=[0.1 0 0]';
xstart(4:6)=[0.2 0 0]';
opts=odeset('RelTol',1e-3);
tend=100;
[t,x] = ode45(@mlcnRHS,[0 tend],xstart,opts,p);
xplot=x; tplot=t;

FileName="mlcn3";
clf
plot(xplot(:,2),xplot(:,5),'-r')
grid on;
set(gcf,'PaperUnits','inch','PaperPosition',[0.25 0.25 4.5 1.5])
epsfilename=sprintf('%stx1x2.eps',FileName);
print('-depsc',epsfilename);
