laguerre(-1,1) returns 0 in Sage but returns exp(1) in Maple and Mathematica.
The docs say: "Return the Laguerre polynomial for integers n>−1."
Is this expected behaviour to return an innocent value instead of giving the user a warning?
Revision history [back]
 | 1 | initial version |
A bug in Laguerre-L?
| | 2 | No.2 Revision |
A bug in Laguerre-L?
laguerre(-1,1) returns 0 in Sage but returns exp(1) in Maple and Mathematica.
The docs say: "Return the Laguerre polynomial for integers n>−1."
Is this expected behaviour to return an innocent value instead of giving the user a warning?
Edit: I give an example how this behaviour can bite: Look at the Inverse Lah transform of 1,2,3,... which starts 1, -2, -1, 2, 17, 94, 487, ... This can be computed by
factorial(n)*(laguerre(n,1)-2*laguerre(n-1,1))
This expression gives for n=0 the expected value '1' and thus hides the fact that the formula is true only for n>0!
