関係ないはずの変数をいじっても出力値が変わってしまう、、、

Question

至急質問があります。
以下のソースコードは偏微分方程式を解いているものなのですが、C言語について質問させていただきます。
下記の
//rho = 6782 - 0.27 * t[i][0];
//cp = 436 - 0.0813 * (t[i][0] - tm);
//lam = 94.61 + 4.41e-2 * (t[i][0] - tm);
//a = lam / rho / cp;
はコメントアウトにしていますが、このコメントアウトを外すと出力されるt[5][1]が変わってきます。
この式の意味や変数の意味は置いといて、このコメントアウト文以降でこれらの変数は使っていないのですが、このようなことは起こり得るのでしょうか。
また、このようなことが起こる理由について教えていただけますでしょうか。

環境

visual studio2019
windows 10

#include < stdio.h >
#define _USE_MATH_DEFINES
#include < math.h >

#define NX 10
#define NT 100
#define Tmax 1

void main(void)
{
	double t[NX + 2][2], dt, dx, r, a, time, pp, cp, lam, rho, tm;
	int i, imax;
	FILE* fp;
	fopen_s(&fp, "u0.xls", "w");

	dt = 1 / (double)NT;
	dx = 1 / (double)NX;
	imax = NX;
	a = 1.0;
	r = dt / dx / dx;
	time = 0;
	tm = 1850 + 273;

	fprintf(fp, "\t");
	for (i = 0; i <= imax; i++)
		fprintf(fp, "%lf\t", i * dx);
	fprintf(fp, "\n");

	//初期条件
	fprintf(fp, "0\t");
	for (i = 0; i <= imax; i++)
	{
		//三角形
		/*if (i <= imax / 2)
			t[i][0] = 2.0 / (dotble)imax * (double)i;
		else
			t[i][0] = -2.0 / (double)imax * (double)i + 2.0;*/

		//sin(x);
		t[i][0]=sin(M_PI*(double)i*dx);
		fprintf(fp, "%7.4lf\t", t[i][0]);

	}
	fprintf(fp, "\n");

	printf(" dx=%7.4lf dt=%7.4lf \n", dx, dt);
	printf(" r=%7.3lf \n", r);

	//時間積分
	do
	{
		//pp = 100 * exp(-pow((time - 0.5), 2) / 0.3);
		for (i = 1; i < imax; i++)
			//if(i==1)
			//rho = 6782 - 0.27 * t[i][0];
			//cp = 436 - 0.0813 * (t[i][0] - tm);
			//lam = 94.61 + 4.41e-2 * (t[i][0] - tm);
			//a = lam / rho / cp;
			t[i][1] = t[i][0] + (t[i + 1][0] - 2 * t[i][0] + t[i - 1][0]) * r;
			//else
				//t[i][1] = t[i][0] + (t[i + 1][0] - 2 * t[i][0] + t[i - 1][0]) * r;

		//境界条件
		//t[0][1] = t[1][1];
		//t[imax][1] =t[imax-1][1];
		t[0][1] = 0;
		t[imax][1] =0;

		time = time + dt;
		printf("t=%7.4lf\n", time);

		fprintf(fp, "%7.4lf\t", time);
		for (i = 0; i <= imax; i++)
			fprintf(fp, "%8.5lf\t", t[i][1]);
		fprintf(fp, "\n");

		//更新
		for (i = 0; i <= imax; i++)
			t[i][0] = t[i][1];

		printf("%8.3e\n",t[5][1]);
	} while (Tmax > time);
}

Answer 1

for (i = 1; i < imax; i++)

このforには{}がついていないため、その次の行がループの対象になります。
コメントアウトの有無でどこが対象になるか見比べればわかると思います。

Comments

[tomonari]

ありがとうございます。
非常に初歩的なミスで失礼しました。

Answer 2

forに{}がないからfor直下の式のみが繰り返されます

Comments

[tomonari]

ありがとうございます。
非常に初歩的なミスで失礼しました。