6.3. Cálculo numérico

6.3.1. Resolución de ecuaciones en una variable (Newton)

Podemos utilizar la función Newton para resolver ecuaciones en una variable. Se utilizará el método de Newton:

In> Newton(x^3-x^2+x-3,x,1,0.001);
Out> 1.5747430742;
In>

Los argumentos indican lo siguiente:

• La ecuación de la que se quiere encontrar una solución.

• La variable, sólo puede haber una.

• Punto inicial.

• Precisión.

6.3.2. Resolución sistemas de ecuaciones

Podemos utilizar Solve (Sección 6.2.2) para resolver sistemas de ecuaciones:

In> Solve({x+y+z==1,x-y-z==2,x-y+z==3},{x,y,z});
Out> {{3/2,-1,1/2}};
In>
	  

Para formar las ecuaciones utilizaremos el operador == y no =.

6.3.3. Cálculo de polinomios interpoladores

Podemos calcular el polinomio interpolador utilizando LagrangeInterpolant. Esta función utiliza el método de Lagrange para calcular el polinomio interpolador:

In> f(x):=LagrangeInterpolant({1,2,3},{5,8,-5},x);
Out> True;
In> f(1);
Out> 5;
In> Simplify(f(x));
Out> -8*x^2+27*x-14;
In>