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
______________________________________________
R-help@r-project.org 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.
