Bisection.cpp

Go to the documentation of this file.

// to compile:  g++ -I. -lm Bisection.cpp<br />
// to compile with debugging: g++ -g -I. -lm Bisection.cpp<br />
#include "Bisection.h"
#include <math.h>

class secant_functor {
public:
  double _offset;
  
  secant_functor( double offset) : _offset(offset) { }

  double f( double x) {
    return - (exp( x) - _offset);
  }
  
};

int main() {
  double x0=1.0, x1=10.0, fx0, fx1;
  secant_functor sf( 5);

  stradle_value( x0, x1, fx0, fx1, 0.0, sf, &secant_functor::f);

  Bisection< secant_functor, double > 
    foobar( 1.0, 10.0, fx0, fx1, 0.0, 1e-8, 1e-8, 100, sf, &secant_functor::f);

  foobar.do_iteration( true);

  cout << "***************************" << endl;

  cout << foobar.get_x_mid() << endl;

  return 0;
}

---