SymbolicNumericIntegration.jl

SymbolicNumericIntegration.jl copied to clipboard

julia> using Symbolics, SymbolicNumericIntegration
julia> @variables α;

Noticing that the case

julia> integrate(2 * α, α)
(α^2, 0, 0)

is correct, going element-wise works essentially as intended

julia> exp2 = [1, 2 * α];
julia> integrate.(exp2, α)
2-element Vector{Tuple{Any, Int64, Int64}}:
 (α, 0, 0)
 (α^2, 0, 0)

but otherwise just integrating exp2 on its own is problematic, yielding

julia> integrate(exp2, α)
(Num[α, 2(α^2)], 0, 0)

with the expectation of [α, α^2 / 2]. It seems that the vector is treated as constant wrt. α and α * [1, 2 * α] is just being simplified behind the scenes.

Since this is not yet supported, should this be explicitly disallowed/ given with a warning, etc?