Scribd
Cargar
294 vistas2 páginas

Método de Falsa Posición en MATLAB

Este documento describe el método de falsa posición para encontrar raíces de una función mediante el uso de MATLAB. Proporciona instrucciones para ingresar la función, los intervalos iniciales a y b, y el error de tolerancia. Luego, itera calculando un nuevo punto xi dentro del intervalo y evaluando la función allí hasta que la diferencia entre a y b sea menor que el error de tolerancia especificado, encontrando así la raíz de la función dentro de ese intervalo.

Cargado por

Gisela Mariana Jimenez
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
294 vistas2 páginas

Método de Falsa Posición en MATLAB

Este documento describe el método de falsa posición para encontrar raíces de una función mediante el uso de MATLAB. Proporciona instrucciones para ingresar la función, los intervalos iniciales a y b, y el error de tolerancia. Luego, itera calculando un nuevo punto xi dentro del intervalo y evaluando la función allí hasta que la diferencia entre a y b sea menor que el error de tolerancia especificado, encontrando así la raíz de la función dentro de ese intervalo.

Cargado por

Gisela Mariana Jimenez
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

%METODO DE FALSA POSICION- METODOS NUMERICOS MATLAB

clc;
syms x
Fx=input('Ingrese la funcion: f(x): ','s');
a=input ('Ingrese el intérvalo a: ');
b=input('Ingrese el intérvalo b: ');
Tol=input('Ingrese el error e: ');
x= a;
Fa=eval(Fx);
x=b;
Fb=eval(Fx);
tol=100;
%mostramos los datos de cabecera
fprintf('\n %6s %7s %8s %10s %8s %8s %8s
\n','A','B','Xi','F(a)','F(b)','FXi','|b-a|');

while abs(b-a)>Tol
%Fórmula de la Falsa posición
xi=((a*Fb)-(b*Fa))/(Fb-Fa);
%evaluamos la funcion en este nuevo punto
x=xi;
Fxi=eval(Fx);
fprintf('\n %8.4f %8.4f %8.4f %8.4f %8.4f %8.4f
%8.4f,\n',a,b,xi,Fa,Fb,Fxi,abs(b-a));
if abs(Fb)<Tol
break;
else
if Fa*Fb<=0
b=xi;
Fb=Fxi;
else
a=xi;
Fa=Fxi;
end
end
end
fprintf('\nLa raiz es: %.4f \n',xi);
ezplot(Fx)%graficar la función
grid on;

También podría gustarte