Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Homework Help: Solve system of nonlinear equation

  1. Jul 23, 2010 #1
    Hello guys,

    I am not good in advance mathematics.
    I have system of nonlinear equation and I want to solve it analytically, but I face some difficulties.

    1. The problem statement, all variables and given/known data

    The system is :


    2. Relevant equations

    where: [tex]a,b,c,d,e,[/tex]and [tex]f[/tex] are constants.

    3. The attempt at a solution

    any suggestion and help would be appreciated.
  2. jcsd
  3. Jul 23, 2010 #2


    User Avatar
    Science Advisor

    That looks rather messy, especially for general coefficients like that. I think what I would try would be to solve one of the equations, the first, say, since it is linear in x, for x as a function of y, then put that value back into the second equation:
    [tex](\epsilon y^2+ 2f)x= (a- b)y^2+ c[/tex]
    [tex]x= \frac{(a-b)y^2+ c}{y^2+ 2f}[/tex]

    Replacing x by that in the second equation will result, eventually, in a sixth degree equation for y. There are no general methods of solving such an equation.
  4. Jul 23, 2010 #3
    When you can't (practially) solve it, then aren't you just curious what the answer might be? A CAS is such a wonderful tool in mathematics I believe and extends one's "reach" significantly in my opinion. Now you said you needed an "analytic" solution to it right? And since that can't be done "exactly", why not just try to solve it in Mathematica?

    You could do:

    Code (Text):
    Clear[a, b, c, d, e, f, x, y]
    mysol = Solve[{(a - b) y^2 - e x y^2 - 2 f x == -c,
       2 f x^2 - 2 c x - 2 f y^2 == -d}, {x, y}]
    Surprisingly, Mathematica returns a "solution expression" in mysol. It's the six solutions indicated by Hallsofivy above and is too large for Mathemtica to display on the screen conveniently. In some ways, the expression returned by Mathematica can be considered the "analytic" solution although I assume it contains "root" objects which are then numerically computed when an explicit result is desired. Just for a test, I tried back-substituting the first solution into the left side of the first equation with some simple values for the constants. The result is the expected negative of the value of c which is -4.

    Code (Text):
    N[(a - b) y^2 - e x y^2 - 2 f x //. {a -> 1, b -> 2, e -> 3, c -> 4,
       d -> 5, e -> 6, f -> 7, mysol[[1, 1]], mysol[[1, 2]]}]

    Out[19]= -4.
    Now, if you wanted to study this further, you could then study the particular format Mathematica is using to represent the solutions.

    Work like this offers quite a unique and new perspective into Mathematics education which is beneficial to the student's success in the field in my opinion. They should do it regularly.
    Last edited: Jul 23, 2010
  5. Jul 23, 2010 #4
    Thanks guys,
    I also tried to simplify the system such as:
    [tex] (a-b)(2fx^2 - 2cx + d) - ex(2fx^2 - 2cx + d) - 4f^2x = -2fc[/tex]

    I searched for the cubic formula http://en.wikipedia.org/wiki/Cubic_formula#General_formula_of_roots". I found that there are some restriction (limitation or may say conditions) that gives different solution which I really did not understand it.
    could you explain to me little bit please or any suggestion about solution using cubic formula.
    Last edited by a moderator: Apr 25, 2017
  6. Jul 23, 2010 #5


    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    Since y only appears as y2, you could change variables to u=y2, which leaves you with equations for two conic sections. Perhaps that might help you understand what the solutions are.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook