function [Cmadot_h,Cmadot]=Cmadot2(Cmadot_w,CLa_h,eta_h,V_h,Xh,c_h,deda)


Cmadot_h=-2*CLa_h*eta_h*V_h*Xh/c_h*deda; %eqn6.5
Cmadot = Cmadot_w+Cmadot_h;


return