[수치해석] Fixed Point Iteration Matlab Code

function [root, ea, iter] = FixedPoint(func, xr, es, maxiter)

if nargin < 3
    es = 1.0;
end

if nargin < 4
    maxiter = 50;
end

iter = 1;
fprintf('iter = %d: xr = %.4f\n', iter, xr);

while(1)
    xrold = xr;
    xr = func(xrold);
    ea = abs((xr - xrold) / xr) *100;
    
    iter = iter + 1;
    fprintf('iter = %d: xr = %.4f, es = %.2f\n', iter, xr, ea);

    if (ea <= es) || (iter > maxiter)
        break;
    end
end

root = xr;

FixedPoint.m

0.00MB