/*******************************************************/
/* Rozklad LU $VER: 1.0 (29.03.2002)                   */
/* Rozklad macierzy kwadaratowej A na iloczyn dwoch    */
/* macierzy trojkatnych L i U -- A=L*U                 */
/* (C) 2002 Paweł Marciniak <pmarciniak@lodz.home.pl> */
/*******************************************************/

#include <stdio.h>

#define MATRIXSIZE 3

double t[9]={3,2,3,6,8,10,9,22,30};

main()
{
	int i, j, k, suma=0;

	for( j=0; j<MATRIXSIZE; j++ )
	{
		for( i=0; i<MATRIXSIZE; i++ )
		{
			if(i>j) /* Macierz U */
			{
				suma=0;
				for( k=0; k<j; k++ )
				{
					if(i!=k)
						suma+=*(t+k*MATRIXSIZE+j) * *(t+i*MATRIXSIZE+k);
					else
						suma+=*(t+k*MATRIXSIZE+j);
				}
					*(t+i*MATRIXSIZE+j)=(*(t+i*MATRIXSIZE+j)-suma)/ *(t+j*MATRIXSIZE+j);
			}
			else /* Macierz L */
			{
				suma=0;
				for( k=0; k<i; k++ )
				{
					if(i!=k)
						suma+=*(t+k*MATRIXSIZE+j) * *(t+i*MATRIXSIZE+k);
					else
						suma+=*(t+k*MATRIXSIZE+j);
				}
				*(t+i*MATRIXSIZE+j)=(*(t+i*MATRIXSIZE+j)-suma);
			}
		}
	}

	/* wyswietlam wyniki */
	for(i=0; i<MATRIXSIZE; i++)
	{
		for(j=0; j<MATRIXSIZE; j++)
			printf("%f,",*(t+i*MATRIXSIZE+j));
		printf("\n");
	}

}