A Graphically Enhanced Method for Computing Real Roots of Nonlinear Functions

Requires a Wolfram Notebook System

Interact on desktop, mobile and cloud with the free Wolfram CDF Player or other Wolfram Language products.

Requires a Wolfram Notebook System

Edit on desktop, mobile and cloud with any Wolfram Language product.

Mathematica's built-in function FindRoot is the preferred method for finding a zero of a nonlinear function . However, determining all the zeros of in a specified interval using FindRoot can be challenging, as this requires having good initial guesses for all the roots. In some cases, the initial guesses can be found by examining the asymptotic structure of as , but this is not always possible.

[more]

This Demonstration shows how the computation of the zeros for a nonlinear equation can readily be accomplished by augmenting the calculation with a graphical technique to obtain the initial guesses in a specified interval . The functions considered in this Demonstration are (the Bessel function of the first kind) or a linear combinations of Bessel functions such as or . Another possible problem involves , in a numerical solution of the Lorenz equation. The general idea of calculating approximate roots of a function using graphical techniques is discussed in [1].

You can specify the order of the Bessel functions using a slider. For the Lorenz function, you can specify the parameters in the Lorenz equations. For certain values of the Lorenz parameters, the solution is chaotic, and thus the zeros of are not predictable. The interval for the roots can be specified using the sliders for and , where . Clicking the "function plot" button shows a plot of the selected function with the calculated roots superimposed on the plot as magenta dots. You can view the numerical values of the zeros by clicking the "first few roots" button.

[less]

Contributed by: Housam Binous and Brian G. Higgins (February 2012)
Open content licensed under CC BY-NC-SA


Snapshots


Details

The Lorenz equations in chaotic dynamics are

,

,

,

.

The zeros of the function satisfy .

Reference

[1] S. Wagon, Mathematica® in Action: Problem Solving through Visualization and Computation, 3rd ed., Berlin: Springer-Verlag, 2010.



Feedback (field required)
Email (field required) Name
Occupation Organization
Note: Your message & contact information may be shared with the author of any specific Demonstration for which you give feedback.
Send