// exec('TP_plante_1.sce')
// paramètres
//puis=1; //  gamma
puis=0.5; //  gamma
beta=0.9  // probabilité de survie
r=1.2;
T=10; 
x0=1;



getf('TP_plante_1.sci');
strategie=list();
strategie(1)=strategie_M;
strategie(2)=strategie_R;
strategie(3)=strategie_C;

correspondance=list();
correspondance(1)=string("mourir");
correspondance(2)=string("aleatoire");
correspondance(3)=string("croitre");

strategie(4)=strategie_O;
correspondance(4)=string("optimale");

xp=(beta*r*puis)^{1/(1-puis)};
xm=(xp/r)^{1/puis};

// dans ce qui suit, nous adoptons les notations des arguments 
// de la fonction traj_feedback
etat_initial=1;
dynamique=dyn_plante_com;
alea=ceil(beta-rand(1,T));
// T réalisations indépendantes d'une binomiale B(beta,1)
cout_inst=rejetons;
cout_fin=cout_fin_zero;

for i=1:4 do
[b,v]=traj_feedback...
          (etat_initial,dynamique,strategie(i),alea);
[f]=cout_total(b,v,cout_inst,cout_fin);

 xset("window",3*(i-1));xbasc();plot2d2(1:T+1,f,rect=[0,0,T+2,max(f)+1]);
 xtitle("stratégie "+correspondance(i)+" : fitness");
 
 xset("window",3*(i-1)+1);xbasc();plot2d2(1:T,v,rect=[0,0,T+1,max(v)+1]);
 xtitle("stratégie "+correspondance(i)+" : commande");
 
 xset("window",3*(i-1)+2);xbasc();plot2d2(1:T+1,b,rect=[0,0,T+2,max(b)+1]);
 xtitle("stratégie "+correspondance(i)+" : état"); 
 
 printf("la fitness totale pour la strategie "...
        +correspondance(i)+" est : %"+" f\n", f($))
 
 halt();
 xdel((3*(i-1)):(3*(i-1)+2));
end


P=0.1:0.2:0.9;
for j=1:prod(size(P)) do
  puis=P(j);
  xp=(beta*r*puis)^{1/(1-puis)};
  xm=(xp/r)^{1/puis};
x0=xm/10;
// taille initiale inférieure au seuil xm
printf("gamma=%"+" f\n", P(j));
     for i=1:4 do
     [y,v]=traj_feedback(x0,dynamique,strategie(i),alea);
     [f]=cout_total(y,v,cout_inst,cout_fin);
     printf("la fitness totale pour la stratégie "...
        +correspondance(i)+" est : %"+" f\n", f($))
     end
end