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

/*************** 表5.2の数値解を求めるプログラム ****************

	・式(5.35)を用いて解を求めます．

	・ただし，i=0,1,2,3,4 とします．
		  l=0,1,2,3,4

*****************************************************************/

#define eps pow(10,-8)		//収束判定条件
#define kmax 1000		//収束判定条件との比較回数の上限

int main(void){
	int i, l, k;
	double L2n, sum1, sum2;
	double a[5][5] = {{10,3,1,2,1},{1,19,2,-1,5},{-1,1,30,1,10},						{-2,0,1,20,5},{-3,5,1,-2,25}};
	double b[5] = {-22,27,89,-73,22};
	double x1[5] = {0,0,0,0,0}, x2[5];

	for(k=1; k<kmax+1; k++){
		for(i=0; i<5; i++){
			sum1 = 0;
			sum2 = 0;
			for(l=0; l<i; l++)
				sum1 += -a[i][l]*x2[l];
			for(l=i+1; l<5; l++)
				sum2 += -a[i][l]*x1[l];
			x2[i] = (sum1+sum2+b[i])/a[i][i];
		}

		for(i=0; i<5; i++){
			printf(" k=%2d, x2[%2d]=%13.6e\n"							, k, i+1, x2[i]);
		}

		L2n = 0;
		for(i=0; i<5; i++)
			L2n += pow((x2[i]-x1[i]),2);

		if(fabs(L2n) < eps)		return 0;

		for(i=0; i<5; i++)
			x1[i] = x2[i];
	}
	return 1;
}