M_MPC = M-N*AA(1,:);
	N_MPC = sum(BB(1,:))*N;
	C_MPC = C;
	D_MPC = sum(BB(1,:))*D;  % a D = 0