#include "stdio.h"
void main()
{
    float A[15][15],L[15][15],U[15][15];
    int i,j,k;
    float s;
    A[1][1]=2;
    A[1][2]=1;
    A[1][3]=1;
    A[2][1]=1;
    A[2][2]=3;
    A[2][3]=2;
    A[3][1]=1;
    A[3][2]=2;
    A[3][3]=2;
   for(j=1;j<=3;j++)
       U[1][j]=A[1][j];
    for(i=2;i<=3;i++)
       L[i][1]=A[i][1]/U[1][1];
    for(i=2;i<=3;i++)
    {
        for(j=i;j<=3;j++)
        {
            s=0;
            for(k=1;k<=i-1;k++)
               s=s+L[i][k]*U[k][j];
            U[i][j]=A[i][j]-s;
        }
        for(j=i+1;j<=3;j++)
        {
            s=0;
            for(k=1;k<=i-1;k++)
               s=s+L[j][k]*U[k][i];
            L[j][i]=(A[j][i]-s)/U[i][i];
        }
    }
    for(i=1;i<=3;i++)
    {
       for(j=1;j<i;j++)
          printf("%4.1f",L[i][j]);
       printf("  1\n");
    }
    for(i=1;i<=3;i++)
    {
           for(j=1;j<=3;j++)
           printf("%4.1f",U[i][j]);
        printf("1\n");
    }
   

}