Segmentation fault en c

slt a tous j'ai un soucis avec mon programme,ça compile bien mais le problème est au niveau de l'exécution,il m'affiche segmentation fault.Je vous laisse le code source
#include <math.h>
#include <stdio.h>
int main()
{
int i,j,Nz,Nt;
long double Na[1000][1000],Ne[1000][1000];
double N,deltaz,deltat,E,mua,mue,d,T;
FILE* ecrit1;
ecrit1= fopen("sortie.txt","w");

/*Données physiques*/

/*mua=3.0e-3;*/
/*mue=3.0e-2;*/

/*paramètres du calcul*/
E=1.0e3;
N=50.0;
d=1.0;
Nt=100;
Nz=100;
T=10.0*d/(mue*E);
deltaz=d/Nz*1.0;
deltat=T/Nt*1.0;

for(i=0;i<=Nz;i++)
{

Na[i][0]=N;
Ne[i][0]=N;

}


for(j=0;j<=Nt;j++)
{

Na[0][j]=0;
Ne[1][j]=0;

}

for(j=0;j<=Nt;j++)
{
for(i=0;i<=Nz;i++)
{
Ne[i][j+1]=0;
Na[i][j+1]=0;
Ne[i][j+1]=(1-(E*deltat)/deltaz)*Ne[i][j]+ (E*deltat/deltaz)*Ne[i+1][j];
Na[i][j+1]=(1-(E*deltat)/deltaz)*Ne[i][j]+ (E*deltat/deltaz)*Ne[i-1][j];
}
}

fprintf( ecrit1,"%i %i %e %e\n",i,j,Na[i][j],Ne[i][j]);
fclose(ecrit1);
return 0;