JuliaMath / Cubature.jl

One- and multi-dimensional adaptive integration routines for the Julia language
Other
126 stars 21 forks source link

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

Open a2468834 opened 5 years ago

a2468834 commented 5 years ago

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])