pathname=get_absolute_file_path('602.sce')
filename=pathname+filesep()+'602data.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;
    Tf(i)=2*16245/1000;
   i = i+1; 
end 
xlabel("Velocity (m/s)");
ylabel("Thrust (kN)");
plot2d(V,Thrust,3);
plot2d(V,Tf,4);
disp("As Thrust required equals Thrust provided by two turbofan at Velocity 297 m/s approx(intersection point of both graphs.)so it will be Vmax")
Vmax=297;
printf("\Answer:\n")
printf("\n\maximum velocity: %f m/s\n\n",Vmax)