/* fixedpoint.c  -- program to perform fixed point iteration */
#include <stdio.h>
#include <math.h>
#include "myfunc.h"

int main(void)
{
	double p0,pnew;
	int i=1;
	int status=ITERATING;

	printf("Fixed Point Iteration\n");
	printf("Please enter an initial guess:\n");
	scanf("%lf",&p0);
	printf("\nIteration\tx\t\tg(x)\n");
	while (i <= ITMAX && status != CONVERGED)
	{
		pnew = g(p0);
		printf("%d\t\t%6.8lf\t\t %6.8lf\n",i,p0,pnew);

		if((fabs(pnew-p0))<=TOL)
		   status = CONVERGED;
		p0 = pnew;
		i++;

	}  //end of while block

	if(status == CONVERGED)
	   {
		   printf("\nMethod converged in %d iterations.",i);
		   printf("\nRoot is %6.8lf",pnew);
	   }
	else
	   printf("\nMethod did not converge after %d iterations.\n",ITMAX);
	return 0;

}