meufatI <- function(K) {
    fattmp = 1
    i = K
    while (i > 0) {
      fattmp <- i*fattmp
      i = i - 1
  }
  return(fattmp)
}