% simul.m : simulation of the model
%%% Commentary: This function is used in the
premiermodele.m file that has
%%% to be launched.
%%% Code:
function
[y,t]=simul(t,xinit)
% [y,t]=simul(t,xinit)
% t : duration
of the simulation
% xinit:
initial state
% y: output
global Ncell
% simulation
vectol=zeros(1,8*Ncell+1)+1e-4;
options = odeset('RelTol',1e-4,'AbsTol',vectol);
[t,y] = ode23(@modele,t,xinit,options);% This function computes the output of the ode
previously defined.