Beetle Population with three states L is the amount of larva population P is the amount of pupa population A is the amount of adult population restart;with(plots): Parameters b birth factor cl crowding factor for births of larva from larva population ca 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 larva dp death rate of pupa da death rate of adults b := 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 adults l0 := 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: