clear functions

net=readNet('reversible.net');

Tmax=1;
X1minmax=[0,20];
X2minmax=[0,20];
X3minmax=[10,30];
nMC=400;

figure(1);clf
[exact,approx,MC]=plotCompare(net,3,5,Tmax,...
    {{'derivative matching',X1minmax},{'MCsample',1,'X1','c-'},{'dm',2,'X1','b-'}},...
    {{'zero cumulants',X1minmax},{'zc',1,'X1','k-'},{'zc',2,'X1','b-'}},...
    {{'low dispersion',X1minmax},{'ld',1,'X1','k-'},{'ld',2,'X1','b-'}},...
    {{'quasi deterministic',X1minmax},{'qd',1,'X1','k-'},{'qd',2,'X1','b-'}},...
    {{'distribution (Monte Carlo)',X1minmax},'MCdistribution',nMC,'X1'},...
    {{'derivative matching',X2minmax},{'MCsample',1,'X2','c-'},{'dm',2,'X2','b-'}},...
    {{'zero cumulants',X2minmax},{'zc',1,'X2','k-'},{'zc',2,'X2','b-'}},...
    {{'low dispersion',X2minmax},{'ld',1,'X2','k-'},{'ld',2,'X2','b-'}},...
    {{'quasi deterministic',X2minmax},{'qd',1,'X2','k-'},{'qd',2,'X2','b-'}},...
    {{'distribution (Monte Carlo)',X2minmax},'MCdistribution',nMC,'X2'},...
    {{'derivative matching',X3minmax},{'MCsample',1,'X3','c-'},{'dm',2,'X3','b-'}},...
    {{'zero cumulants',X3minmax},{'zc',1,'X3','k-'},{'zc',2,'X3','b-'}},...
    {{'low dispersion',X3minmax},{'ld',1,'X3','k-'},{'ld',2,'X3','b-'}},...
    {{'quasi deterministic',X3minmax},{'qd',1,'X3','k-'},{'qd',2,'X3','b-'}},...
    {{'distribution (Monte Carlo)',X3minmax},'MCdistribution',nMC,'X3'});


