Nonlinear equations

Introduction

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Residual Form

Typically when solving a nonlinear equation we solve for its roots - i.e. where the equation equals zero. In some cases we want to solve an equation for where $f(x)=\alpha$. In such a case, we must rewrite the equation in residual form:

Original Equation

Residual Form

$f(x)=\alpha$

$r(x)=f(x)-\alpha$

For example, if we wanted to solve $y=x^2\sin(x)$ for the value of $x$ where $y=3$, we would solve for the roots of the function $r(x)=x^2\sin(x)-3$.

Solving a Single Nonlinear Equation

There are two general classes of techniques for solving nonlinear equations:

• Closed domain methods - These methods look for a root inside a specified interval by iteratively shrinking the interval containing the root. These methods are typically more robust than open domain methods, but are somewhat slower to converge, and require that you be able to bracket the root.

• Open domain methods - These methods look for a root near a specified initial guess, but are not constrained in the region that they search for a root. These methods typically converge faster than closed domain methods, if the initial guess is "close enough" to the root.

Closed Domain Methods

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

==== Bisection ====

Open Domain Methods

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Solving a System of Nonlinear Equations

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Solution Approaches: Single Nonlinear Equation

Excel

To solve a nonlinear equation in Excel, we have to options:

• Goal Seek is a simple way to solve a single nonlinear equation.

• Solver is a more powerful way to solve a nonlinear equation. It can also perform minimization, maximization, and can solve systems of nonlinear equations as well.

Goal Seek

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Solver

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Matlab

To solve a single nonlinear equation in Matlab, we use the fzero function. If, however, we are solving for the roots of a polynomial, we can use the roots function. This will solve for all of the polynomial roots (both real and imaginary).

ROOTS

The roots function can be used to obtain all of the roots (both real and imaginary) of a polynomial

  xroots = roots( coefs );

• coefs are the polynomial coefficients, in descending order.
• xroots is a vector containing all of the polynomial roots.

For example, consider the quadratic equation $y=3x^2-2x+1$. From the quadratic formula, we can identify the roots as

$x = \frac{1 \pm i\sqrt{2}}{3}$

where $i\equiv\sqrt{-1}$ is the imaginary number.

Using the roots function we have

 x = roots( [3 -2 1] );

which gives

 x =
  0.3333 + 0.4714i
  0.3333 - 0.4714i

This is equivalent to the answer above obtained via the quadratic formula.

FZERO

In Matlab, fzero is the most flexible way to find the roots of a nonlinear equation. Its syntax is slightly different depending on the type of function we are using:

• For a function stored in an m-file named myFunction.m we use

x = fzero( 'myFunction', xguess );

• For an anonymous function named fun>/tt>

x = fzero( fun, xguess );

• For a built in function (like <tt>sin, exp, etc.)

x = fzero( @sin, xguess );

Note that fzero can only find real roots. Therefore, if we tried to use it on the quadratic function $y=3x^2-2x+1$, it would fail. To demonstrate its usage, let's solve for the roots of the function

$y=3x^2-2x-1$

From the quadratic formula, or by factoring this equation, we know that its roots should be $x=-1/3$ and $x=1$.

To use fzero to solve this, we could use an anonymous function:

f = @(x)(3*x.^2-2*x-1);
x1 = fzero( f, -2 ); % start looking for the root at -2.0
x2 = fzero( f,  2 ); % start looking for the root at 2.0

This results in x1=-0.33333 and x2=1. Note that here we found the two roots by using different initial guesses. Of course it helps to know where the roots are so that you can supply decent initial guesses!

To solve this using a m-file, first create the function you will use:

myQuadratic.m

function y = myQuadratic(x) y = 3*x.^2 - 2*x -1

Save this in a file myQuadratic.m and then use fzero:

x1 = fzero( 'myQuadratic', -2.0 );  % start looking for the root at -2.0
x2 = fzero( 'myQuadratic',  2.0 );  % start looking for the root at 2.0

FMINSEARCH

Occasionally we want to maximize or minimize a function. Matlab provides a tool called fminsearch to do this. It searches for the minimum of a function near a starting guess. As with fzero, you use this in different ways depending on what kind of function you are minimizing:

 x = fminsearch( 'myFun', xo );  % when using a function in an m-file
 x = fminsearch( myFun,   xo );  % when using an anonymous function
 x = fminsearch( @myFun,  xo );  % when using a built-in function

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Solution Approaches: System of Nonlinear Equations

This section is a stub and needs to be expanded. If you can provide information or finish this section you're welcome to do so and then remove this message afterwards.

Excel

Solver...

Matlab

fsolve