% SI marine model % % SI model for a single population % % The model considers a number of individuals in a population % The model considers a number of individuals in a (surface) % % Units: % time is in days % populations are in number of individuals (/m^2) % % variables are % y(1): S, susceptible, uninfected individuals [number/m^2] % y(2): I, infected individuals [number/m^2] % PAR is a structure containing the parameters for the model % % PAR_SI.m sets the values of the various model parameters % global PAR % array indexes of variables nVar=2; % number of variables in the model iS=1;iI=2; PAR=PAR_SI; % define model parameters y0=zeros(nVar,1); % initial conditions y0(iS)=100; y0(iI)=1; tspan=[0 400]; % time span, simulation time %tspan=[0 350]; % time span sumulation time [t,y]=ode45(@RHS_SI,tspan,y0); % rename model results S=y(:,1);I=y(:,2);