Modify the model derived in class to assess the following situation. Suppose an epidemic (which shall remain nameless) raises all group death rates for given period of time (i.e. it works much like the pulse function from last week & problem 2). This epidemic also lowers b2 birthrate. Assume migration is constant for group P2.

SETUP:

parameters= [ b2,d1,d2,d3,d4,t12,t23,t34] where

b2, birth rate in group 2

di, death rates in group Pi

tij, transition from group Pi to group Pj due to ageing

Epieffect= [b2Epifactor,drateEpifactor,tStartEpi,tEndEpi], where

b2Epifactor, b2Epi= b2*b2Epifactor

drateEpifactor, d1Epi= d1*drateEpifactor, etc.

M2ann, constant migration rate for P2

initconds= [ Pinit1,Pinit2,Pinit3,Pinit4]

timedata= [tinit,tfinal,delt]

RESULTS:

P1,P2,P3,P4, group time histories

This epidemic also lowers b2 birthrate. Assume migration is constant for group P2. SETUP: - parameters =[b2,d1,d2,d3,d4,t12,t23,t34] where - b2, birth rate in group 2 - di, death rates in group Pi - tij, transition from group Pi to group Pj due to ageing - Epieffect= [b2Epifactor,drateEpifactor,tStartEpi,tEndEpi], where - b2Epifactor, b2Epi= b2*b2Epifactor - drateEpifactor, d1Epi= d1*drateEpifactor, etc. - M2ann, constant migration rate for P2 - initconds =[ Pinit1,Pinit2,Pinit 3, Pinit4 ] - timedata =[ tinit,ttinal,delt ] RESULTS: - P1,P2,P3,P4, group time histories Function 8 1 function [P1,P2,P3,P4]= Demographics_3_fcn(parameters,Epieffect,M2ann,initconds,timedata) end Sode to call your function (?)