Vc=0.1:0.1:0.5 ;
V_tip=180 ;
rho=1.225 ;
w=180 ;
sigma=0.047193 ;
Cd_avg=0.01 ;
R=2.235 ;
Rv=80.529 ;
A=pi*R.^2 ;
Vh=sqrt(w./(2*rho*A)) ;
Pi=w*sqrt(w./(2*rho*A)) ;
Po=(Cd_avg*rho*sigma*A.*(V_tip).^3)./8 ;
Ph=Pi+Po ;
C_ph=Ph./((rho*A.*(Rv*R).^3)) ;
Cp=((-Vc./2.*Vh)+sqrt((Vc./2.*Vh).^2+1)).*C_ph ;
Ct=Cp.*V_tip ;
plot(Ct,Vc./Vh,"r",'LineWidth',4)