Skip to main content
2 of 4
added 587 characters in body
Artes
  • 57.9k
  • 13
  • 159
  • 247

If we are looking for real solutions, there are none.

Reduce[{ x - Exp[-(1 - x - y - z)] == 0, 
         y - Exp[-(1 - x/2 - y - z)] == 0, 
         z - Exp[-(1 - x - 4/5 y - z)] == 0}, {x, y, z}, Reals]
False

Without domain specification (it means we are looking for complex solutions) Reduce has been working for several minutes without yielding results. This is not surprising. Instead we take a pedestrian approach, although it might be automized a bit more.

From the first equation we see that $x=a$ where we introduced $\quad a = \exp(-1+x+y+z)$. Now we have from the second equation we get $y=a \exp(-x/2) $ and from the third equation we get $z=a\exp(-y/5)$. Next we get an equation for x by eliminating $y$ and $z$ i.e. $x=e^{-1+x+x e^{-\frac{x}{2}} + x e^{-\frac{x}{5} e^{-\frac{x}{2}}}}$. We need to restrict appropriately the complex domain since knowing properties of exponential functions we expect there are infinitely many solutions.

Reduce[ x == Exp[-1 + x (1 + Exp[-x/2] + Exp[-(x/5) Exp[-(x/2)]])] && 
        Abs[x] < 2, x]
sx = {ToRules @ % } // N;
 x == Root[{-E^(
   E^(-(#1/2) - 
     1/5 E^(-(#1/2)) #1) (-E^(#1/2 + 1/5 E^(-(#1/2)) #1) + 
      E^(#1/2) #1 + E^(1/5 E^(-(#1/2)) #1) #1 + 
      E^(#1/2 + 1/5 E^(-(#1/2)) #1) #1)) + #1 &, 
0.055803204968110806464128827202522530324263607292 - 
 0.510206373369089749181274412607230169611834407013 I}] || 
x == Root[{-E^(
   E^(-(#1/2) - 
     1/5 E^(-(#1/2)) #1) (-E^(#1/2 + 1/5 E^(-(#1/2)) #1) + 
      E^(#1/2) #1 + E^(1/5 E^(-(#1/2)) #1) #1 + 
      E^(#1/2 + 1/5 E^(-(#1/2)) #1) #1)) + #1 &, 
0.055803204968110806464128827202522530324263607292 + 
 0.510206373369089749181274412607230169611834407013 I}]

The above we've got symbolic representation of exact solutions. For example taking restriction Abs[z]<5 after roughly 3 minutes I got 133 complex solutions for x. The results for y and z might be expressed also in the symbolic form, nontheless for brevity we demonstrate them only in the numerical form

sy = y -> x Exp[-x/2] /. sx

{y -> 0.177717 - 0.466416 I, y -> 0.177717 + 0.466416 I}

z -> x  Exp[-y/5] /. Transpose[{sx // Flatten, sy}]
{z -> 0.0994856 - 0.485233 I, z -> 0.0994856 + 0.485233 I}

Reduce couldn't find solutions without approprate restriction of domains for x, y, z separately. This is crucial point to get any solutions.

Artes
  • 57.9k
  • 13
  • 159
  • 247