function [ClrCL_M, dClrGam,Clr_w,Clr_v,Clr]=Clr2(a,AR_w,CL_w,B,b_w,b_v,l_v,Z_v,Lamq_w,ClrCL_M0,Gam_w,dClrTheta,Theta,dClradelfdelf,adelf,delf,CyB_v)

num=(1+ AR_w*(1-B^2)/(2*b_w*(AR_w*B+2*cos(Lamq_w))) + (AR_w*B+2*cos(Lamq_w)) / (AR_w*B+4*cos(Lamq_w)) * (tan(Lamq_w))^2/8)
den=(1+ (AR_w+2*cos(Lamq_w))/(AR_w+4*cos(Lamq_w)) * (tan(Lamq_w))^2/8);
ClrCL_M=num/den*ClrCL_M0;		%eqn 9.5
dClrGam=(1/12)*(pi*AR_w*sin(Lamq_w))/(AR_w+4*cos(Lamq_w));								%eqn 9.7

Clr_w=CL_w*ClrCL_M + dClrGam*Gam_w + dClrTheta*Theta+dClradelfdelf*adelf*delf;		%eqn 9.4

%for vert-tail
Clr_v=-2/(b_v^2)*(l_v*cos(a)+Z_v*sin(a))*(Z_v*cos(a)-l_v*sin(a)) * CyB_v;			%eqn 9.8

Clr=Clr_w+Clr_v;		%eqn 9.3
