function CDu=CDu(M,dCDdM)

CDu=M*dCDdM;

return