上の続き（最後）

for(t=0;t<tmax;t++) {

/*境界条件として回りの温度0*/
for(j=0;j<N;j++) {
T0[j][j][0]=0.0;
T0[j][j][N-1]=0.0;
T0[0][j][j]=0.0;
T0[N-1][j][j]=0.0;
T0[j][0][j]=0.0;
T0[j][N-1][j]=0.0;


}

for(j=1;j<N-1;j++) {
for(k=1;k<N-1;k++){
for(l=1;l<N-1;l++) {
T1[j][k][l]=T0[j][k][l]+0.25*(T0[j+1][k][l]+T0[j-1][k][l]+T0[j][k+1][l]+T0[j][k-1][l]
+T0[j][k][l+1]+T0[j][k][l-1]-6.0*T0[j][k][l]);
}
}
}