On 11/04/2008, Harald Geyer <harald at lefant.net> wrote:
>
> Hi!
>
>
> > Hi there: I am trying to solve Airy's equation by power series
> > I found a manual that says that this can done in DOE-macsyma with
> > load(series)$
>
> <snip>
> I don't know how to do this with current maxima, but it's a nice
> problem and probably not too difficult. Anyone?
>
>
maybe something along the lines of
powerdisp : true$
eq: diff(y,x,2)-x*y$
coeffs: makelist(ratcoef(eq,x,i),i,0,9);
ans: ev(linsolve(reverse(coeffs), makelist(a[i],i,0,10)),
linsolve_params=false);
subst(ans, y);
subst(a[0]/6,a[3], %);
I've not checked to see if that is the right answer btw, but I hope it shows
the general idea.
Cheers,
Richard