Ошибка: Unhandled exception at 0x772115de in TimeShift.exe: 0xC0000005: Access violation reading location 0x00041000.

T_complex ** calc_T ( T_complex *fr, // используемые функиции
                                 double fr_begin, double fr_end, double fr_step,    // границы и шаг по p
                                 double l_begin, double l_end, double l_step )      // границы и шаг по x
 
{
    double x, p, ndeltatetta, nq, l, t, nkappa;
    double j1=1.0, j2=1.0, j3=1.0, j4=1.0;
    int i, j;
    T_complex inter,  e_ipx;
    //int NNN_fr, NNN_R_BEAM;
    int max_NNN;
    max_NNN = (int) floor ( ( fr_end - fr_begin ) / fr_step) + 2;
    T_complex ** R;
    R = (T_complex**)malloc(max_NNN * sizeof(T_complex*));
    for(int i=0; i < max_NNN; i=i++){
        R[i] = (T_complex*)malloc(max_NNN*sizeof(T_complex));
    }
    l=1.0;
    i=0;
    for (nkappa = 0; nkappa < 0.5; nkappa = nkappa + fr_step )
    {
        j=0;
        R[i][j].re=0;
        R[i][j].im=0;
            //for (  x = 0; x < 10; x = x + x_step)
            //{
                inter.im = 0;
                inter.re = 0;
            //  NNN_fr = 0;
            //  e_ipx.im = 0;
            //  e_ipx.re = 0;
                for ( p = fr_begin; p < fr_end; p = p + fr_step )
                {
                    ndeltatetta = p - 1;
                
                    if ( ( ndeltatetta * ndeltatetta ) < ( nkappa * nkappa ) )
                    {
                        nq = sqrt( nkappa * nkappa - ndeltatetta * ndeltatetta );
                        if ( ndeltatetta * sinh ( nq * l ) < 0 )
                    { 
                        j1 = -1.0;
                    }
                    if ( ndeltatetta * sinh ( nq * l ) > 0 )
                    { 
                        j1 = 1.0;
                    }
                    R[i][j].re = ( nq * nq * cosh ( nq * l ) *  cos ( pi*l ) + nq*ndeltatetta*sin(pi*l)*sinh(nq*l))/( nq * nq * cosh ( nq * l) * cosh ( nq * l) + ndeltatetta * ndeltatetta * sinh( nq*l)* sinh(nq*l) );
                    R[i][j].im = ( nq * ndeltatetta * sinh ( nq * l ) * cos ( pi*l ) - nq*nq*sin(pi*l)*cosh(nq*l))/( nq * nq * cosh ( nq * l) * cosh ( nq * l) + ndeltatetta * ndeltatetta * sinh( nq*l) * sinh(nq*l) );
                    }
                    j=j+1;
                }
                i=i+1;
            }
        return (R);
}

соut<<"i= "<<i<<"  "<<... и тд. 
system("pause");