def test_prbnormndpc():

    rho2 = np.random.rand(2)
    a2 = np.zeros(2)
    b2 = np.repeat(np.inf, 2)