разложение матрицы на произведение двух треугольных в и н

[Monomah]

необходимо разложить матрицу на две треугольные(верхнюю и нижнюю(так чтобы их главная диагональ состояла из едениц) причем их произведение должно давать исходную матрицу.....
С верхней треугольной(b[][]) все нормально...а вот нижняя прблема(с[][])
при выводе какой-то бред выводит...порой её торкает и выводит матрицу b[][], а так обычно выводит первый элемент исходной потом какоето 100500значное число а потом одни нули...че делать ваще не понимаю
ЗЫ код на Си

#include <stdio.h>
#include <conio.h>
void main(void)
{
int i,j,k,n;
double z,a[25],b[5][5],c[5][5],a1[25];

scanf("%d",&n);

for (i=0;i<n;i++)
    scanf("%lf",&a[i]);	
for(i=1;i<n;i++)
	for(j=1;j<n;j++)
		{
		if (i<j) c[i][j]=0;
		else
				 if (i==j) {b[i][j]=1;c[i][j]=0;}
				 else b[i][j]=0;
		}
		
for (i=0;i<n;i++)
	for (j=0;j<n;j++)
		if (i<j) if (c[i][i]==0) return 0;
                 else 
                 {
                      b[i][j]=a[i*n+j];
                      for (k=0;k<(i-2);k++)
                          b[i][j]=a[i*n+j]-c[i][k]*b[k][j];
                      b[i][j]=b[i][j]/c[i][i];
                 }
        else 
        {
            c[i][j]=a[i*n+j];
            for (k=0;k<(j-2);k++)
                c[i][j]=c[i][j]-c[i][k]*b[k][j];
        }
b[0][0]=1;
for (i=0;i<n;i++)
    for (j=0;j<n;j++)    
        printf("%lf\n",c[i][j]);
for (i=0;i<n;i++)
    for (j=0;j<n;j++)    
        printf("%lf\n",b[i][j]);
getch();
}