Вложенные циклы

T_complex ** calc_reflected_beam ( T_complex *fr, T_complex **R,                        
                                 double fr_begin, double fr_end, double fr_step,    
                                 double f_begin, double f_end, double f_step,  
                                 double l_begin, double l_end, double l_step)
{
    double x, p, nkappa;
    int  i, k, j, n, l;
    int max_NNN_R_BEAM;
 
    max_NNN_R_BEAM = (int) floor ( ( l_end - ( - l_end) ) / l_step)+2;
    T_complex ** R_BEAM;
    R_BEAM = (T_complex **)malloc(max_NNN_R_BEAM * sizeof(T_complex*));
    for(int i=0; i < max_NNN_R_BEAM; i=i++){
    R_BEAM[i] = (T_complex*)malloc(max_NNN_R_BEAM*sizeof(T_complex));
    }
    i = 0;
    n = 0;
    for (nkappa = l_begin; nkappa < l_end; nkappa = nkappa + l_step )
    {
        k=0;
        for (x = f_begin; x < f_end; x = x + f_step)
        {
            R_BEAM[i][k].re = 0;
            R_BEAM[i][k].im = 0;
            j=0;
            l=0;
            for ( p = fr_begin; p < fr_end; p = p + fr_step)
            {
                R_BEAM[i][k].im += (fr[l].re*R[n][j].im * cos ( p * x ) + fr[l].im*R[n][j].re * cos ( p * x ) + fr[l].re*R[n][j].re * sin ( p * x ) - fr[l].im*R[n][j].im * sin ( p * x ))* fr_step/(dsqrt_2_pi);
                R_BEAM[i][k].re += (fr[l].re*R[n][j].re * cos ( p * x ) - fr[l].im*R[n][j].im * cos ( p * x ) - fr[l].re*R[n][j].im * sin ( p * x ) - fr[l].im*R[n][j].re * sin ( p * x ))* fr_step/(dsqrt_2_pi);
                j += 1; 
                l += 1;
            }
            k += 1;
        }
        n += 1;
        i += 1;
    }
    return ( R_BEAM );
}