Cubature.jl icon indicating copy to clipboard operation
Cubature.jl copied to clipboard

Suggest a example, which is easy to understand, of function 'hcubature'.

Open a2468834 opened this issue 5 years ago • 0 comments

The integrand in the example now is written in the form of inline function.

But, I believe it's easier to understand how to use 'hcubature' when the integrand is written in "caller" style. Because it's more often using multi-dimensional integrals of more complex function than f(x, y)=(x^3)*y.

I mean that maybe the example of 'hcubature' could use the following codes.

# ORIGINAL example
hcubature(x -> begin println(x[1],",",x[2]); x[1]^3*x[2]; end, [0,0],[1,1])
# SUGGESTION example
function f( x::Vector{Float64} ) # the integrand i.e., f(x, y)=x^3 y
    println(x[1], ",", x[2])
    return (x[1]^3) * x[2]
end
hcubature(f, [0,0],[1,1])

a2468834 avatar May 01 '19 14:05 a2468834