format long
x = [1960 1970 1990 2000];
y = [3039585530 3707475887 5281653820 6079603571];
A = [ones(4,1) x']

A =

           1        1960
           1        1970
           1        1990
           1        2000

c = A\log(y)'

c =

 -12.26234902764795
   0.01740324629807

exp(-12.26234902764795+0.01740324629807*1980)-4452584592

ans =

    -9.109867659928703e+07

diary off