Creates matrix population models for use in population ecological
analyses. Specializes on the construction of historical matrices,
which are 2d matrices comprising 3 consecutive times of demographic
information. Estimates both raw and function-based forms of
historical and standard ahistorical matrices. It also estimates
function-based age-by-stage matrices. Methodology based on
Ehrlen (2000)