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

#include "nr.h"

int main()
{
	int i,j;
	double d;
	double **input;
	int *index;

	/*   メモリ確保   */
	input = matrix(DIM,DIM);
	index = ivector(DIM);

	/*   入力行列の作成   */
	for (i=0;i<DIM;i++){
		for (j=0;j<DIM;j++){
			input[i][j] = 2;
		}
	}
	input[0][0] = 4;
	input[2][1] = 3;

	/*   入力行列の表示   */
	printf("元の行列表示\n");
	for (i=0;i<DIM;i++){
		for (j=0;j<DIM;j++){
			printf("%10lf",input[i][j]);
		}
		printf("\n");
	}
	printf("\n");

	/*   LU分解関数実行   */
	ludcmp(input, DIM, index, &d);
	
	/*   行列式計算と表示   */
	for (i=0;i<DIM;i++){
		d *= input[i][i];
	}
	printf("行列式：%lf\n",d);

	/*   結果行列の表示   */
	printf("変換後行列の表示\n");
	for (i=0;i<DIM;i++){
		for (j=0;j<DIM;j++){
			printf("%10lf",input[i][j]);
		}
		printf("\n");
	}

	/*   メモリ開放   */
	free_matrix(input);
	free_ivector(index);

	return 0;
}