'algorithm for real arguments changes at 1.0001
test against analytical result for m=2, n=1'

def test_lqmn_gt1(self):

    x0 = 1.0001
    delta = 2e-05
    for x in ((x0 - delta), (x0 + delta)):
        lq = special.lqmn(2, 1, x)[0][((-1), (-1))]
        expected = (2 / ((x * x) - 1))
        assert_almost_equal(lq, expected)
