function LS = malsqsim(p)
%Sum of square errors: Malthusian growth model
t = [0:10:220];
y1 = [3.929 5.308 7.24 9.638 12.866 17.069 23.192 31.443 39.818];
y2 = [50.189 62.948 76.212 92.228 106.022 122.775 132.165 150.697];
y = [y1,y2,179.323 203.302 226.546 248.71 281.422 308.746];  
LS = sum((p(1)*exp(p(2)*t)-y).^2);
end

