All Roots of Two Systems of Algebraic Equations

This Demonstration illustrates how one can use NDSolve's WhenEvent option to determine the roots of a system of nonlinear algebraic equations: and . An attractive feature of the method is that one does not need not need to know at the outset an approximate guess for the roots. All one needs to do is specify the solution domain , .
Consider the two test problems consisting of two-dimensional systems of algebraic equations:
You can vary the parameters and .
This Demonstration finds all the roots of these two problems in the range (for test problem 1) and (for test problem 2).
The first contour plot is the graph of and , and the second shows the level curves of , which is minimized at a solution to the system.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The following describes the solution methodology.
First, define . By construction . If is a root of the system , then and . Thus, a necessary condition is given by .
You can solve the differential-algebraic system of equations:
Eq. (1),
Eq. (2).
In doing so, we are solving the ODE given by Eq. (1) along the curve defined in the x-y plane by Eq. (2). A solution to such a problem is always guaranteed for the particular case when the DAE system has index 1, which is equivalent to . We use the built-in Mathematica function WhenEvent to obtain the values of (i.e., the candidate roots) that verify the necessary condition . For each value of , we check if the system of equations is indeed verified by calculating . If , that particular candidate solution is simply rejected.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-Step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2018 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+