REGLA FALSA
clc
clear
clear all
format long
ingresefuncion=input('ingrese la funcion: ', 's'); %Aqui el usuario debe ingresar la funcion para calcularle
la raiz
x0=input('ingrese el primer valor inicial: '); %Se debe ingresar la primera aproximacion
x1=input('ingrese el segundo valor inicial: '); %Se debe ingresar la segunda aproximacion
tol=input('ingrese la tolerancia deseada: '); %Se debe ingresar la tolerancia deseada por el usuario
iter=input('ingrese el numero de iteraciones: '); %Se debe ingresar el maximo numero de iteraciones a
realizar
funcion=inline(ingresefuncion); %Aqui se convierte la funcion en una serie de caracteres para poder evaluar
valores
fx0=funcion(x0); %Aqui se evalua la funcion en x0
fx1=funcion(x1); %Aqui se evalua la funcion en x1
if fx1==0 %Si al evaluar la funcion en x1 da 0 entonces indica que la raiz es x1
fprintf('La raiz es : %g \n',x1);
else
if fx0==0
fprintf('La raiz es : %g \n',x0); %Si al evaluar la funcion en x0 da 0 entonces indica que la
raiz es x0
else
if fx0*fx1<0 %Si la funcion en x0 y x1 tiene cambio de signo entonces se utiliza este intervalo
x2= x0-((fx0*(x1-x0))/(fx1-fx0)); %Encuentra el valor del punto medio
fx2=funcion(x2); %Evalua la funcion en el punto medio encontrado
i=1; %Se inicia el contador
error=tol+1; %Se inicia el error
while error>tol && fx2~=0 && i<=iter %Mientras el error sea mayor a la tolerancia, fxm
diferente de 0 y la iteracion actual sea menor o igual a iter, evalue las siguientes condiciones
if fx0*fx2<0 %Si la funcion evaluada en a y xm cambia de signo entonces hacer limite
superior xm y evalue la funcion
x1=x2; %Se iguala x1 con x2
fx1=fx2; %Como las variables son iguales. Las respectivas funciones también lo
seran
