set.seed(1) x0 = sort(rnorm(10)) f1 = function(x) x^2 f2 = dnorm f3 = function(x) x^3 listaFs = list(f1, f2, f3) results = do.call(cbind, lapply(listaFs, function(myf, myarg) myf(myarg), x0)) matplot(x0, results, type='l')