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

float x1,x2;
float d1,d2;
float h1,h2;
float a,e,score;
void JieFp();
void JieF();
void JieA();
void JieX2();


void main()
{
	printf("請輸入x1，x2，精度 e\n");
	scanf("%f %f %f",&x1,&x2,&e);

	h1=x1;
	h2=x2;
	do
	{
		x1=h1;
		x2=h2;
		JieFp();
		JieA();
		JieF();
		JieX2();
		if((fabs(h1-x1) < e) || (fabs(h2-x2)<e)) break;        //如果符合精度，則退出
		else 
		{x1=h1;
		 x2=h2;
		}
	;}while(1);

	printf("函數值為%e,x1為%e,x2為%e\n",score,h1,h2);  
}
void JieFp()
{
	d1=2*x1-x2-10;
	d2=2*x2-x1-4;
	d1=-d1;
	d2=-d2;
}
void JieA()
{
	float up,down;
	up=2*d1*x1+2*d2*x2-d2*x1-d1*x2-10*d1-4*d2;
	down=(2*d1*d1+2*d2*d2-2*d1*d2);
	a=up/down;
	a=-a;
}
void JieX2()
{
	h1=x1+a*d1;
	h2=x2+a*d2;
}
void JieF()
{
	score=h1*h1+h2*h2-h1*h2-10*h1-4*h2+60;
}