Aerospace Toolbox 

This example shows how to perform glide calculations for a Cessna 172 following Example 9.1 in reference 1 using the Aerospace Toolbox™ software.
Best glide calculations provide values (velocity and glide angle) that minimize drag and maximize liftdrag ratio (also called the glide ratio).
On this page… 

The aircraft parameters are declared as follows.
W = 2400; % weight, lbf S = 174; % wing reference area, ft^2; A = 7.38; % wing aspect ratio C_D0 = 0.037; % flaps up parasite drag coefficient e = 0.72; % airplane efficiency factor
Set the current aircraft conditions. The bank angle (phi) is zero for this case.
h = 4000; % altitude, ft phi = 0; % bank angle, deg
Convert altitude to meters using convlength. The atmospheric calculations in the next step require values in metric units.
h_m = convlength(h,'ft','m');
Calculate atmospheric parameters based on altitude using atmoscoesa:
[T, a, P, rho] = atmoscoesa(h_m, 'Warning');
Convert density from metric to English units using convdensity:
rho = convdensity(rho,'kg/m^3','slug/ft^3');
Best glide velocity is calculated using the following equation. TAS (true airspeed in feet per second) is the velocity of the aircraft relative to the surrounding air mass.
TAS_bg = sqrt((2*W) / (rho*S))... *(1./(4*C_D0.^2 + C_D0.*pi*e*A*cos(phi)^2)).^(1/4); % TAS, fps
Convert velocity from fps to kts using convvel. KTAS is true airspeed in knots.
KTAS_bg = convvel(TAS_bg,'ft/s','kts')';
Convert KTAS to KCAS using correctairspeed. KCAS (calibrated airspeed in knots) is the velocity corrected for instrument error and position error. This position error comes from inaccuracies in static pressure measurements at different points in the flight envelope.
KCAS_bg = correctairspeed(KTAS_bg,a,P,'TAS','CAS')';
Best glide angle is calculated using:
This is the angle between the flight path and the ground that provides the highest L/D ratio.
gamma_bg_rad = asin( sqrt((4.*C_D0')./(pi*e*A*cos(phi)^2 + 4.*C_D0')) );
Convert glide angle from radians to degrees using convang:
gamma_bg = convang(gamma_bg_rad,'rad','deg');
Best glide drag is calculated using:
D_bg = W*sin(gamma_bg_rad);
Best glide lift is calculated using:
L_bg = W*cos(gamma_bg_rad);
Calculate dynamic pressure using dpressure:
qbar = dpressure([TAS_bg' zeros(size(TAS_bg,2),2)], rho);
Calculate drag and lift coefficients using:
C_D_bg = D_bg./(qbar*S); C_L_bg = L_bg./(qbar*S);
Here are the best glide values:
These plots show drag and liftdrag ratio plots for the aircraft as a function of KCAS. The plots are used to verify the best glide calculations.
Set range of airspeeds and convert to KCAS using convvel and correctairspeed:
TAS = (70:200)'; % true airspeed, fps KTAS = convvel(TAS,'ft/s','kts')'; % true airspeed, kts KCAS = correctairspeed(KTAS,a,P,'TAS','CAS')'; % corrected airspeed, kts
Calculate dynamic pressure for new airspeeds using dpressure:
qbar = dpressure([TAS zeros(size(TAS,1),2)], rho);
Calculate parasite drag using:
Dp = qbar*S.*C_D0;
Calculate induced drag using:
Di = (2*W^2)/(rho*S*pi*e*A).*(TAS.^2);
Calculate total drag using:
D = Dp + Di;
Approximate lift as weight (assuming small glide angle and small angle of attack). At this speed, assuming
and using
from above, the angle of attack is about 7 degrees. Adding the flight path angle (i.e. best glide angle) from above shows the fuselage pitch (attitude angle theta) to be about 2 degrees.
L = W;
Plot L/D versus KCAS
As expected, the maximum L/D occurs at approximately the best glide velocity calculated above.
h1 = figure; plot(KCAS,L./D); title('L/D vs. KCAS'); xlabel('KCAS'); ylabel('L/D'); hold on plot(KCAS_bg,L_bg/D_bg,'Marker','o','MarkerFaceColor','black',... 'MarkerEdgeColor','black','Color','white'); hold off legend('L/D','L_{bg}/D_{bg}','Location','Best'); annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','TAS_{bg}');
Plot parasite, induced, and total drag curves
Notice the minimum total drag (i.e. D_bg) occurs at approximately the same best glide velocity calculated above.
h2 = figure; plot(KCAS,Dp,KCAS,Di,KCAS,D); title('Parasite, induced, and total drag curves'); xlabel('KCAS'); ylabel('Drag, lbf'); hold on plot(KCAS_bg,D_bg,'Marker','o','MarkerFaceColor','black',... 'MarkerEdgeColor','black','Color','white'); hold off legend('Parasite, D_p','Induced, D_i','Total, D','D_{bg}','Location','Best'); annotation('textarrow',[0.49 0.49],[0.23 0.12],'String','TAS_{bg}');
close(h1,h2);
[1] Lowry, J. T., "Performance of Light Aircraft", AIAA(R) Education Series, Washington, DC, 1999.