The following code works fine: # g is the function that returns the square of a number g <- function(y) y^2
# f1 is a function that takes one function # as argument, and returns another function f1 <- function(f) function(x) f(x+1) - f(x) # h(x) is g(x+1) - g(x) or 2x + 1 h <- f1(g) # h(1) = 3 # h(2) = 5 # h(3) = 7 So far, so good. But why: h shows: function(x) f(x+1)-f(x) <environment: 0264BE84> I don't get it. h should show function(x) g(x+1)-g(x) or something like that. Alberto Monteiro ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.