Scribd
Cargar
18 vistas1 página

Codigo

Este documento presenta tres funciones matemáticas: una función f(c) que calcula un polinomio de cuarto grado, una función f(x) que calcula un polinomio de tercer grado y una función NewtonRaphson que implementa el método de Newton-Raphson para encontrar raíces aproximadas iterativamente.

Cargado por

Bryan Juarez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
18 vistas1 página

Codigo

Este documento presenta tres funciones matemáticas: una función f(c) que calcula un polinomio de cuarto grado, una función f(x) que calcula un polinomio de tercer grado y una función NewtonRaphson que implementa el método de Newton-Raphson para encontrar raíces aproximadas iterativamente.

Cargado por

Bryan Juarez
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

function retval = f(c)

retval = (3*(c)^4)-(4*(c)^3)-(10*(c))-1
endfunction

function retval = f(x)


retval = (12*(x)^3)-(12*(x)^2)-10
endfunction

function retval = NewtonRaphson (x0, imax, ea)

fprintf("Aproximación inicial de x: %d\n", x0);


fprintf("Error requerido: %d\n", ea);
fprintf("Cantidad de iteraciones maxima: %i\n", imax);
x_r = x0;
iteracion = 0;

do
fprintf("Iteracion %i \n", iteracion);
x_r_old = x_r;

x_r = x_r_old - (f(x_r_old) / f_prima(x_r_old));


fprintf(" Nuevo valor de x_r: %d \n", x_r);

iteracion = iteracion + 1;

if (x_r != 0)
es = abs((x_r - x_r_old) / x_r) * 100;
fprintf(" Error aproximado: %d \n", es)
endif

until (es < ea && iteracion <= imax)

retval = x_r;

endfunction

También podría gustarte