
I have a function, say:
f<-function(x)  exp(x)

and I would like to obtain the integration of the function while adding a few operations within the integration (I need point-by-point integration), say:

t<-seq(0,40, by=1)
z<-array(0,length(t))
for (i in 1:40){ 
z[i] <- integrate ( f *  (t [i+1] - t)   ,  (i-1)  ,  i   ) 

In R, I can only have the function in 'integrate'
i.e.  
z[i] <- integrate ( f  ,  (i-1)  ,  i   ) 

but I want to add additional operations within the 'integrate' command, which I cannot just add to the function f.


any input is appreciated.

Ayman






      
	[[alternative HTML version deleted]]

