Demonstrate LVLH angular rate for an eccentric orbit

The nominal rotation is constant about Z. The inclination adds a Y component. For an eccentric orbit there is also a periodic change in the components of omega.

%--------------------------------------------------------------------------
%   Copyright (c) 2011, 2012 Princeton Satellite Systems, Inc.
%   All rights reserved.
%--------------------------------------------------------------------------

[r,v] = RVFromKepler([7500 0.0 0 0 0.0 0]);
omega = OmegaLVLH( r, v );
Plot2D(1:size(omega,2),omega,'Sample','Omega','LVLH Rate Vectors - Circular');

[r,v] = RVFromKepler([7500 0.2 0 0 0.0 0]);
omega = OmegaLVLH( r, v );
Plot2D(1:size(omega,2),omega,'Sample','Omega','LVLH Rate Vectors - Inclined');

[r,v] = RVFromKepler([7500 0.2 0 0 0.1 0]);
omega = OmegaLVLH( r, v );
Plot2D(1:size(omega,2),omega,'Sample','Omega','LVLH Rate Vectors - Eccentric');

%--------------------------------------