//
// Simulation d'une file M/M/1/R
//
l = 1;                          // taux d'arrivee 
m = 1.001;                        // taux de depart
R = 50;                        // capacite maximale
l = [0,l,zeros(1,R-1)];
m = [0,m,zeros(1,R-1)];
L = toeplitz(m,l);
L = L-diag(sum(L,"c"));         // generateur
x0 = 1;                         // etat initial
n = 20000;                      // nombre de pas
getf("saut.sci");
//
// Simulation par la chaine incluse
//
traj = saut_incluse(L,x0,n);
dates = [0,cumsum(traj(2,:))]'; // dates de saut
etats = [x0,traj(1,:)]'-1;      // etats visites
xset("window",0); xbasc();      // graphique
plot2d2("gnn",dates,etats,1,"111",..
          "Nombre de clients",[0,0,max(dates),max(etats)]);
//
// Simulation par la chaine harmonisee
//
traj = saut_harmonisee(L,x0,n);
dates = [0,cumsum(traj(2,:))]'; // dates de saut
etats = [x0,traj(1,:)]'-1;      // etats visites
xset("window",1); xbasc();      // graphique
plot2d2("gnn",dates,etats,1,"111",..
          "Chaine harmonisee",[0,0,max(dates),max(etats)]);

