Thankyou Brian > integrandtotest <- function(x) colSums(outer(t, x, "^")) > This worked fine >in fact this integration can be done analytically (it is a > sum of exponentials). I know but i have another five terms in the integrand as well as the one that was causing the trouble Many, many thanks Alan