function g=grav_lat(phi) 
phi=phi*pi/180; 
ge=9.7803267715; 
g = ge*(1 + 0.0052790414*sin(phi).^2 + 0.0000232718*sin(phi).^4 + 0.0000001262*sin(phi).^6 + 0.0000000007*sin(phi).^8) ; 