def METHOD_NAME(nIter, size_R, cols, rows, J, dN, dS, dE, dW, c, r1, r2, c1, c2, Lambda):
    for ii in range(nIter):
        Sum = 0.
        Sum2 = 0.
        for i in range(r1, r2+1):
            for j in range(c1, c2+1):
                tmp = J[i * cols + j]
                Sum += tmp
                Sum2 += tmp*tmp
        meanROI = Sum / size_R
        varROI = (Sum2 / size_R) - meanROI*meanROI
        q0sqr = varROI / (meanROI*meanROI)
        srad_1(q0sqr, cols, rows, J, dN, dS, dE, dW, c)
        srad_2(cols, rows, J, dN, dS, dE, dW, c, Lambda)