Mathematica Solving systems of nonlinear ordinary equations in Mathematica.

AI Thread Summary
The discussion focuses on solving a complex system of five nonlinear equations of degree seven in Mathematica, which is taking an excessively long time to compute. Suggestions include simplifying the equations by substituting common expressions, such as replacing (e-q) with a new variable to reduce complexity. Participants highlight the importance of identifying and factoring out large common subexpressions to improve processing time. It is noted that attempting to solve such complicated multivariate polynomials may not yield results in a reasonable timeframe. Overall, simplification and strategic substitutions are key to enhancing the efficiency of the solution process.
fab333
Messages
3
Reaction score
0
I'm trying to solve a system of five nonlinear equations (degree 7) in Mathematica:

Code:
Solve[{q (a^2 + 
       b^2 (3 q^2 + r^2 + 8 t^2)) ((e - q)^2 - R^2 - 4 T^2)^2 + 
    q ((e - q)^2 + R^2 + 4 T^2) - 2 (e - q) (r R + 4 t T) == 0, 
  R (a^2 + b^2 (3 R^2 + (e - q)^2 + 8 T^2)) ((e - q)^2 - R^2 - 
        4 T^2)^2 + R (R (R - r) - (e - 3 q) (e - q) + 4 T (T - 2 T)) -
     r ((e - q)^2 + 4 T^2) == 0, 
  r (a^2 + b^2 (3 r^2 + q^2 + 8 t^2)) ((e - q)^2 - R^2 - 4 T^2)^2 - 
    r (R^2 + (e - q)^2 + 4 T ^2) + 2 R (q (e - q) - 4 t T) == 0, 
  T (a^2 + 2 b^2 (2 T^2 + (e - q)^2 + R^2)) ((e - q)^2 - R^2 - 
        4 T^2)^2 + T (4 T (T - t) - (e - 3 q) (e - q) + R (R - 2 r)) -
     t ((e - q)^2 + R^2) == 0, 
  t (a^2 + 2 b^2 (2 t^2 + q^2 + r^2)) ((e - q)^2 - R^2 - 4 T^2)^2 - 
    t (4 T^2 + (e - q)^2 - R ^2) + 2 T (q (e - q) - rR) == 0
  }, {q, r, t, R, T}]

but the program is running for an hour...
there is hope that it get to the end?

How do I have to change it to hope to arrive at the conclusion?

Thanks
 
Physics news on Phys.org
If it isn't just a simple scrape and paste error then this worries me

) - rR) == 0

Perhaps you want a space between r and R there.

Sometimes when you have very complicated expressions you can make a simplification that can speed things up. For example, you have dozens of (e-q) in your equations. If you could substitute emq for (e-q) everywhere, patch the few isolated remaining q with (e-emq) and see whether this resulted in enough simplification then it might help speed things up.

I would also look for large common subexpressions and see whether you might be able to find a way to substantially simplify this before throwing it at Solve.

In general I would expect throwing 5 seventh order complicated multivariate polynomials at Solve to never finish in an acceptable time or with an acceptable result.
 
Assumptions may help as well.
 

Similar threads

Replies
2
Views
1K
Replies
2
Views
2K
Replies
1
Views
2K
Replies
19
Views
2K
Replies
2
Views
835
Back
Top