Question

 constant <- function(x){1}

 integrate(constant,lower=1,upper=10)

this is probably more than ridiculously simple, but any reason as to why this isn’t working?

is there a different way I should be writing the function for constants?

Answer 1

You can use the Vectorize function to convert a non-vectorized function into the sort of function that integrate requires:

> constant <- function(x){1}
> Vconst <- Vectorize(constant)
>  integrate(Vconst,lower=1,upper=10)
9 with absolute error < 1e-13