/* lattice_init.c */

#define RANGEX 9
#define RANGEY 120

int main(void) {
  int r[RANGEX];
  int s[RANGEY];
  int t[RANGEY];

    int k, i;
    int d[RANGEX][RANGEY];
    int u[RANGEX][RANGEY];

    for (i = 0; i < RANGEX; i++) {
      for (k = 0; k < RANGEY; k++) {
 d[i][k] = 0;
 u[i][k] = 0;
      }
    }

    for (i = 0; i < RANGEX; i++) {
      r[i] = i;
    }

    for (i = 0; i < RANGEY; i++) {
      s[i] = RANGEY - i;
      t[i] = i;
    }

    for (k = 1; k < RANGEY; k++)
    {

 d[0][k] = s[k];
 u[0][k] = s[k];
 for (i = 1; i < RANGEX; i++)
 {
     d[i][k] = d[i - 1][k] + r[i] * u[i - 1][k - 1];
     u[i][k] = u[i - 1][k - 1] + r[i] * d[i - 1][k];
 }
 t[k] = d[RANGEX-1][k];
 printf("t[%d] = %d\n",k,t[k]);
    }
    return(0);
}
 
 

 back