function [CnpCL_MO, CnpCL,Cnp_w,Cnp_v,Cnp]=Cnp2(AR_w,Lamq_w,Xbar,c_w,B,Clp_w,a,Clp,CL_w,dCnpTheta,Theta,dCnp_add,adelf,delf,b_v,l_v,Z_v,CyB_v)

%for wing
CnpCL_M0=-1/6*(AR_w+6*(AR_w+cos(Lamq_w))*(Xbar/c_w*(tan(Lamq_w)/AR_w+(tan(Lamq_w))^2/12))) / (AR_w+4*cos(Lamq_w)); %eqn 8.10

CnpCL=(AR_w+4*cos(Lamq_w))/(AR_w*B+cos(Lamq_w)) * (AR_w*B+.5*(AR_w*B+cos(Lamq_w))*(tan(Lamq_w))^2)/(AR_w+.5*(AR_w+cos(Lamq_w))*(tan(Lamq_w))^2) * CnpCL_M0;
										%eqn 8.8

Cnp_w=-Clp_w*tan(a)-(-Clp*tan(a)-CnpCL*CL_w)+dCnpTheta*Theta+dCnp_add*adelf*delf;	
										%eqn 8.7
%for vert-tail
Cnp_v=-2/b_v*(l_v*cos(a)+Z_v*sin(a)) * (Z_v*cos(a)-l_v*sin(a))/b_v * CyB_v;
										%eqn 8.11

%add them up
Cnp=Cnp_w+Cnp_v;					%eqn 8.6