disp(cdfnor("PQ", -1.5, 0, 1), "P{error|message is 1}")
disp(1-cdfnor("PQ", 2.5, 0, 1), "P{error|message is 0 }")
