R help - Jul 2008 - Using integrate

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]]

On 21/07/2008 4:56 PM, Ayman Oweida wrote:
> 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.
But it is very simple to define new functions.  So you could do

z[i] <- integrate( function(x) f(x) * (t [i+1] - t[i]), i-1, i )

Duncan Murdoch