> From: Radovan Omorjan
>
> Billinghurst, David (RTATECH) wrote:
> >> From: Radovan Omorjan
> >> Hello,
> >>
> >
> > %i1) load('contrib_ode)$
> >
> > (%i2) eq2:x^2*'diff(y,x,2)+x*'diff(y,x)-(x^2+4)*y=0;
> > 2
> > 2 d y dy 2
> > (%o2) x --- + x -- + (- x - 4) y = 0
> > 2 dx
> > dx
> > (%i3) contrib_ode(eq2,y,x);
> > (%o3) [y = bessel_y(2, - %i x) %k2 + bessel_j(2, - %i x) %k1]
> >
> Thanks for the answer, I checked it out.
> By the way, why do not use bessel_k and bessel_i - they are
> also defined in Maxima.
> I guess you find it more suitable this way.
>
For second order ODEs, contrib_ode just tries ode2 then odelin. You
would have
to ask Barton Willis why he wrote it that way. I personally think it is
a fine
answer, but maxima could simplify the Bessel functions of complex
arguments to
modified Bessel functions. There has been significant progress in
complex
simplification recently, and this is now probably trivial to add.
David
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.