On Thu, 11 Mar 2021 at 18:54, Rohan Gupta <[email protected]> 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 -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/CAHVvXxQKPoO7LH1Mg%3DzWKWN_jbG2GjQm0y%3DV3%2B1b1WCX2i9GYg%40mail.gmail.com.
