Asked by Veronica
on 27 Apr 2013

From input file (the length of CL is known--26):

i = 1:length(CL); in.C_L = CL(i); in.C_D = CD(i);

From EOM Function file:

C_L = in.C_L; C_D = in.C_D;

i = 1; % this runs case 1 only B = mass./(C_D(i)*area); % ballistic coefficient rho = rho_0*exp(-h/H); % exponential atmosphere r = R + h; % radius of planet plus height L = 0.5*rho*V^2.*C_L(i)*area; % lift D = 0.5*rho*V^2.*C_D(i)*area; % drag L_D = L/D; % lift-to-drag coefficient

% Planar equations of motion Vdot = -(rho*V^2)./(2*B) - g*sin(gamma); gammadot = (((V^2*cos(gamma))/r + ((rho.*V^2)./(2*B)).*(L_D) - g*cos(gamma))/V); hdot = V.*sin(gamma); EOM = [Vdot; gammadot; hdot];

Calling ode45 from main file with constant values for y0:

[t,y] = ode45(@HW7_P4_EOM, [0:96], y0, options);

I can extract all the information needed for C_L and C_D. But I'm having trouble running each case of the 26-element vector through ode45 and saving the results from each case for plotting. I was thinking of running a for loop with the ode45 command, and using hold all for plotting, but I don't know how to vary loop through ode45 using the different parameters.

Thanks!

*No products are associated with this question.*

Opportunities for recent engineering grads.

## 0 Comments