P=[1650 1750 1860 2070 2300 2560 3040 3710 4450 5280 6080];
Y=[1900 1910 1920 1930 1940 1950 1960 1970 1980 1990 2000];
logP=log(P)
hold off
plot(Y,logP,'*')
xlabel('Year')
ylabel('Log of World Population')
title('Growth of World Population')
pause
BestLinearFit=polyfit(Y,logP,1)
hold
plot(Y,polyval(BestLinearFit,Y))