function LS = mallstsqc(p0,r,t,y)
%Least Squares sum of square errors to Malthusian growth model
LS = sum((p0*exp(r*t)-y).^2);
end

