[Date Index]
[Thread Index]
[Author Index]
Re: Methods of solving nonlinear equation
*To*: mathgroup at smc.vnet.net
*Subject*: [mg2489] Re: Methods of solving nonlinear equation
*From*: "Brian J. Albright" <albright at physics.ucla.edu>
*Date*: Wed, 15 Nov 1995 01:57:06 -0500
*Organization*: UCLA Department of Physics and Astronomy
rdm at hpbs2500.boi.hp.com (Bob Morrison) wrote:
>
>I'm a non-expert with Mathematica and mathematics, so this is probably
>kind of a stupid question: when I tried to solve a common nonlinear
>equation D^2 y[t] + k/(y[t]^2) == 0, I seem to be unable to get solutions.
>On paper, there is a degenerate solution at (a +/-b t)^(2/3), and
>also on paper, you can derive a polynomial solution (a0 + a1 t +
> a2 t^2 + ....) which shows two degrees of freedom (a0 and a1 are
>independently unconstrained). However, the polynomial solution's
>region of convergence is too small to be very useful to me, and the
>boundary conditions of my problem disallow the degenerate solution.
> I tried using Mathematica to get solutions using DSolve and NDSolve,
>but am too inexperienced to get anything useful. It turns out if you
>assume that y is of form (f(t))^(2/3) you get another nonlinear equation
>for f(t) which appears to be just as unsolvable. I tried representing
>y in Exp[f(t)] form, and failed here. I tried separation
>of variables and didn't get anywhere (to be honest, I can't remember
>why, it's been a little while since I tried it). Can anybody
>provide a reference or maybe a little hint how to do this? Ideally,
>I would like a closed form solution, of course, but a limiting function
>would also be useful. I created a difference equation which correctly
>shows a coarse version of the solution, but in the region of interest,
>the dt portion is too large to correctly solve the problem. I think
>the solutions will involve an elliptic integral, but I'm just guessing
>here. I freely admit this is not my area of expertise and may have
>overlooked the obvious solution, but I don't see it...
>
>Thanks, Bob Morrison, rdm at hpbs670.boi.hp.com
>
Hi Bob.
I think the standard trick for doing something like y'' + k / y^2 = 0
is to multiply both sides by y' and then integrate w.r.t. x. This gives you
something like y' = [ c + 2 k / y ]^(1/2). (c is just an integration
constant.)
Rearranging, you get x = x0 + Integrate[ 1 / Sqrt[ c + 2 k / y ], y],
which you can either do by hand or have mma do for you. This will prob.
give you some kind of trancendental relationship betw. x and y which you
can look at in whatever limit you desire. The solution will be parametrized
by the two integration constants, x0 and c.
Hope this helps.
-Brian
--
Brian J. Albright |
Department of Physics and Astronomy, UCLA | To err is human...
albright at physics.ucla.edu | to err really big
http://bohm.physics.ucla.edu/~albright | is government.
Prev by Date:
**Re: 2-d graphs to 3-d**
Next by Date:
**Re: Request for MMA benchmarks & platform comparisons**
Previous by thread:
**Re: Methods of solving nonlinear equation**
Next by thread:
**Re: Imaginary ...**
| |