%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%lat_correction.m
%this function deucts the gravity value computed by the latitude from the
%given gravity values.
%input: gravity values, latitudes (in radians)
%output: corrected gravity values
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function lat_g = lat_correction (g, rad_latitudes)
	theory = 978032.* (1 + 5.2789e-3*(sin(rad_latitudes)).^2 - 2.35e-6*(sin(rad_latitudes)).^4);
	temp = g - theory';
	temp = temp + 975000;
	lat_g = temp;