In the docs of special functions [1] there is a warning:
"SciPy’s versions are poorly documented and seem less accurate than the Maxima and PARI versions."
Less accurate? Where does Sage draw the border to unusable code?
I certainly do not expect such incorrect results produced with Sage:
A000262 = lambda n: (-1)^(n-1)*hypergeometric_U(1-n, 2, -1, "scipy")
[round(A000262(n).n(100)) for n in (0..19)]
[1, 1, 3, 13, 73, 501, 4051, 37633, 394353, 4596553, 58941091,
824073141, 12470162233, 202976401213, 3535017524400, 65573803186900,
1290434218670000, 26846616451200000, 588633468315000064,
13564373693600000000]
Peter
[1] http://www.sagemath.org/doc/reference/functions/sage/functions/special.html
[2] https://oeis.org/A000262
