MathGroup Archive 1996

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

Search the Archive

Re: differential equation

In article <51c4mm$9no at>,
   Ralph Gensheimer <ralphg at> wrote:
->dear group,
->i have the folowing problem:
->in an ordinary differential equation system,
->(solution functions are  y1(r), y2(r), y3(r) ) is an integration in 
another variable x.
->and in this integral there is also the solution function y3(r).
->r and x are independent.
->the problem has the following structure:
->NDSolve[{y1'[r] == y2[r],
->         y2'[r] == -(2/r)*y2[r]+
->         y1[r]*y3'[r] == y2[r],
->                                 y1[0.001]==-1000,
->                                 y2[0.001]==0,
->                                 y3[0.001]==0 },{y1,y2,y3},{r,0.001,0.2}]
->can you give me ideas to solve this problem ?

Maybe a successive approximation approach?  Replace the instance of y3[r] 
in the integrand (only) with a new function y4[r].  Initially define 
y4[r_]:=0.  Run NDSolve (still solving only for y1, y2 and y3).  Redefine 
y4[r_]:= y3[r].  Iterate ad nauseum, and hope it converges?

Caveat:  I don't do differential equations, so I have no idea if this will 
in fact converge.  It might just be a new way to waste cpu cycles.

-- Paul

* Paul A. Rubin                                  Phone: (517) 432-3509   *
* Department of Management                       Fax:   (517) 432-1111   *
* Eli Broad Graduate School of Management        Net:   RUBIN at MSU.EDU    *
* Michigan State University                                              *
* East Lansing, MI  48824-1122  (USA)                                    *
Mathematicians are like Frenchmen:  whenever you say something to them,
they translate it into their own language, and at once it is something
entirely different.                                    J. W. v. GOETHE


  • Prev by Date: Re: functional code
  • Next by Date: Re: Display a List of Graphics
  • Previous by thread: differential equation
  • Next by thread: equ. of 'ls' in 'Context' ?