A work around is to call odelin directly
(%i5) load(odelin)$
(%i6) odelin(x^2*'diff(y,x,2)+x*'diff(y,x)+(a^2*x^2-b^2)*y=0,y,x);
(%o6) {bessel_j(b,-a*x),bessel_y(b,-a*x)}
I think ode2 is supposed to call odelin automatically; apparently something
in this
process is broken. Maybe you could file a bug report.
Barton
-----maxima-bounces at math.utexas.edu wrote: -----
>I?am?trying?to?solve?the?following?diff?eq:
>
>x^2*diff(y,x,2)+x*diff(y,x)+(a^2*x^2-b^2)*y=0
>
>which?should?lead?to?the?following?solution:
>
>y=c1*bessel_j(b,?a*x)?+?c2*bessel_y(b,a*x)
>
>If?I?use?ode2?to?solve,?I?get?the?message?"false",?meaning?it?couldn't
>solve?it,?and?if?I?use?odelin?or?contrib_ode,?I?get?the?message
>
>EVAL:?too?many?arguments?given?to?SRATSIMP:?(SRATSIMP?P?X)
>Automatically?continuing.
>To?reenable?the?Lisp?debugger?set?*debugger-hook*?to?nil.
>
>However,?if?I?ask?ode2?to?solve:
>
>x^2*diff(y,x,2)+x*diff(y,x)+(x^2-1)*y=0
>
>I?get:
>
>y=c1*bessel_j(1,x)?+?c2*bessel_y(1,x)
>
>which?is,?of?course,?the?correct?solution.??Clearly,?I?am?not?doing
>something?correctly?with?the?constants?a?and?b.??Might?there?be?someone
>who?could?provide?me?with?some?insight???Thank?you?all?for?your?help!
>
>
>Cheers,
>Carl
>_______________________________________________
>Maxima?mailing?list
>Maxima at math.utexas.edu
>http://www.math.utexas.edu/mailman/listinfo/maxima