разложение матрицы

[Monomah]

необходимо разложить матрицу на две треугольные(верхнюю и нижнюю(и их главная диагональ состояла из едениц) причем так чтобы их произведение довало исходную матрицу.....
С верхней треугольной(b[][]) все нормально...а вот нижняя прблема(с[][])
при выводе какой-то бред выводит...на главной диагонали откуда-то нули берутся либо дикие 9-27значные числа...хотя по формуле все делаю верно( помогите найти ошибку

#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();
}

[Д_М]

n - размер матрицы?
a - исходная матрица?

for (i=0;i<n;i++)
    scanf("%lf",&a[i]);

Здесь вводится только первая строка матрицы, остальные строки не инициализированы.

И почему a[25], а не a[5][5]? Неудобно же так.

[Monomah]

Сообщение от Д_М

n - размер матрицы?
a - исходная матрица?

for (i=0;i<n;i++)
    scanf("%lf",&a[i]);

Здесь вводится только первая строка матрицы, остальные строки не инициализированы.

И почему a[25], а не a[5][5]? Неудобно же так.

Потому что в дальнейшем в задаче необходимо будет динамическое
выделение памяти что в двумерном массиве не очень удобно.... поэтому удобнее использовать одномерный массив