% Loop invoked by ex_6_5.m kk = 1; warning off for ii = 0:5; for jj = 0:5; % Trim parameters Al = ii*(0.349/5); % from 0 to 0.349 Ma = 2 + jj*(2/5); % from 2 to 4 % Computation of equilibrium points X0 = [0;0;Al;0]; [X,U,Y,dX] = trim('missile53',X0,0,0,3); % Computation of linearized models [A,B,C,D] = linmod('missile53',X,U); B = B/B(1,1); C = C*B(1,1); % Storage of values for state-space model interpolation abcd_data{kk} = [A B;C D]; al_ma_data{kk} = [Al,Ma]; % Storage of values for equilibrium surface interpolation q_dp_data{kk} = [X(4) U]; % values of q and dp kk = kk + 1; end; end; warning on