Oops, missed pasting a line...
On 11/04/2008, Richard Putman <rcj.putman at googlemail.com> wrote:
>
>
> >
> maybe something along the lines of
>
> powerdisp : true$
y: sum(a[i]*x^i, i, 0, 10);
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
>
>