> ;
# read `C:/Users/cgrie/Dynam Models Bio/Homeworks/HW5/M5.txt`;
# RecToSeq([1, 2, 4, 11], [0, -5, 0, -6], 1000)[1000];
# print("GrowthC");
# GrowthC([1, 2, 2, 3, 8, 2, 1, 9, 5, 7], [1, 1, 1, 1, 1, 1, 1, 1, 1, 1], 60);
# print("GrowthCe");
# GrowthCe([1, 1, 1, 1, 1, 1, 1, 1, 1, 1]);
# with(LinearAlgebra);
# f0 := [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0];
# f1 := [0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5];
# f2 := [0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25];
# fen := [op(f0), op(f1), op(f2), 0];
# survivalRates := [0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99, 0.99];
# stuff := LeslieMod(survivalRates, fen);
# GrowthCe(Stuff);
# M := Matrix(LeslieMat(survivalRates, fen));
# Largest := Eigenvalues(M);
# print("Growth rate of salmon using GrowthCe");
# GrowthCe([0, 0, 0, 0.16, 0.41]);
# PlantGseq := proc(alpha, beta, gamma, sigma, INI, K) local a, b; option remember; a := alpha*sigma*gamma; b := beta*sigma*sigma*(1 - alpha)*gamma; if K = 0 then INI[1]; elif K = 1 then INI[2]; elif K = 2 then INI[3]; else a*PlantGseq(alpha, beta, gamma, sigma, INI, K - 1) + b*PlantGseq(alpha, beta, gamma, sigma, INI, K - 2); end if; end proc;
# print("output of PlantGseq");
# seq(PlantGseq(0.5, 0.25, 2.0, 0.8, [100, 80, 80], n), n = 0 .. 20);
# print("This is in pretty good agreement with the table on page 18, but declines slightly less quickly");
# PlantGCe := proc(alpha, beta, gamma, sigma) local a, b; a := alpha*sigma*gamma; b := beta*sigma*sigma*(1 - alpha)*gamma; GrowthCe(a, b); end proc;
# PlantGCe(0.5, 0.25, 2.0, 0.8);
> NULL;
> NULL;