Fixed-point iteration Method for Solving non-linear equations in MATLAB(mfile)

% Fixed-point Algorithm 

 % Find the fixed point of y = cos(x).

 g = @(x) cos(x);

 p0 = input('Please enter initial approximation, p0:  ');

 n = input('Please enter no. of ierations, n: ');

 tol = input('Please enter tolerance, tol: ');

 i = 1;

 while i <= n

     p = g(p0);

     if abs(p-p0) < tol

        fprintf('\nApproximate solution p = %11.8f \n\n', p)

        break;

     else

        i = i+1;

        p0 = p;

     end

 end

Output:

>> Fixed_point_method_m

Please enter initial approximation, p0:  1

Please enter no. of ierations, n: 30

Please enter tolerance, tol: 0.001

Approximate solution p =  0.73876032 

>>