comprend rien :-(

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

#define N 3
#define M 3

float C[M][N];
int i,j,k,l;

main(){
  i=1;
  j=1;
  while (i<N+1){
    while (j<M+1){
      C[i][j] = (i*j*(j+1)) / ((i+j-1)*(1- pow((2.0/3.0),(j+1))));
      j++;
    }
    i++;
    j=1;
  }
  for(i=1;i<N+1;i++){
    for(j=1;j<M+1;j++){
      printf("%f  ", C[i][j]);
    }
    printf("\n");
    }
}

Je compile, tout se passe bien...
Mais lorque je lance le "programme", j'ai pour le tableau C:
3.600000  4.263158  4.984615
3.600000  5.684210  7.476923
3.600000  6.394737  0.000000

Pourquoi est ce que la valeur C[3][3] est nulle??? alors que lorsque j'écris 
le programme en pascal j'ai:
8.972307692307693?!

Merci d'avance pour toute piste...