Metodo Punto Fijo Matlab

El método de punto fijo consiste en una forma iterativa de resolver una ecuación de la forma f(x)=x . El método consiste en elegir una aproximación inicial X0 y realizar la iteración:

Xk+1 = f(Xk)

function puntofijo
global fun
fprintf('Método del punto fijo:\n');
fx = input('Ingrese la función original:\n','s');
fun=input('Ingrese la función despejada:\n','s');
x0=input('Ingrese el punto inicial:\n');
n=input('Ingrese el numero de iteraciones:\n');
it=0;
fprintf('  it       x0              x1          x0-x1');
while(it<n)
    it = it+1;
    x=x0;
    x1=eval(fun);
    fprintf('\n%3.0f%15.10f%15.10f%15.10f\n',it,x0,x1,abs(x1-x0));
    x0=x1;
end
fprintf('\n el punto fijo aproximado es=%10.6f\n',x1);
clf
hold on
fplot(fx,[-2-x1,x1+2]);
%ezplot('x');
%ezplot fun
%legend('y=x','fun')
plot(x1,0,'-*g')
grid on;
text(x1+0.2,0,'raiz')
hold off

Descargar archivo m