Código en Fortran 90

Código en Fortran 90[editar · editar código]
Programa escrito en Fortran 90 correspondiente al ejemplo f(x) = x3 + 2x2 + 10x - 20
PROGRAM Metodo_Secante
IMPLICIT NONE
REAL (KIND = 8) ::x0, x1, x, f0, f1, f, tol = 1.0E-3
INTEGER (KIND = 1) :: i, ITER_MAX = 25
EXTERNAL f

x0 = 0. ! aproximación inicial 1
x1 = 1. !aproximación inicial 2
f0 = f(x0)
f1 = f(x1)

DO i = 2, ITER_MAX
x = x1 - (x1 - x0)*f1/(f1 - f0)
IF (ABS(x - x1).LT.tol) THENPRINT*, 'La raíz es:', x, 'en iteración No:', i
EXIT
END IF
x0 = x1
x1 = x
f0 = f1f1 = f(x)
END DO

IF (i.GT.ITER_MAX) THEN
PRINT*, 'No se halló raíz: cambiar aprox. iniciales o aumentar ITER_MAX'
END IF

END PROGRAMMetodo_Secante

FUNCTION f(x) RESULT (funcion)
IMPLICIT NONE
REAL (KIND = 8) :: funcion, x
INTRINSIC :: cos

funcion = x**3 + 2.*x**2 + 10.*x - 20. ! funciónEND FUNCTION f
Para compilar en GNU/Linux con compilador de GNU, se escribe en una terminal:
$ gfortran programa.f90 -o programa
$ ./programa

Código en Matlab[editar · editar código]
Programaescrito en Matlab para ejecutar el método de la secante.
% Una implementación del método de la secante para búsqueda de raices en
% funciones continuas dentro de un intervalo.
%
% Por GerardoTinoco Guerrero
%
% Ejemplo:
% Ejecutar las siguientes lineas dentro de la ventana de comandos:
%
% ff = @(x)(x.^2-4)
% x = secante(ff, 2, 5, 0.0001);
%
% Se buscará la raíz de la función(x^2)-4 tomando como puntos iniciales para
% el método de la secante a = 2 y b = 5, con una tolerancia tol = 0.0001.

function xs = secante(fun,a,b,tol)
fprintf('Método de la secante\n\n');
i = 1;...