rhosl= .002378;
rho15=;
rho25=;
rho35=;
W=;
S=;
b=;
AR= b^2/S;
e=;
CDo=;
Tsl=;
for v=250:1250:200
CL= W*2/rho15/v^2/S;
k=1/(3.1416*e*AR);
Treq15[v]= rho15/2*v^2*S*(CDo+ k*CL^2);
Tav15[v]=rho15/rhosl*Tsl;
V=[v];
end
for v=250:1250:200
CL= W*2/rho25/v^2/S;
k=1/(3.1416*e*AR);
Treq25[v]= rho25/2*v^2*S*(CDo+ k*CL^2);
Tav25[v]=rho25/rhosl*Tsl;
V=[v];
end
for v=250:1250:200
CL= W*2/rho35/v^2/S;
k=1/(3.1416*e*AR);
Treq35[v]= rho/2*v^2*S*(CDo+ k*CL^2);
Tav35[v]=rho35/rhosl*Tsl;
V=[v];
end

pl2=plot(t1,180/pi*z1(: ,1),'k-',t1,180/pi*z1(: ,1),'ko','EraseMode','none');
	
pl1=plot(t1,180/pi*z1(: ,1),'k-','EraseMode','none');
set
