MathGroup Archive 1999

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: O.D.E in Power Series


Hi Alan,

if you want a general infinite power series you should try
the SpecialFunction.m package from MathSource.
If a solution to the first order terms is sufficient
for you than

deqn=y''[x]-2(x+3)y'[x]-3y[x]==0;
order=10;
Solve[LogicalExpand[
    deqn /.Thread[ {#,D[#,x],D[#,{x,2}]} & /@ (
            y[x]->Sum[a[i]*x^i,{i,0,order}]+O[x]^(order+1))]],
  Table[a[i],{i,2,order}]]

will do it.

Regards
  Jens

ALAN S BLOOM wrote:
> 
> I am trying to solve this equation below:
> y''-2(x+3)y'-3y=0 with the DSolve command.
> 
> DSolve [y''[x]-2(x+3)y'[x]-3y[x]==0,y[x],x]
> 
> How can I make it work?
> 
> Alan


  • Prev by Date: Re: CounterBox to different notebook
  • Next by Date: newbie Mathematica question (astrophysics problem)
  • Previous by thread: Re: O.D.E in Power Series
  • Next by thread: Re: O.D.E in Power Series