for (i in X) {
  for (j in coeffic) {
      Wy_xR1[i][j] = 1/(SigmY * sqrt(1 - R[1]^2)* sqrt(2 * pi)) *
        exp(1)^(-((((My_xR1 - coeffic[j] * SigmY) - My_xR1)^2)/(2 * Dy_xR1)))
      Wy_xR2[i][j] = 1/(SigmY * sqrt(1 - R[2]^2)* sqrt(2 * pi)) *
        exp(1)^(-((((My_xR2 - coeffic[j] * SigmY) - My_xR2)^2)/(2 * Dy_xR2)))
  }
}