BOUNDING_PHASE_BOOM.cpp

/*


BOUNDING PHASE METHOD
Sudharsan Neelamegam


**** NOTE :Below float fun(float x), put the funation as given in the question in terms of x

*/
#include <bits/stdc++.h>
#include <math.h>
#include <iomanip>

using namespace std;


float fun(float x)
{ float result;
  result=pow(100-x,2);//ENTER FUNCTION IN TERMS OF X AS GIVEN IN QUESTION PROPERLY HERE
  return result;
    
}
float fun();

int main()
{
  float x0,k=0,xk,xkp,n,fxk,fxkp,fmm,fpp,d,f0,v=0;
  float c,con=0;
  cout<<"Enter the initial guess value x0 and increment value :";
  cin>>x0>>d;
  
  
  
  cout<<"k   xk      x(k+1)      fxk     fx(k+1)     cond"<<endl;
    fmm=x0-d;
    f0=x0;
    xk=f0;
    fpp=x0+d;
    xkp=fpp;
  while(v==0)
  {
      
      float a=fun(fmm),b=fun(xk),c=fun(xkp); 
      if((a>=b) && (b>=c))
      {
          d=d*1;
      }
      else
      {
          d=d*(-1);
      }
      if(c>b)
      { con++;
        cout<<k<<"  "<<xk<<"  "<<xkp<<"  "<<b<<"  "<<c<<"  "<<con<<endl;
        break;
      }
     cout<<k<<"  "<<xk<<"  "<<xkp<<"  "<<b<<"  "<<c<<"  "<<con<<endl;
     
     xk=xkp;
     xkp=xk+((pow(2,k+1))*d);
     k++;
     
  }
  cout<<"Req range is: ("<<(xk+(pow(2,k-1)*d))<<","<<xkp<<")";

return 0; 
}