// False Position 방법

#include <stdio.h>
#include <math.h>

#define true    1
#define false   0

double func(double x);

void main()
{
   int i, ok;
   double q0, q1, q, p;
   double p0 = 0.1;       // 왼쪽 경계치
   double p1 = 1.0;       // 오른쪽 경계치
   double tolerance = 0.0001;   // 오차의 한계
   double num = 40;             // 계산의 횟수
            
   printf("This is the Method of False Position.\n");

   q0 = func(p0); // f0 에 좌측 경계 함수값 입력
   q1 = func(p1); // f1 에 우측 경계 함수값 입력

   printf("  i            p              Func(p)\n");
   
   i = 1;

   while((i < num) && ok) {

      p = p1 - q1 * (p1 - p0) / (q1 - q0);
      q = func(p);

      printf("%3d   %15.9e   %15.9e \n", i+1, p, q);

      // 오차한계보다 작게 나오면 계산 끝냄. 
      if ((fabs(p - p1) < tolerance)) {
         printf("\nApproximate solution P = %12.9f \n", p);
         printf("with Func(P) = %12.9f\n", q);
         printf("Number of iterations = %3d\n", i+1);
         printf("Tolerance = %12.9e\n", tolerance);

	 ok = false;

      } else {

         i++;
         if((q*q1) < 0.) {
	    p0 = p1; q0 = q1;
         }
	 p1 = p;  q1 = q;
      }
      
   }

   if (ok) {
      printf("\nIteration number %3d", num);
      printf(" gave approximation %12.9f\n", p);
      printf("Func(P) = %12.9f not within tolerance : %15.9e\n", q, tolerance);
   }
}

// 해를 구하려는 함수
double func(double x)
{
   return (600*x*x*x*x-550*x*x*x+200*x*x-20*x-1);
}