read `DMB.txt`: 
print(`This is Proj4.txt, a Maple package accompanying Project 4 of Dr. Z.'s Math 336 (Dynamical Models in Biology), Fall 2025`):
print(` Project Leader: Daniyal Chaudhry `):
print(``):
print(`Other members : Katryn Keating, Zach Ali, Ezra Chechik`):
print(``):
print(`For a list of the procedures done by this team type: Help4(); For Help with a specific procedure, type: Help4(ProcedureName);`):

Help4 := proc()
if nargs = 0 then
print(`The available procedures are: HWgUAA, HWgUAa, HWgUaa, HWgU `):

elif nargs = 1 and args[1] = HWgU then
print(`HWgU(a,b,c): Uses a generalized version of the Hardy-Weinberg equation where:`):
print(`a represents universal mating preference for the AA genotype,`):
print(`b represents universal mating preference for the Aa genotype,`):
print(`c represents universal mating preference for the aa genotype,`):
print(`such that the matrix used in calculation is [[a,b,c],[a,b,c],[a,b,c]]`):
print(`Outputs a vector with three components, which are the Stable Steady State values of u, v, and w, in that order.`):
print(`Try:`):
print(`HWgU(2,1,1);`):

elif nargs = 1 and args[1] = HWgUAA then
print(`HWgUAA(a,b,c): Uses a generalized version of the Hardy-Weinberg equation where:`):
print(`a represents universal mating preference for the AA genotype,`):
print(`b represents universal mating preference for the Aa genotype,`):
print(`c represents universal mating preference for the aa genotype,`):
print(`such that the matrix used in calculation is [[a,b,c],[a,b,c],[a,b,c]]`):
print(`Outputs a single value, which is the Stable Steady State value of u (AA proportion).`):
print(`Try:`):
print(`HWgUAA(2,1,1);`):

elif nargs = 1 and args[1] = HWgUAa then
print(`HWgUAa(a,b,c): Uses a generalized version of the Hardy-Weinberg equation where:`):
print(`a represents universal mating preference for the AA genotype,`):
print(`b represents universal mating preference for the Aa genotype,`):
print(`c represents universal mating preference for the aa genotype,`):
print(`such that the matrix used in calculation is [[a,b,c],[a,b,c],[a,b,c]]`):
print(`Outputs a single value, which is the Stable Steady State value of v (Aa proportion).`):
print(`Try:`):
print(`HWgUAa(2,1,1);`):

elif nargs = 1 and args[1] = HWgUaa then
print(`HWgUaa(a,b,c): Uses a generalized version of the Hardy-Weinberg equation where:`):
print(`a represents universal mating preference for the AA genotype,`):
print(`b represents universal mating preference for the Aa genotype,`):
print(`c represents universal mating preference for the aa genotype,`):
print(`such that the matrix used in calculation is [[a,b,c],[a,b,c],[a,b,c]]`):
print(`Outputs a single value, which is the Stable Steady State value of w (aa proportion).`):
print(`Try:`):
print(`HWgUaa(2,1,1);`):

else
print(`There is no Help for`, args):
fi:
end:

HWgUAA := proc(a, b, c) SSSgN(HWg(u, v, [[a, b, c], [a, b, c], [a, b, c]]), [u, v])[1]; end: 
HWgUAa := proc(a, b, c) SSSgN(HWg(u, v, [[a, b, c], [a, b, c], [a, b, c]]), [u, v])[2]; end: 
HWgUaa := proc(a, b, c) 1 - HWgUAA(a, b, c) - HWgUAa(a, b, c); end: 
HWgU := proc(a, b, c) [HWgUAA(a, b, c), HWgUAa(a, b, c), HWgUaa(a, b, c)]; end: