f_rfc = np.zeros((N, N))