On Thu, 11 Mar 2021 at 18:54, Rohan Gupta <
rohan...@gmail.com> wrote:
>
> I'm interested in working on numerical integration techniques for functions.
Is this listed on the ideas page somewhere?
In general numerical techniques are out of scope for sympy and should
be implemented in mpmath or numpy/scipy etc. The idea is that sympy
should use the routines from e.g. mpmath rather than implement its own
numerical algorithms.
There are some algorithms already implemented in mpmath though that
are not used in sympy. For example in sympy you can do:
In [164]: Integral(x, (x, 0, 1)).evalf()
Out[164]: 0.500000000000000
Under the hood this calls mpmath's quad function (I think). This
doesn't work for multiple integrals though:
In [165]: Integral(x*y, (x, 0, 1), (y, 0, 1)).evalf()
Out[165]:
1 1
⌠ ⌠
⎮ ⎮ x⋅y dx dy
⌡ ⌡
0 0
It should be relatively straight-forward to make that work because
mpmath's quad function can handle multiple integrals:
https://mpmath.org/doc/current/calculus/integration.html#standard-quadrature-quad
--
Oscar