快速弦截法解方程

/**
***快速弦截法***
 
  公式：x(k+1) = x(k) - ( x(k) - x(k-1) ) * f( x(k) ) / f( x(k) - f(x(k-1) )

  属性：两步迭代法

  描述：计算x(k+1)时需要利用前两步信息x(k),x(k-1).免去了Newton法中需要求解一阶导函数的繁琐

  《数值分析简明教程》-2 Editon -高等教育出版社 -page 140 -无算法流程图

  代码维护：2005.6.14  DragonLord
**/

#include<iostream.h>
#include<math.h>
#include<stdio.h>

//范例程序中方程为：x*exp(x)-1=0

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

int main()
{
 double x0,x1,x2,e,temp;
 int N,k;
 while(cin>>x0>>x1>>e>>N)
 {
  k=1;
   
  loop: if(f(x1)-f(x0)==0){cout<<"函数异常！"<<endl;break;}