Dear R community,

I am trying to code a fairly complex equation for optim().  My current
approach is too slow for optim().

I have a function that takes a double integral (hopefully correctly) across
two terms, e.g.,

 doubleint <- function(c1,c2) {integrate(function(y) {
   sapply(y, function(y) {
   integrate(function(x) {a*(c1+x)+b*(c2+x)}, boundsx[1], boundsx[2])$value
   })
   }, boundsy[1],boundsy[2])
   }

I would like to rapidly evaluate this function over vectors c1 and c2 of
equal length where the double integral is calculated for each (matching)
element c1_i & c2_i.  At present I get length mismatch errors.  Furthermore,
mapply() takes too long.  Can I expedite this evaluation or do I need to
reformulate my approach?

Thanks in advance for your help,
Scott C. Stark
University of Arizona

	[[alternative HTML version deleted]]

