 Hi..
  
     i was trying to integrate the indicator funtion but had  problems when limits where negative or equal to the indicator condition
  
  my function is 
  ________________________
  fun1<-function(x){
      as.numeric(x>=2) 
  }
  _________________________
  
  which should be   Ind(x>=2)*x
  
  seems to work for the following two cases
  ------------------------------------------------------------
  > integrate(fun1,3,5)
  2 with absolute error < 2.2e-14
  
  > integrate(fun1,5,100)
  95 with absolute error < 1.1e-12
  ----------------------------------------------------------
   Does not work for the following
  
  > integrate(fun1,0,2)
  0 with absolute error < 0     ( i was expecting  = 2)
  
  > integrate(fun1,-1,5)
  3 with absolute error < 3.3e-14   (i was expecting =5)
  
  > integrate(fun1,-2,5)
  3 with absolute error < 5.3e-15    (i was expecting =5)
  
  Any suggestions?
  
  Thanks.
  
  Harsh,
  
  
  

		
---------------------------------

	[[alternative HTML version deleted]]


