//*************************************
//程序功能：Steffenson加速收斂法  求線性方程的解
//
//作    者：杭永東
//
//日    期: 2002.03.26
//*************************************

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

//*************************************
//
//注:更改此宏定義的表達式可以解不同的線性方程
//   根據宏定義要求,x必須以(x)形式出現式中
//
//*************************************

#define fi(x) atan((x))

long double x0,e;

long double Dd(long double x0)
{
	printf("\n%.9f",x0);
	long double y=fi(x0);
	long double z=fi(y);
	long double x1=x0-pow((y-x0),2.0L)/(z-2.0L*y+x0);

	if((((x1-x0)>0)? (x1-x0): (x0-x1))<e)
		return x1;
	return Dd(x1);
}

main()
{
	cout<<"Please input x0,e:";
	cin>>x0>>e;
	printf("\n%.9f",Dd(x0));
}