Labels2 = [1790:10:1860]';
Labels = [1790:10:2010]';
Pe= 330.0
a1= 1.23065
b2= 2110.5
p3= [2.913 3.929 5.308 7.239 9.638 12.866 17.069 23.191 39.818 50.189 62.947 76.212 92.228 106.021 122.775 132.164 150.697 179.323 203.302 226.545 248.709 281.421 308.745]';
Data = p3;
C = (Data(1:22,1)./Data(2:23,1))-1;
D = Data(1:22,1).*(1-(330.0./Data(2:23,1)));
Y = C\D;
b6 = Y(1,1);
a6 = 1 + 330.0./b;
P0 = Data(1,1);
FF = zeros(23,1);
FF(1,1) = P0;
for i= 2:23
EE(i,1) = (a1.*p3(i-1,1))./(1+(p3(i-1,1)./b2));
end
t = (1790:10:2010)';
a1=(1790:10:2010)';
b2= [3.929 5.308 7.239 9.638 12.866 17.069 23.191 31.443 39.818 50.189 62.947 76.212 92.228 106.021 122.775 132.164 150.697 179.323 203.302 226.545 248.709 281.421 308.745]';
plot(Labels,Data,'b-*',Labels,FF,'r-*');
title('Beverton-Holt Model, 1790-2010')
legend('Beverton-Holt Model','Census Data')
axis([1790 2010 0 500]);
ylabel('Population (Millions)');
xlabel ('Years');
hold on
SST = sum((Data-mean(Data)).^2);
SSE = sum((Data-FF).^2);
R2 = 1- (SSE/SST)