Ошибка памяти

#pragma hdrstop
#include<stdio.h>
#include<conio.h>
 
//---------------------------------------------------------------------------
 
#pragma argsused
int main(int argc, char* argv[])
{
  FILE *tr;
  FILE *tr_an;
  char str[256],str_a[256];
  int N=100,i,j,T=100;                              //N-number of spatial cells, T - number of time cells  
  float h,tau,kappa=0.25,a=1,L=1;             //a-speed,L-length of colclation area, h - spatial step 
  float f_num[100]={0},f_num1[100]={0};  //numeral solution 
  float f_an[100]={0},f_an1[100]={0};      //analitic solution
  float x=0,t=0,ta=0.375;
  h=L/N;
  tau=kappa*h/a;
 
 
  //tr=fopen("transfer_eq.txt","w");
  //tr_an=fopen("transfer_eq_a.txt","w");
  for(i=0;i<N;i++) {                  //initial conditions
    if ((x>=0)&&(x<=L/2)) {
      f_num[i]=2;
      f_an[i]=2;
    }
    else if ((L/2<x)&&(x<=1)) {
      f_num[i]=1;
      f_an[i]=1;
    }
    x+=h;
  }
  for(j=0;j<1;j++) {                  //boundery conditions
    f_num[j*N]=2;
    f_an[j*N]=2;
    t+=tau;
    if ((0<=t)&&(t<L/2/a)) {
      f_num[j*N+(N-1)]=1;
      f_an[j*N+(N-1)]=1;
    }
    else if (L/2/a<=t) {
      f_num[j*N+(N-1)]=2;
      f_an[j*N+(N-1)]=2;
    }
  }
 
  for(j=0;j<T;j++) {                       //loop of time coordinates
    sprintf(str,"data_num/a_num%d.txt",j); //recording to files that
    sprintf(str_a,"data_an/a_an%d.txt",j);
    tr=fopen(str,"w");                     //are in the folder 'data'
    tr_an=fopen(str_a,"w");
    for(i=1;i<N;i++) {                     //loop of spatial coordinates
      f_num1[i]=-kappa*(f_num[i]-f_num[i-1])+f_num[i];
      fprintf(tr,"%.2f  %.2f\n",(float)i*h,f_num1[i]);
      f_an1[i]=f_an[i-a*ta];
      fprintf(tr_an,"%.2f  %.2f\n",(float)i*h,f_an1[i]);
    }
    for(i=1;i<N;i++) {
      f_num[i]=f_num1[i];
      f_an[i]=f_an1[i];
    }
  }
 fclose(tr);
 fclose(tr_an);
  getch();
  return 0;
}

    tr=fopen(str,"w");                     //are in the folder 'data'
    tr_an=fopen(str_a,"w");