用弦截法求方程的解

#include <iostream>
#include <iomanip>
#include <cmath>
double f(double x);
using namespace std;
int main()
{
 double x0=0.5,x1=1,x2;
 cout<<x0<<endl;
 cout<<x1<<endl;
 for(;fabs(x1-x0)>=0.000001;)
 {
  x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0));
     x0=x1;
  x1=x2;
        cout<<x2<<endl;
 }
 return 0;
}

double f(double x)
{
 double y;
 y=x*x*x*exp(x)-2;
 return(y);
}