% SIPDB marine model % % SIPDB model for a single population % Variables S, I, D, P % Infection occurs when S contact with contact I, D and P % PARAMETERS: 3 transmission rates (contact with I, D, P) Recruitment, natural mortality, disease mortality, dead decay, dead % remove, particle release, particle removal. % % 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] % y(3): D, Dead individuals [number/m^2] % y(4): P, infectious particles (free-living pathogens, particles [number/m^3] % PAR is a structure containing the parameters for the model % % PAR_SIPDB.m sets the values of the various model parameters % global PAR % array indexes of variables nVar=4; % number of variables in the model iS=1;iI=2;iD=3;iP=4; PAR=PAR_SIPDB; % define model parameters y0=zeros(nVar,1); % initial conditions y0(iS)=100; y0(iI)=1; y0(iD)=0; y0(iP)=0; tspan=[0 400]; % time span, simulation time [t,y]=ode45(@RHS_SIPDB,tspan,y0); % Rename model results. New names for the plots. S=y(:,1);I=y(:,2);D=y(:,3);P=y(:,4);