    % Roll angle from horizontal guess (deg)
    thetag = 0;     % Pitch angle guess (deg)
    xguess = [Vtg 0 0 phig thetag 0 0 0 0 0 0 hg 0];

    % u = [throttle elevator aileron rudder]
    uguess = [.2 0 0 0];

    % Orientation for Linearization
    % 1:    Wings Level (gamma = 0)
    % 2:    Wings Level (gamma <> 0)
    % 3:    Constant Altitude Turn
    % 4:    Steady Pull Up
    orient = 4;
    inputs = [xguess(1), xguess(12), 0, 0, 0];
   
    [xequil,uequil] = trimmerFun(xguess,uguess,orient,inpu