Thanks, guys... but it seems these suggestions won't work.

Let me try to be more specific with a simple example:

Y<-list()
Y[[1]]<-function(u) sqrt(u)
Y[[2]]<-function(u) sin(u)
Y[[3]]<-function(u) 1/2*u

I wanted something equivalent to

Ybar<-function(u){
   1/3*(Y[[1]](u) + Y[[2]](u) + Y[[3]](u))
}

but with arbitrary length(Y) and without using any loops. Also, I can't
allow for discretization, since I must be able to evaluate Ybar at any u, as
I'm going to integrate it with the function "integrate".

Thanks again,

Eduardo Horta

	[[alternative HTML version deleted]]

