Ask Your Question
1

A bug in Laguerre-L?

asked 10 years ago

Peter Luschny gravatar image

updated 10 years ago

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!

Preview: (hide)

1 Answer

Sort by » oldest newest most voted
2

answered 10 years ago

kcrisman gravatar image

Great point. Those orthogonal polynomials haven't had much TLC for a long time. I've opened http://trac.sagemath.org/ticket/17192 for this.

Preview: (hide)
link

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 10 years ago

Seen: 397 times

Last updated: Oct 22 '14