// 弦切法求根.cpp : 定义控制台应用程序的入口点。
//
#include "stdafx.h"
#include<iostream>
#include<cmath>
using namespace std;
double f(double x)
{
double y;
//y = x*x*x - 5*x*x + 16*x -80;
y = x*x*x + x*x - 4;
return y;
}
double xpoint(double x1 , double x2)
{
double y1 , y2,x;
y1 = f(x1);
y2 = f(x2);
x = (x1*f(x2) - x2*f(x1))/(f(x2)-f(x1));
return x;
}
double root(double x1 , double x2)
{
double y,x,y1;
x = xpoint(x1 , x2);
y = f(x);
y1 = f(x1);
if(fabs( y ) <= 0.0001)
return x;
if(y*y1 < 0)
x2 = x;
else
x1 = x;
return root(x1 , x2);
}
int _tmain(int argc, _TCHAR* argv[])
{
double x1 , x2 , x ,f1 , f2;
do{
cout<< "请输入x1,x2的值"<<endl;
cin>>x1>>x2;
f1 = f(x1);
f2 = f(x2);
}while(f1*f2 > 0);
x = root(x1 , x2);
cout<< x <<endl;
return 0;
}