Beetle Population with three statesL is the amount of larva populationP is the amount of pupa population
A is the amount of adult populationrestart;with(plots):Parametersb birth factorcl crowding factor for births of larva from larva populationca crowding factor for births of larva from adults population
caa crowding factor for transition from pupa to adults from adult papulation dl death rate of larvadp death rate of pupada death rate of adultsb := 7.48:
cl := 0.012:
ca := 0.009:
caa := 0.004:
dl := 0.267:
dp := 0.0:
da := 0.0036:Iinitial conditions
l0 larva
p0 pupa
a0 adultsl0 := 100:
p0 := 100:
a0 := 100:n0: Initially there are n0 iterates that are not displayed. n1: Then n1 more iterates are calculated with the values of the the iterates are displayed.n0 := 1000:
n1 := 100:f := proc(l,p,a)
local l1,p1,a1:
l1 := b*a*exp(-ca*a-cl*l):
p1 := l*(1-dl):
a1 := p*(1-dp)*exp(-caa*a) + a*(1-da):
(l1,p1,a1);
end:for i from 1 to n0 do
(l0,p0,a0) := f(l0,p0,a0):
od:printf(`\t i \t\t larva \t pupa \t adult`);
for i from 1 to n1 do
(l0,p0,a0) := f(l0,p0,a0):
printf(`\t%3d\t\t%4.2f\t\t%4.2f\t\t%4.2f\n`,i,l0,p0,a0);
od: