Неправильно выполняется формула в цикле

#include "stdafx.h"
long temp, timer;
double T1[20][20][20], T2[20][20][20];
double dT1[20][20][20], dT2[20][20][20];
float dtime=0.01, a1, a2;
int timer0;
 
 
int _tmain(int argc, _TCHAR* argv[])
{
temp=100;
timer=10;
 
int x=10;
int y=10;
int z=10;
int X=5;
int Y=5;
int Z=5;
int X1=2;
int Y1=2;
int Z1=2;
 
 for(X=1;X<=x;X++)
 for(Y=1;Y<=y;Y++)
 for(Z=1;Z<=z;Z++)
 {
    if (x=1)
      T1[X][Y][Z]=0; 
      dT1[X][Y][Z]=0;
       T2[X][Y][Z]=0;
       dT2[X][Y][Z]=0;
 }
 
 double PXYZ=dtime*dtime;
 
 for (timer0=1; timer0*dtime<=timer; timer0++)
{
 
    for(X=2;X<=(x-1);X++)
    for(Y=2;Y<=(y-1);Y++)
    for(Z=2;Z<=(z-1);Z++)
    {
      a1=0.005;
       a2=0.0009;
       dT1[X][Y][Z]=(T1[X-1][Y][Z]-2*T1[X][Y][Z]+T1[X+1][Y][Z])/PXYZ;
       dT1[X][Y][Z]=dT1[X][Y][Z]+(T1[X][Y-1][Z]-2*T1[X][Y][Z]+T1[X][Y+1][Z])/PXYZ;
       dT1[X][Y][Z]=dT1[X][Y][Z]+(T1[X][Y][Z-1]-2*T1[X][Y][Z]+T1[X][Y][Z+1])/PXYZ;
       dT1[X][Y][Z]=a1*dtime*dT1[X][Y][Z];
       dT2[X][Y][Z]=(T2[X-1][Y][Z]-2*T2[X][Y][Z]+T2[X+1][Y][Z])/PXYZ;
       dT2[X][Y][Z]=dT2[X][Y][Z]+(T2[X][Y-1][Z]-2*T2[X][Y][Z]+T2[X][Y+1][Z])/PXYZ;
       dT2[X][Y][Z]=dT2[X][Y][Z]+(T2[X][Y][Z-1]-2*T2[X][Y][Z]+T2[X][Y][Z+1])/PXYZ;
       dT2[X][Y][Z]=a2*dtime*dT2[X][Y][Z];
    
 
    }
        printf("  %d \n ", dT1);
}
 
 
 
 getchar();
    return 0;
}