найти ошибку в коде С задача на обратную матрицу

[Monomah]

Задача найти обратную матрицу методом разложения её на произведение двух треугольных

void main (void)
{
double z,a[1024],b[32][32],c[32][32],a1[1024];
int i,j,k,p,m,s,l;
scanf("%d",n);
for (i=1;i<n;i++)
{scanf("%lf",z);
a[i]=z;}
p=0;m=0;
for (i=1;i<n;i++)
for (j=1;j<n;j++)
for (i=k;k<(j-1);k++)
c[i][j]=a[i+i*n+j]-(p+c[i][k]*b[k][j]);
for (i=1;i<n;i++)
for (j=1;j<n;j++)
for (i=k;k<(i-1);k++)
b[i][j]=(a[i+i*n+j]-(m+c[i][k]*b[k][j]))/c[i][i];
p=0;m=0;s=0;
for (i=1;i<n;i++)
for (j=1;j<n;j++)
for (i=k;k<(j-1);k++)
if (i>j) a1[i]=(-1)*(p+a1[i+i*n+k]*c[k][j])/c[j][j];
for (i=1;i<n;i++)
for (j=1;j<n;j++)
for (i=k;k<(i+1);k++)
{if (i==j) a1[i+i*n+i]=(1-(m+a1[i+i*n+k]*c[k][i]))/c[i][i]; 
if (i<j) a1[i+i*n+j]=(-1)*(b[i][k]*a1[k+k*n]+m);}
for(i=1;i<n;i++)
{a[i]=z;
printf("%lf",z);}
}