Group Abstract Group Abstract

Message Boards Message Boards

0
|
4.6K Views
|
18 Replies
|
8 Total Likes
View groups...
Share
Share this post:

System of 2 nonlinear equations: no result

Posted 10 years ago
POSTED BY: max ma
18 Replies
POSTED BY: Daniel Lichtblau

Cheating or not MMA has problems with this equation. I hope that development team fix this issue, will improve in a future release.

Thanks for Reply.

POSTED BY: Mariusz Iwaniuk
In[1]:= eq1 = 
  FullSimplify[(\[Lambda] \[Alpha] + \[Lambda] \[Mu] y) (p - \[Theta] \
x^2 - (1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 \[Theta] x) == 0];
eq2 = FullSimplify[(\[Lambda] (1 - \[Alpha]) + \[Lambda] \[Mu] x) (p \
- \[Theta] x^2 - (1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 (1 - \[Theta]) y) == 0];

In[3]:= AbsoluteTiming[Solve[Join[eq1, eq2], {x, y}];]

Out[3]= {0.627934, Null}

Are you running out of memory? My PC has 16 GBytes of RAM.

POSTED BY: Frank Kampas
Posted 10 years ago
POSTED BY: max ma

Mr. Frank Kampas you made ??a mistake.

We want to solve these equations:

eq1 = FullSimplify[1/2 ? y^(1/2) x^(-(1/2))(p - 1/2 x^2 - 1/2 y^2) + (a - k p + ? (x y)^(1/2)) (-x) == 0]
eq2 = FullSimplify[1/2 ? x^(1/2) y^(-(1/2))(p - 1/2 x^2 - 1/2 y^2) + (a - k p + ? (x y)^(1/2)) (-y) == 0]

not those:

eq1 = FullSimplify[(? ? + ? ? y) (p - ? x^2 - (1 - ?) y^2 - b - d) + (a - k p + ? (? x + (1 - ?) y + ? ? y)) (-2 ? x) == 0]
eq2 = FullSimplify[(? (1 - ?) + ? ? x) (p - ? x^2 - (1 - ?) y^2 - b - d) + (a - k p + ? (? x + (1 - ?) y + ? ? y)) (-2 (1 - ?) y) == 0]

My Laptop has 4 GBytes of RAM.With these 2 equations Kernel level of consumption of RAM is only 55MB and did not increase for the next hours calculating.

See file attachments:

Attachments:
POSTED BY: Mariusz Iwaniuk
POSTED BY: Frank Kampas
Posted 10 years ago
POSTED BY: max ma

Try with this ?

eq1 = FullSimplify[1/2 ? y^(1/2) x^(-(1/2))(p - 1/2 x^2 - 1/2 y^2) + (a - k p + ? (x y)^(1/2)) (-x) == 0]
eq2 = FullSimplify[1/2 ? x^(1/2) y^(-(1/2))(p - 1/2 x^2 - 1/2 y^2) + (a - k p + ? (x y)^(1/2)) (-y) == 0]
Solve[Join[eq1, eq2], {x, y}]

I am waiting and nothing.

Maples solutions after few seconds. enter image description here

POSTED BY: Mariusz Iwaniuk
Posted 10 years ago

Yes, you are right, transferring Windows version is not helpful. Let me try Maples. Thanks a lot.

POSTED BY: max ma
POSTED BY: Mariusz Iwaniuk
Posted 10 years ago

Thanks for your help. I tried this method, but still no result, could you guide me your code? Thanks a lot.

POSTED BY: max ma

Maybe you can start MMA with a fresh kernel.I'm using Mathematica 10.2 and works fine.

Attachments:
POSTED BY: Mariusz Iwaniuk
Posted 10 years ago

Really appreciate your help. My MMA version is 10.2 too, but it is Mac version, not Windows version. I reduced the symbolic coefficients as much as possible, and still not work, should I transfer Windows version also?

Attachments:
POSTED BY: max ma

Mac or Windows, both versions are the same.MMA has difficulty with this equations.After several hours of computation I'm not saw any result.I do not know why it has such problems. Only numericall is a quick solution with MMA. Maple 2015 solves these equations very quickly.

Attachments:
POSTED BY: Mariusz Iwaniuk

Your system is of high degree with symbolic coefficients. A full symbolic solution may be too complicated. If you lower your expectations, you get quick explicit solutions for specific numeric values of the parameters, for example

With[{\[Lambda] = 1, \[Alpha] = -1, \[Mu] = 1, p = 1, \[Theta] = 1, 
  b = 1, d = 1, a = 1, k = 1},
 Solve[{(\[Lambda] \[Alpha] + \[Lambda] \[Mu] y) (p - \[Theta] x^2 - \
(1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 \[Theta] x) == 
    0, (\[Lambda] (1 - \[Alpha]) + \[Lambda] \[Mu] x) (p - \[Theta] \
x^2 - (1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 (1 - \[Theta]) y) == 0}, {x, y}]]

You can go interactive too:

Manipulate[
 ContourPlot[{(\[Lambda] \[Alpha] + \[Lambda] \[Mu] y) (p - \[Theta] \
x^2 - (1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 \[Theta] x) == 
    0, (\[Lambda] (1 - \[Alpha]) + \[Lambda] \[Mu] x) (p - \[Theta] \
x^2 - (1 - \[Theta]) y^2 - b - d) + (a - 
        k p + \[Lambda] (\[Alpha] x + (1 - \[Alpha]) y + \[Lambda] \
\[Mu] y)) (-2 (1 - \[Theta]) y) == 0}, {x, -1, 1}, {y, -1, 1}],
 {\[Lambda], -1, 1}, {\[Alpha], -1, 1}, {\[Mu], -1, 1}, {p, -1, 
  1}, {\[Theta], -1, 1}, {b, -1, 1}, {d, -1, 1}, {a, -1, 1}, {k, -1, 
  1}]
POSTED BY: Gianluca Gorni
Posted 10 years ago

Thanks for your help. Trying to reduce the variables numbers, but can not reduce all...but know the reason why it can not run. Thanks.

POSTED BY: max ma
Posted 10 years ago

Added code and the .nb file. Thanks.

Solve[{(? ? + ? ? y) (p - ? x^2 - (1 - ?) y^2 - b - d) + 
(a - k p + ? (? x + (1 - ?) y + ? ? y)) (-2 ? x) == 0, 
(? (1 - ?) + ? ? x) (p - ? x^2 - (1 - ?) y^2 - b - d) + 
(a - k p + ? (? x + (1 - ?) y + ? ? y)) (-2 (1 - ?) y) == 0}, {x, y}]
Attachments:
POSTED BY: max ma

Please post CODE not the image or attach the notebook.

POSTED BY: EDITORIAL BOARD
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard