ts,printOn);
end
if(nargin==2)
    printOn=true;
end;

[A, B, C, D] = jacobFun(xequil, uequil, printOn);

% y = [ Az q alpha theta Vt Ay p r beta phi ]T
C([2:4 7:10],:) = deg2rad(C([2:4 7:10],:));
D([2:4 7:10],:) = deg2rad(D([2:4 7:10],:));

% Build Default Linear SS Model
lin_f16 = ss(A,B,C,D);
lin_f16.stateName = {'Vt', 'alpha', 'beta', ...
    'phi','theta','psi', ...
    'p', 'q', 'r',....
    'pn','pe','alt',...
    'pow'};
lin_f16.stateUnit = {'ft/s', 'rad', 'rad',...
    'rad','rad','rad',...
    'ra