simplified_main

clc, clear, close all

% SIMPLIFICATIONS
% 3 stages fixed (is it?)
% same diameter for each stage
% same thickness for each stage
% same length for each stage

% Design variables
phi = [4, 4, 4];    % Rocket stage diameter (3rd, 2nd, 1st) [m]
th = [0.005, 0.005, 0.005]; % Rocket stage wall thickness (3rd, 2nd, 1st) [m]
L = [10, 10, 10];     % Rocket stage length (3rd, 2nd, 1st) [m]
n_st = 3;   % Number of stages
n_e = 6;   % Number of engines (in the 1st stage)

% Define constant parameters
parameters;

% Objective function
f = structural_mass(phi,th,L,n_st,n_e);

% Constraints
[g,h] = constraints(phi,th,L,n_st,n_e);

% Sensitivity analysis
% gfd_mass(phi,th,L,n_st,n_e);

% Function analysis
range_phi = 12; % [m] range for variation of 1st stage diameter
step_phi = 0.01*range_phi; % [m] step for variation of 1st stage diameter
range_th = 0.05; % [m] range for variation of 1st stage thickness
step_th = 0.01*range_th; % [m] step for variation of 1st stage thickness
range_L = 50; % [m] range for variation of 1st stage length
step_L = 0.01*range_L; % [m] step for variation of 1st stage length
n_st_range = 3; % discrete number of stages (MAX IS 3)
n_e_range = 15; % discrete number of engines 
fun_analysis_mass(phi,range_phi, step_phi, th, range_th, step_th, L, range_L, step_L, n_st, n_st_range ,n_e, n_e_range)