どうしても分からなくて悩んでいます。
以下に示すプログラムをビルドしますと　debug assertion failed!という警告が出て実
行することができません。
どうしてでしょうか？どなたかご教授願えませんか。
プログラミング環境はwindows vc++ 6.0です。
#include stdafx.h
#include <stdio.h>
#define M 4
main(int argc, char *argv[])
{
double a[M][M];
double b[M];
double c[M];
double l[M][M];
double u[M][M];
double x[M];
int i,j,k;
FILE *fp;
fp = fopen(argv[1], r); /* ファイルより入力行列を読み込む*/
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
fscanf(fp,%lf,&a[i][j]);
}
fscanf(fp,%lf,&b[i]);
}
fclose(fp);

for(i = 0; i < M; i++){ /* L行列、U行列を1と0で初期化*/
for(j = 0; j < M; j++){
u[i][j] = 0;
if(i == j)
l[i][j] = 1;
else
l[i][j] = 0;
}
}
for(i = 0; i < M; i++){
for(j = i; j < M; j++){ /* U行列の生成*/
u[i][j] = a[i][j];
for(k = 0; k < i; k++){
u[i][j] -= u[k][j] * l[i][k];
}
}
for(j = i+1; j < M; j++){ /* L行列の生成*/
l[j][i] = a[j][i];
for(k = 0; k < i; k++){
l[j][i] -= u[k][i] * l[j][k];
}
l[j][i] /= u[i][i];
}
}

for(i = 0; i < M; i++){ /* c行列の生成*/
c[i] = b[i];
for(j = 0; j < i; j++){
c[i] -= l[i][j] * c[j];
}
}
for(i = M-1; i >= 0; i--){ /* x行列の生成*/
x[i] = c[i];
for(j = M-1; j > i; j--){
x[i] -= u[i][j] * x[j];
}
x[i] /= u[i][i];
}
printf(入力行列\n); /* 入力行列の出力*/
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
printf(%10.5lf ,a[i][j]);
}
printf(%10.5lf\n,b[i]);
}
printf(\nL行列\n); /* L行列の出力*/
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
printf(%10.5lf ,l[i][j]);
}
printf(\n);
}
printf(\nU行列\n); /* U行列の出力*/
for(i = 0; i < M; i++){
for(j = 0; j < M; j++){
printf(%10.5lf ,u[i][j]);
}
printf(\n);
}
printf(\n);
for(i = 0; i < M; i++){ /* 解の出力*/
printf(x%d = %10.5lf\n,i,x[i]);
}
}