Team:Paris/Modeling/simul

From 2008.igem.org

Revision as of 15:46, 8 August 2008 by YannLC (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)

To go back to the previous page : Clic here.


% simul

% 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.