/////////////////////////
// 
// reacteur biologique: 
// calcul de la quantité de pbioogaez produite
//
/////////////////////////
//  
//
kappa=evstr(x_dialog('valeur de kappa','0.5'))
uinit=evstr(x_dialog('quantité initiale de micro-organisme','0.1'))
N=1000;  // nombre de points calculés
T=8;     // temps final
//
h=T/N;t=0:h:T;
//
u0=uinit;u=[uinit];
S=0;
for i=1:N
   uold=u0; 
   u0=u0+h*u0/kappa*(1-u0-S)
   S=S+h/2*(uold+u0)
   u=[u,u0];
end
//
figure(1)
plot(t,u,'black','LineWidth',3)
//////////////////////
//
disp('production totale de biogaz:');disp(S)
