function [ClB_wb,ClB_hb,ClB_h,ClB_v,ClB]=ClB2(CL_wb,ClB_CL2,KML,Kf,ClB_CL_A,Gam_w,ClB_Gam,KMG,dClB_Gam,dClB_Zw,theta_w,Lamq_w,dClB_thetatanL4,S_w,b_w,CL_hb,ClB_CL2_h,KML_h,Kf_h,ClB_CL_A_h,Gam_h,ClB_Gam_h,KMG_h,dClB_Gam_h,dClB_Zh,theta_h,Lamq_h,dClB_thetatanL4_h,S_h,b_h,CyB_v,Z_v,a_v,l_v,b_v)


ClB_wb=57.3*(CL_wb*(ClB_CL2*KML*Kf+ClB_CL_A)+Gam_w*(ClB_Gam*KMG+dClB_Gam)+dClB_Zw+theta_w*tan(Lamq_w)*dClB_thetatanL4); %eqn 7.9

ClB_hb=57.3*(CL_hb*(ClB_CL2_h*KML_h*Kf_h+ClB_CL_A_h)+Gam_h*(ClB_Gam_h*KMG_h+dClB_Gam_h)+dClB_Zh+theta_h*tan(Lamq_h)*dClB_thetatanL4_h); %eqn 7.9
ClB_h=ClB_hb*S_h*b_h/(S_w*b_w); %eqn 7.13

ClB_v=CyB_v*(Z_v*cos(a_v)-l_v*sin(a_v))/b_v;	%eqn 7.14

ClB=ClB_wb+ClB_h+ClB_v; %eqn 7.8