Opened 11 years ago
Last modified 4 years ago
#2787 new enhancement
make the interface to integrate() like the (new consistent) interface to diff()
Reported by: | jason | Owned by: | was |
---|---|---|---|
Priority: | major | Milestone: | sage-6.4 |
Component: | calculus | Keywords: | |
Cc: | kcrisman, ktkohl | Merged in: | |
Authors: | Reviewers: | ||
Report Upstream: | N/A | Work issues: | |
Branch: | Commit: | ||
Dependencies: | Stopgaps: |
Description (last modified by )
It would be nice if the following worked, if f was a function:
sage: integrate(f, x) sage: # Double indefinite integral wrt x sage: integrate(f, x, x) sage: # limits and a double integral: x is the inner integral, y is the outer integral (note that this order is backwards from Mma...they think of nested integrals as int dx int dy function, so the first parameter is the outer integral in Mma.) sage: integrate(f, (x, 0, 1), (y, 1, 2)) sage: integrate(f, (x, 0, y), (y, 1, 2)) sage: # Double integral, x is the inner integral, y is the outer integral sage: integrate(f, x, y)
Change History (13)
comment:1 Changed 10 years ago by
comment:2 Changed 10 years ago by
What do you propose to do about the following case:
# Current behavior, definite integral with x and y as bounds sage: var('x,y') sage: integrate(x^2,x,x,y) y^3/3 - x^3/3 # Triple Integral, twice with respect to x, once with respect to y sage: integrate(x^2,x,x,y) x^4*y/12 # i.e. same as sage: (x^2).integrate(x).integrate(x).integrate(y)
comment:3 Changed 10 years ago by
As noted in #1221, I would vote with the consensus from the devel list that the syntax integrate(f, x, a, b) be deprecated in favor of integrate(f, (x, a, b)). However, you're right that this would introduce a backwards-incompatible change.
comment:4 Changed 8 years ago by
- Cc kcrisman added
- Priority changed from major to minor
- Report Upstream set to N/A
comment:5 Changed 6 years ago by
- Milestone changed from sage-5.11 to sage-5.12
comment:6 Changed 5 years ago by
- Milestone changed from sage-6.1 to sage-6.2
comment:7 Changed 5 years ago by
- Milestone changed from sage-6.2 to sage-6.3
comment:8 Changed 5 years ago by
- Milestone changed from sage-6.3 to sage-6.4
comment:9 Changed 4 years ago by
- Description modified (diff)
- Priority changed from minor to major
comment:10 Changed 4 years ago by
- Description modified (diff)
comment:11 Changed 4 years ago by
- Cc ktkohl added
comment:12 Changed 4 years ago by
I would say that deprecating integral(f(x),x,a,b)
for a definite integral will be too large a backward incompatible change to pull of.
The gain for multiple integrals is also a little tricky: The proposed order for integral(f(x,y),x,y) where x is inner and y is outer is definitely consistent with mathematical notation int int f(x,y) dx dy
but it's opposite to python's iteration order [ (i,j) for i in range(3) for j in range(3) ]
, where j iterates in the inner loop.
Also, the proposed notation is for multiple integrals, which in mathematics are usually written as multiple (nested) integrals anyway. Single integral notation is usually reserved for higher dimensional integrals (Lebesgue, usually), in which case no particular "integration order" is implies, and indeed the domain of integration might not allow for a direct translation into nested integrals.
The fact that integral has to both accommodate definite integrals and taking antiderivatives breaks the complete symmetry with the diff
interface. It would be nice if the interface were more symmetrical, but I think there are good reasons why it's not. If a complete converse to diff()
is required, perhaps we need antiderivative()
, which doesn't need to accommodate definite integrals (although of course no-one would find its existence).
comment:13 Changed 4 years ago by
I think one big meta-motivation behind this is that some competitor programs may be more consistent or at least allow multiple integrals this way? Certainly it seems much more cumbersome in Sage - plus, you have to do them in the "right order" to get nice answers sometimes, I would guess, just like in calculus textbooks :)
See also #1221. It would be nice if the interface here was also consistent with the interface to plot, at least when a range is specified, since both take a function, a variable, and a range.