function [Dactual,Jactual,nactual,RPMactual,Etaactual,Poutactual,Ppinactual,ifig]=...
DesignProp(rho,Pop,Vop,Jstar,CTstar,CPstar,Etastar,J,CT,CP,ifig)

Dstar=sqrt(Pop/(Vop*rho*CTstar*(Vop/Jstar)^2)); % Propeller diameter (ft)
Dstarinches=Dstar*12;  % Propeller diameter (inches)
nstar=Vop/(Jstar*Dstar); % propeller speed (Hz)
RPMstar=nstar*60; % Propeller revolutions per minute
disp(' '); disp('PRELIMINARY PROPELLER ANALYSIS');
string2=['  If the propeller of diameter ',num2str(Dstarinches),' inches'];
string3=[' spins at ',num2str(RPMstar),' rpm,'];
string4=['  the propeller efficiency will be maximum at Eta*= ',num2str(Etastar),','];
string5=['  at the airspeed of ',num2str(Vop),' ft/sec,'];
string6=['  and the prop will be generating the required ',num2str(Pop),' ft-lbf/sec of power.'];
disp([string2,string3]);disp(string4);disp(string5);disp(string6),disp(' ');
Dinches=input('What propeller diameter in inches do you want to use? Default=D*  >>>>');
if isempty(Dinches)
	Dactual=Dstar;   % Actual prop diameter is the one for maximum efficiency (ft)
	Jactual=Jstar;  % Actual propeller advance ratio
	Poutactual=Pop;    % Actual propeller power, ft-lbf/sec
	nactual=nstar;  % Actual propeller frequency, Hz
	RPMactual=nactual*60; % Actual propeller RPM
	Etaactual=Etastar; % Actual propeller efficiency
	Ppinactual=rho*nactual^3*Dactual^5*CPstar; % Propeller input power (from shaft) ft-lbf/sec
else
	Dactual=Dinches/12; % Actual diameter is the manually input one converted to ft.
	nj=length(J);
	Pprop(2:nj)=Vop*rho*(Vop./(J(2:nj)*Dactual)).^2*Dactual^4.*CT(2:nj);
	%ifig=ifig+1; figure(ifig); plot(J(2:nj),Pprop(2:nj)); xlabel('J');ylabel('Propeller power ft-lf/sec')
	%string7=['Propeller power as a function of J for the prop of ',num2str(Dinches),' inches.'];
	%title(string7)
	
	%YI = INTERP1(X,Y,XI)
	Jactual=interp1(Pprop(2:nj),J(2:nj),Pop);  % Actual propeller advance ratio
	Poutactual=interp1(J,Pprop,Jactual); % Actual propeller power, ft-lbf/sec
	nactual=Vop/(Jactual*Dactual); % Actual propeller frequency, Hz
	RPMactual=nactual*60; % Actual propeller RPM
	%hold on
	%plot([Jactual,Jactual],[0,Poutactual],'r:',[J(2),Jactual],[Poutactual,Poutactual],':')
	%hold off
	CPoutactual=interp1(J,CP,Jactual);
	Ppinactual=rho*nactual^3*Dactual^5*CPoutactual; % Propeller input power (from shaft) ft-lbf/sec
	CTactual=interp1(J,CT,Jactual);
	Poutactual2=Vop*rho*nactual^2*Dactual^4*CTactual;
	%Etaactual3=Poutactual/Ppinactual;
	%Etaactual2=Poutactual2/Ppinactual;
	Etaactual=interp1(J,CT.*J./CP,Jactual);
end
disp(' '); disp('PROPELLER SUMMARY')
string8=['  A propeller of ',num2str(Dactual*12),' inches'];
string9=['  spinning at ',num2str(RPMactual),' RPM at a translational speed of ',num2str(Vop),' ft/sec'];
string10=['  produces an output power of ',num2str(Poutactual),' ft-lbf/sec,'];
string10a=['  has an efficiency of ',num2str(Etaactual),','];
string11=['  and absorbes an input power of ',num2str(Ppinactual),' ft-lbf/sec or ',num2str(Ppinactual/550),' HP.'];
disp(string8); disp(string9); disp(string10); disp(string10a); disp(string11);
