#include "f.hh";
#include <math.h>
#include <iostream>
#include <iomanip>
using namespace std;



void secante(double x0, double x1, double e, int MaxSteps,  double zero[])  {
double x, y0,y1;
double  n=0;

 while (fabs(x1-x0)>e){
 y0=f(x0);
 y1=f(x1);

 // Controllo per y1 -y0 troppo piccolo!! 
 // if(fabs(y1-y0)<0.001) y1 += 0.001;
 x = (x0*y1 - x1*y0)/(y1-y0);
  n++;

  //Test di convergenza

    x0=x1;
    x1=x;

    //  cout<< "x0= "<< setprecision(6)<<x0 << setw(15) << "  x1= " <<setprecision(6)<<  x1<< setw(15) << "fabs = "<< fixed <<setprecision(6)<< fabs(x1-x0) << " e = " << setprecision(6) <<  e << endl;

    //  cout<< setprecision (6) << "x0= "<< setprecision(6)<< x0<< setw(15) << setprecision(6)<< " x1= " << x1 << " Npass = "<< n << endl;
  //  if(n>=MaxSteps) break; 
 } 
  zero[0]=x1;
  zero[1]= y1;
  zero[2]= n;


}