clear functions

net=readNet('carlos_iet2011.net');

Tmax=100;
X1minmax=[0,2];
X8minmax=[0,3000];
nMC=500;

figure(1);clf
[exact,approx,MC]=plotCompare(net,2,3,Tmax,...
    {{'derivative matching',X1minmax},...
     {'MCsample',1,'X1','c-'},...
     {'dm',2,'X1','b-'},{'dm',3,'X1','g-'}},...
    {{'zero cumulants',X1minmax},...
     {'zc',1,'X1','k-'},{'zc',2,'X1','b-'},{'zc',3,'X1','g-'}},...
    {{'distribution (Monte Carlo)',X1minmax},...
     'MCdistribution',nMC,'X1'},...
    {{'derivative matching',X8minmax},...
     {'MCsample',1,'X8','c-'},...
     {'dm',2,'X8','b-'},{'dm',3,'X8','g-'}},...
    {{'zero cumulants',X8minmax},...
     {'zc',1,'X8','k-'},{'zc',2,'X8','b-'},{'zc',3,'X8','g-'}},...
    {{'distribution (Monte Carlo)',X8minmax},...
     'MCdistribution',nMC,'X8'});

print -depsc carlos_iet2011.eps