changes to:

@XFAIL
def test_P36():
    M=Matrix([
        [10, 7],
        [7, 17]])
    assert sqrt(M) == Matrix([[3, 1], [1, 4]])