Thread: Mutidimensional Array problem

Code:

  for(trial = 0; trial < N; trial++)
    {
      b[trial] = 0.0;
      nstep = 0;
      t = 0.0;
      dis = 0;
      counter = 0;
      while(t < to)
	{
	  tau = tao * (1 + g * t); //tao, g, t are all entered at the command prompt
	  dt = -log(a) * tau; //a is entered at the command prompt
	  u = gsl_ran_flat (r, 0.0, 1.0);
	  b[nstep * maxstep + dis]++; //maxstep is pre-determined, as is nstep
	  if(u > a)
	    {
	      dis++;
	    }
	  t = t + dt;
	  nstep++;
	}
    }
  
  for(i = 0; i < maxstep; i++)
    {
      for(j = 0; j < maxstep; j++)
	{
	  fprintf(stderr, "%lf ", b[j * maxstep + i]);
	}
	  fprintf(stderr, "\n");
    }

Code:

20.000000 4.000000 0.000000 2.000000 1.000000 0.000000 0.000000 0.000000 
0.000000 7.000000 2.000000 6.000000 5.000000 4.000000 2.000000 0.000000 
0.000000 0.000000 0.000000 10.000000 4.000000 3.000000 2.000000 2.000000 
0.000000 0.000000 0.000000 2.000000 10.000000 9.000000 8.000000 6.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 4.000000 5.000000 7.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 3.000000 4.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 1.000000 
0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000

Supposedly, the sum of each column should be 20 (N = 20), however, column 2 and 3 do not add up to 20...why not?