Bessel differential equation



> From:  esvax
> 
> Maxima can't solve Bessel differential equation:
> 
> (%i2) ode2(x^2*'diff(y,x,2)+x*'diff(y,x)+(x^2*m^2-k^2)*y[x] = 
> 0, y, x);
> %o2) false
> 
> But Mathematica 5.2 solve this quation without problems:

Firstly, there is a syntax error.  Using y[x] is wrong.  Try something
like:

(%i1) ode2( x^2*'diff(y,x,2)+x*'diff(y,x)+(x^2-k^2)*y=0, y, x);
is k an integer?  type y or n. 
y;
(%o1)             y = %k2 bessel_y(k, x) + %k1 bessel_j(k, x)


Your are correct that ode2 cannot solve your ode ...

(%i2) ode2( x^2*'diff(y,x,2)+x*'diff(y,x)+(m^2*x^2-k^2)*y=0, y, x);
(%o2)                                false


... but contrib_ode can

(%i4) load('contrib_ode)$
(%i5) contrib_ode(x^2*'diff(y,x,2)+x*'diff(y,x)+(m^2*x^2-k^2)*y=0, y,
x);
(%o5)        [y = %k2 bessel_y(k, - m x) + %k1 bessel_j(k, - m x)]
(%i6)

This email is confidential and may also be privileged.  If you are not the intended recipient, please notify us immediately and delete this message from your system without first printing or copying it. Any personal data in this email (including any attachments) must be handled in accordance with the Rio Tinto Group Data Protection Policy and all applicable data protection laws.