function phi = MC(a)
    phi = max(0,min((1+a)./2,min(2,2*a)));
end