MathGroup Archive 1998

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

Search the Archive

x == x^1.2/(2*x^3+y^0.7+4*z^2.5)


Sorry, I lost the original message so I am sending this to the whole
group--

Someone wanted to solve the equations

x == x^1.2/(2*x^3+y^0.7+4*z^2.5)
y == y^0.7/(2*x^0.6+y^2+z^2.2),
z == 0.9*z^1.5/(x^0.7+2*y^0.2+z^1.1)

Here is my reply:

The form of your equations suggests trying something like this:

f[{x_,y_,z_}]:={x^1.2/(2*x^3+y^0.7+4*z^2.5),
       y^0.7/(2*x^0.6+y^2+z^2.2),
       0.9*z^1.5/(x^0.7+2*y^0.2+z^1.1)}

FixedPoint[f,{1,1,1}, 50]

You could raise that 50 on a faster machine than the one I have.

I get a solution of approximately {0.666839, 0.204396, 0}

You could try other starting points and do the same thing.

FixedPoint[f,{0,0,1}, 50] gives {0, 0, 0.838953}.

FixedPoint[f,{0, 1, 0}, 50] gives {0, 1, 0}.

FixedPoint[f,{1, 0, 0}, 50] gives a computational underflow.

If this is from a real situation then your intuition should suggest
where to look for a solution.

There is also a three-dimensional version of Newton's method that you
could try if this doesn't seem to be working.  Let me know if you want
anything on that.

Best wishes,

Will Self



  • Prev by Date: Getting Delayed rules from NSolve
  • Next by Date: RE: Add to menus?
  • Previous by thread: Re: Getting Delayed rules from NSolve
  • Next by thread: (q) motifps