Fortunately, the MATLAB fsolve function is easy to use for solving algebraic equations. For a simplified presentation, we use the form
obtained from Equation (2.66) with a fixed p and u. The most commonly used numerical techniques are related to Newton-Raphson iteration. The “guess” for iteration k + 1 is determined from the value at iteration k, using
where f[x(k)] is the vector of function evaluations at iteration k, and J(k) is the Jacobian matrix
The ij element ...