Xk+1 = f(Xk)
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
Este codigo si corre
ResponderBorraraea
ResponderBorrarmanito piola
ResponderBorrarNo puedo descargar el archivo
ResponderBorrar