clear all
t=[1 3 5 7 11 12 17 18]';
N=[263 357 439 483 693 773 1444 1249]';
X=[ones(size(t)) t]
c=X\log(N)
z=exp(X*c)

plot(t,N,'r*',t,z,'-ob')
xlabel('years after 1960')
ylabel('size in thousands')
legend('data','fitted model')