h=0:10:100;
gamma=6.67408e-11;
m=5.9722e+24;
r=6378.;
g=gamma*m*1.e-6./(h+r).^2;

sprintf('%10.3f %10.3f \n' ,[h;g])

%[h;g]