# Parámetros de entrada
N <- as.numeric(readline(prompt = "Introduce el numero total de nutrientes "))
M <- as.numeric(readline(prompt = "Introduce el numero de condiciones experimentales "))
K <- as.numeric(readline(prompt = "Introduce el numero de factores de ajuste "))  
x <- matrix(seq(1,N * M), nrow = N, ncol = M)  # Matriz N×M con valores específicos 
w <- seq(1,K)                                    #de la concentración de cada variable x[i,j]
f <- seq(1,N)
c <- seq(1,M)
d <- seq(1,K)
#No tenemos los datos exactos de cómo son estos vectores, así que simplemente hemos puesto 
#una secuencia. No es lo importante.


P<-0

for (i in 1:N) {
  Prod<-1
  for (j in 1:M) {
    Sum1<-0
    for (k in 1:K) {
      # Cálculo de v[i,j,k] (RESOLUCIÓN 2)
      v<-f[i]*c[j]*d[k]
      Sum1<-Sum1+(w[k]*v)
    }
    Prod<-Prod*(x[i, j] + Sum1)
  }
  P<-P+Prod
}

print(P)
print(v)