548/CH6/EX6.6/6_06.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

pathname=get_absolute_file_path('6_06.sce')
filename=pathname+filesep()+'6_06data.sci'
exec(filename)
clf();

i = 1;
Cl = 0;Cd = 0;Cl_Cd =0;Thrust = 0;
while(i<=length(V)) 
    Cl(i) = 2*W/(D*S*V(i)^2);
    Cd(i) = Cdo + Cl(i)^2/(%pi*e*AR);
    Cl_Cd(i) = Cl(i)/Cd(i);
    Thrust(i) = W/Cl_Cd(i)/1000;
    Power(i)=Thrust(i)*V(i);
    R_C(i)=(Pa-Power(i))*1000*60/W;//rate of climb(R/C in meter per minute)
   i = i+1; 
end 
xlabel("Velocity (m/s)");
ylabel("Rate of climb(m/min) for CP-1");
plot2d(V,R_C,4);