#include<stdio.h>
void main()
{
	
float a[11][12],t;
	int n,i,j,k,r;

    printf("LUx=b,Ux=y\n");
	printf("input unknown numbers<max=10>:");
    scanf("%d",&n);
	printf("input all quotieties:\n");

	for(i=1;i<=n;i++)
		for(j=1;j<=n+1;j++)
			scanf("%f",&a[i][j]);



for(r=1;r<=n;r++)
{
	for(i=r;i<=n;i++)//行
	{t=0;
	for(k=1;k<=r-1;k++)
		t=t+a[r][k]*a[k][i];
	a[r][i]=a[r][i]-t;
	}

	for(i=r+1;i<=n;i++)//列
	{t=0;
	for(k=1;k<=r-1;k++)
		t=t+a[i][k]*a[k][r];
	a[i][r]=(a[i][r]-t)/a[r][r];
	}


    t=0;                  //下三角
	for(k=1;k<r;k++)t=t+a[k][n+1]*a[r][k];
	a[r][n+1]=a[r][n+1]-t;
}

    

printf("y:\n");
for(i=1;i<=n;i++)
{
	for(j=1;j<=n+1;j++)
		printf("%f ",a[i][j]);
	printf("\n");
}






for(i=n;i>=1;i--)             //上三角
{t=0;
for(j=i+1;j<=n;j++)
t=t+a[i][j]*a[j][n+1];
a[i][n+1]=(a[i][n+1]-t)/a[i][i];
}

printf("x:\n");

for(i=1;i<=n;i++)
{
	for(j=1;j<=n+1;j++)
		printf("%f ",a[i][j]);
	printf("\n");
}

printf("The roots are:\n");  //輸出
for(i=1;i<=n;i++)
{printf("X%d=%f\n",i,a[i][n+1]);}


}