Una ecuación lineal es aquella en la que sus incógnitas están solas (con su coeficiente) y tienen exponente 1. La que no cumpla estas condiciones es una ecuación no lineal.

Método de Gauss-Seidel

Es necesario despejar cada variable manualmente y codificarla como una función. El resto del algoritmo es básicamente igual que en sistemas de ecuaciones lineales. Este método no es muy apropiado si la función es relativamente compleja.

Método de Newton-Rapson

Partiendo del desarrollo en serie de Taylor, este método reduce el sistema de ecuaciones no lineales a uno con ecuaciones lineales en las que aparecen derivadas de las funciones originales como nuevos coeficientes y las funciones originales (con signo negativo) como términos independientes. Será necesario codificar las funciones (ecuaciones) originales y sus derivadas (que se calcularán a mano) en funciones individuales.

El algoritmo resuelve el nuevo sistema de ecuaciones iterativamente, de forma que en cada iteración se aproxima más a la solución.

La solucion del algoritmo queda en el vector xvect.

El procedimiento RellenarMat(mat, xvect) será:

Se produce aquí un pequeño cambio de variable. Las incógnitas del nuevo sistema de ecuaciones lineales ya no son las iniciales (supongamos x). De forma que:

O sea, hi es la diferencia entre la variable xi de la iteración anterior y la actual.