using matlab Use the demographics equation derived in class to make a script which finds the population history for $4$ population groups. Assume there is no migration.

SETUP:

parameters$=[$ BR$,$d$1,$d$2,$d$3,$d$4,$t$12,$t$23,$t$34]$ where

BR$,$ birth rate in group $2$

di$,$ death rates in group Pi

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

initconds$=[$ Pinit$1,$Pinit$2,$Pinit$3,$Pinit$4]$

timedata$=[$tinit$,$tfinal,delt$]$

RESULTS:

P$1,$P$2,$P$3,$P$4,$ group time histories as row vectors. If these conventions are not followed an incorrect answer will result. $\%$ Parameters

BR$=\mathrm{.}02$;

d$1=\mathrm{.}004$;

d$2=\mathrm{.}0075$;

d$3=\mathrm{.}012$;

d$4=\mathrm{.}019$;

t$12=\mathrm{.}067$;

t$23=\mathrm{.}04$;

t$34=\mathrm{.}04$;

$\%$ Initial conditions $($in $1000\text{'}$s for convenience$)$

Pinit$1=1200$;

Pinit$2=1000$;

Pinit$3=800$;

Pinit$4=600$;

$\%$ Time parameters

tinit$=0.$;

tfinal$=100.$;

delt$=0\mathrm{.}5$;

$\%\%$ Your code goes here $\%\%$

$\%$ DISPLAY RESULTS SECTION $-----$ THIS SECTION WILL HELP YOU VISUALIZE YOUR RESULTS

figure

subplot$(2,1,1)$ plot$($t$,$P$1($:$),\text{'}-\text{'},$t$,$P$2($:$),\text{'}--\text{'},$'LineWidth',$2)$

title$(\text{'}$Demographics $1$: P$1,$ P$2$ histories, in $1000\text{'}\text{'}$s$\text{'})$

$\%$xlabel$(\text{'}$t$,$ years'$)$

ylabel$(\text{'}$P$1,$ P$2\text{'})$

legend$(\text{'}$P$1\text{'},\text{'}$P$2\text{'})$

grid on

subplot$(2,1,2)$

plot$($t$,$P$1($:$),\text{'}-\text{'},$t$,$P$2($:$),\text{'}--\text{'},$t$,$P$3($:$),\text{'}-.\text{'},$t$,$P$4($:$),\text{'}$:$\text{'},$'LineWidth',$2)$

title$(\text{'}$P$1,$ P$2,$ P$3,$ P$4$ histories, in $1000\text{'}\text{'}$s$\text{'})$

xlabel$(\text{'}$t$,$ years'$)$

ylabel$(\text{'}$P$1,$ P$2,$ P$3,$ P$4\text{'})$

legend$(\text{'}$P$1\text{'},\text{'}$P$2\text{'},\text{'}$P$3\text{'},\text{'}$P$4\text{'})$

grid on