%% Calculate centrifugal effect3
%% Calculate centrifugal effect3
% Calculate the precessing centrifugal effect of gravity for Earth at 15000 m over the equator and 11000 m over the North Pole.

p       = [2412.648e3 -2412.648e3 0; 0 0 3376e3]
% Set julian date to January 1, 2000 at noon GMT
JD      = 2451545
% Calculate precession rate in right ascension in meters
pres_RA = 7.086e-12 + 4.3e-15*(JD - 2451545)/36525
% Calculate the rotational rate in a precessing reference
% frame
Omega   = 7.2921151467e-5 + pres_RA
[gx, gy, gz] = gravitycentrifugal( p, 'custom', Omega )

% Copyright 2015 The MathWorks, Inc.