meufat <- function(K) {
  if (K==0) {
    fattmp = 1
  } else {
    fattmp <- K*meufat(K-1)
  }
  return(fattmp)
}