Hi, I have a problem with bessel functions. First I load(orthopoly) Then I tried Legendre_p(3,x) and got the expression. Then I tried Bessel_j(3,x) and the answer is Bessel_j(3,x) but bessel_j(0,0) gives 1 bessel_j(1,0) gives 0 Can someone help me with this problem; Thanks