Re: Issue 3148 in sympy: Too many constants from dsolve()

2013-05-27 Thread sympy
Comment #11 on issue 3148 by smi...@gmail.com: Too many constants from dsolve() http://code.google.com/p/sympy/issues/detail?id=3148 A simpler example of renumbering doing the wrong thing (?) is constant_renumber(C1/(C1*x + 1), C, 1, 2) - C1/(C2*x + 1) -- You received this message because

Issue 3851 in sympy: separable_reduced dsolve solution incorrect?

2013-05-27 Thread sympy
Status: Valid Owner: Labels: Type-Defect Priority-Medium New issue 3851 by smi...@gmail.com: separable_reduced dsolve solution incorrect? http://code.google.com/p/sympy/issues/detail?id=3851 dsolve(eq, hint='separable_reduced') log(x*f(x)) - log(x*f(x) - 1) == C1 + log(x) sol=_

Re: Issue 3851 in sympy: separable_reduced dsolve solution incorrect?

2013-05-27 Thread sympy
Comment #1 on issue 3851 by smi...@gmail.com: separable_reduced dsolve solution incorrect? http://code.google.com/p/sympy/issues/detail?id=3851 In the above, eq = sin(x)*cos(f(x)) + cos(x)*sin(f(x))*f(x).diff(x) -- You received this message because this project is configured to send all

Re: Issue 3148 in sympy: Too many constants from dsolve()

2013-05-27 Thread sympy
Comment #12 on issue 3148 by smi...@gmail.com: Too many constants from dsolve() http://code.google.com/p/sympy/issues/detail?id=3148 And here's another simpler case where the numbering issue can be seen: Here is a differential equation eq x**2*f(x)**2 + x*Derivative(f(x), x) The

Issue 3852 in sympy: Multiple integrals with piecewise should be smarter

2013-05-27 Thread sympy
Status: Valid Owner: CC: julien.r...@gmail.com Labels: Type-Defect Priority-Medium Integration New issue 3852 by asmeu...@gmail.com: Multiple integrals with piecewise should be smarter http://code.google.com/p/sympy/issues/detail?id=3852 integrate(y*cos(x*y), x, y) ⌠ ⎮ ⎧

Re: Issue 3852 in sympy: Multiple integrals with piecewise should be smarter

2013-05-27 Thread sympy
Comment #1 on issue 3852 by asmeu...@gmail.com: Multiple integrals with piecewise should be smarter http://code.google.com/p/sympy/issues/detail?id=3852 I guess I can just use a polynomial for my tutorial. -- You received this message because this project is configured to send all issue

Issue 3853 in sympy: recognize elliptical integrals

2013-05-27 Thread sympy
Status: Valid Owner: Labels: Type-Defect Priority-Medium New issue 3853 by smi...@gmail.com: recognize elliptical integrals http://code.google.com/p/sympy/issues/detail?id=3853 This requires about 2 minutes Ellipse((0,0),3,1).circumference.n() 13.3648932205553 This is nearly

Re: Issue 3148 in sympy: Too many constants from dsolve()

2013-05-27 Thread sympy
Comment #13 on issue 3148 by asmeu...@gmail.com: Too many constants from dsolve() http://code.google.com/p/sympy/issues/detail?id=3148 Exactly. The bug here is in constantsimp, not checkodesol, which cannot help it if constantsimp returns what is basically a wrong result. -- You

Re: Issue 3851 in sympy: separable_reduced dsolve solution incorrect?

2013-05-27 Thread sympy
Updates: Cc: manojkum...@gmail.com Labels: Solvers Comment #2 on issue 3851 by asmeu...@gmail.com: separable_reduced dsolve solution incorrect? http://code.google.com/p/sympy/issues/detail?id=3851 (No comment was entered for this change.) -- You received this message because

Re: Issue 3851 in sympy: separable_reduced dsolve solution incorrect?

2013-05-27 Thread sympy
Comment #3 on issue 3851 by asmeu...@gmail.com: separable_reduced dsolve solution incorrect? http://code.google.com/p/sympy/issues/detail?id=3851 I didn't know about _random. That's a useful method. We should give it a better name and expose it to the public API (remove the leading

Re: Issue 3853 in sympy: recognize elliptical integrals

2013-05-27 Thread sympy
Updates: Labels: Geometry Comment #2 on issue 3853 by asmeu...@gmail.com: recognize elliptical integrals http://code.google.com/p/sympy/issues/detail?id=3853 (No comment was entered for this change.) -- You received this message because this project is configured to send all issue

Re: Issue 3853 in sympy: recognize elliptical integrals

2013-05-27 Thread sympy
Updates: Labels: Integration Comment #1 on issue 3853 by asmeu...@gmail.com: recognize elliptical integrals http://code.google.com/p/sympy/issues/detail?id=3853 What is the actual integral being computed? -- You received this message because this project is configured to send all

Re: Issue 1769 in sympy: cos(oo) should return nan

2013-05-27 Thread sympy
Comment #9 on issue 1769 by asmeu...@gmail.com: cos(oo) should return nan http://code.google.com/p/sympy/issues/detail?id=1769 This comes up with a bunch of functions, actually. If evaluation at oo is not specifically defined in eval(), it is left alone. This leads to a lot of wrong

Re: Issue 1769 in sympy: cos(oo) should return nan

2013-05-27 Thread sympy
Comment #8 on issue 1769 by asmeu...@gmail.com: cos(oo) should return nan http://code.google.com/p/sympy/issues/detail?id=1769 One issue with sin(oo) is that it is not treated like nan in the core: In [74]: sin(oo)/oo Out[74]: 0 (note that limit(sin(x)/x, x, oo) is 1). This is because

Re: Issue 3851 in sympy: separable_reduced dsolve solution incorrect?

2013-05-27 Thread sympy
Updates: Blockedon: sympy:3148 Comment #4 on issue 3851 by asmeu...@gmail.com: separable_reduced dsolve solution incorrect? http://code.google.com/p/sympy/issues/detail?id=3851 I guess the source of this problem is described in issue 3148. -- You received this message because this

Re: Issue 3148 in sympy: Too many constants from dsolve()

2013-05-27 Thread sympy
Issue 3148: Too many constants from dsolve() http://code.google.com/p/sympy/issues/detail?id=3148 This issue is now blocking issue sympy:3851. See http://code.google.com/p/sympy/issues/detail?id=3851 -- You received this message because you are listed in the owner or CC fields of this issue, or

Re: Issue 3434 in sympy: warning when starting isympy

2013-05-27 Thread sympy
Updates: Status: Fixed Comment #14 on issue 3434 by asmeu...@gmail.com: warning when starting isympy http://code.google.com/p/sympy/issues/detail?id=3434 https://github.com/sympy/sympy/pull/2101 -- You received this message because this project is configured to send all issue

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Comment #11 on issue 3427 by asmeu...@gmail.com: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 https://github.com/sympy/sympy/pull/2101 was merged, but was only a partial fix. -- You received this message because this project is configured to

Re: Issue 3761 in sympy: incorrect integral evaluation

2013-05-27 Thread sympy
Updates: Status: Fixed Comment #4 on issue 3761 by julien.r...@gmail.com: incorrect integral evaluation http://code.google.com/p/sympy/issues/detail?id=3761 (No comment was entered for this change.) -- You received this message because this project is configured to send all issue

Re: Issue 3852 in sympy: Multiple integrals with piecewise should be smarter

2013-05-27 Thread sympy
Comment #4 on issue 3852 by asmeu...@gmail.com: Multiple integrals with piecewise should be smarter http://code.google.com/p/sympy/issues/detail?id=3852 I'm unsure about reducing that piecewise, especially if we adopt the programming type convention that was preferred back in

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Comment #12 on issue 3427 by hacm...@gmail.com: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 In order to get an objective test, it is probably a good idea to look at the import time after clearing the disk cache. Under Linux, this can be done

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Comment #13 on issue 3427 by hacm...@gmail.com: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 I guess I'm suggesting that this may be a dupe of https://code.google.com/p/sympy/issues/detail?id=1291 -- You received this message because this

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Comment #14 on issue 3427 by asmeu...@gmail.com: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 If you import SymPy in a clean environment (i.e., a virtualenv or something similar), and compare it against your normal environment, you can see how

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Issue 3427: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 This issue is now blocking issue sympy:1291. See http://code.google.com/p/sympy/issues/detail?id=1291 -- You received this message because you are listed in the owner or CC fields of this

Re: Issue 1291 in sympy: reduce the import sympy time

2013-05-27 Thread sympy
Updates: Blockedon: sympy:3427 Comment #27 on issue 1291 by asmeu...@gmail.com: reduce the import sympy time http://code.google.com/p/sympy/issues/detail?id=1291 (No comment was entered for this change.) -- You received this message because this project is configured to send all

Re: Issue 3427 in sympy: sympy imports matplotlib on start-up

2013-05-27 Thread sympy
Comment #16 on issue 3427 by asmeu...@gmail.com: sympy imports matplotlib on start-up http://code.google.com/p/sympy/issues/detail?id=3427 There's already some good discussion here, so I'll just mark this issue as blocking that one. -- You received this message because this project is

[sympy] Functional equations

2013-05-27 Thread F. B.
Generic Partial Differential Equations may yield arbitrary functions in their solutions. When matching this generic solutions to initial or boundary conditions, we get a functional equation: that is an equation whose variable is a function (without derivatives). If the function to be found

[sympy] Merge dsolve into solve?

2013-05-27 Thread F. B.
Hi, what about putting a dispatcher into solve to redirect to dsolve in case it detects the system contains a differential equation? I mean something like this: f = Function('f') solve(f(x).diff(x) - f(x), f(x)) # solve as a non-differential equation [f(x)] solve(f(x).diff(x) - f(x), f)

Re: [sympy] Merge dsolve into solve?

2013-05-27 Thread Aaron Meurer
This has been discussed before. We want both the behaviors of solve and dsolve. i.e., we want both solve(f(x).diff(x) - f(x), f(x)) [Derivative(f(x), x)] and dsolve(f(x).diff(x) - f(x), f(x)) f(x) == C1*exp(x) I don't like the idea of using f instead of f(x) as a flag. What happens if

Re: [sympy] Functional equations

2013-05-27 Thread Aaron Meurer
On Mon, May 27, 2013 at 12:14 PM, F. B. franz.bona...@gmail.com wrote: Generic Partial Differential Equations may yield arbitrary functions in their solutions. When matching this generic solutions to initial or boundary conditions, we get a functional equation: that is an equation whose

[sympy] Undoable integral example

2013-05-27 Thread Aaron Meurer
What is a good example of a (preferably simple) integral that SymPy will not likely be able to ever do, because there really aren't any closed forms of it, even in terms of special functions? I need a nice example of when integrate() returns an Integral() in my new tutorial. Either definite or

Re: [sympy] Undoable integral example

2013-05-27 Thread Chris Smith
The circumference of an ellipse has no closed form solution, I believe. -- 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 sympy+unsubscr...@googlegroups.com. To post to

Re: [sympy] Undoable integral example

2013-05-27 Thread Chris Smith
see docstring of Ellipse.circumference for integrand -- 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 sympy+unsubscr...@googlegroups.com. To post to this group, send

Re: [sympy] Undoable integral example

2013-05-27 Thread Tom Bachmann
But that's essentially by definition a complete elliptic integral ... I think any integral which we would describe as simple is simple enough to have attracted someone's attention and been given a special name. I might be wrong, of course. On 27.05.2013 19:06, Chris Smith wrote: The

[sympy] Re: Undoable integral example

2013-05-27 Thread Rathmann
How do you want to define closed form? If you allow Meijer G-functions, isn't pretty much anything integable? There are lots of examples expression without elementary antiderivatives: sin(x)/x , e**(-x**2), etc., but Sympy gives answers for these, of course. On Monday, May 27, 2013 10:53:22

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread Aaron Meurer
Closed-form for me means SymPy returns an answer. simple for me means not too complicated, and preferably not using any special functions. I just want something like If integrate() cannot compute the integral, it returns an unevaluated Integral() object example here for the new tutorial, but

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread Chris Smith
I don't think anybody will have this one as it is of little use according to http://www.sciforums.com/showthread.php?76644-Integrate-X-X -- 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,

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread Chris Smith
I don't get that from master. -- 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 sympy+unsubscr...@googlegroups.com. To post to this group, send email to

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread Aaron Meurer
That's a good idea. WolframAlpha doesn't return a solution. Maybe someone could try it in Axiom and elsewhere (what other system besides Mathematica is good with special functions?). Aaron Meurer On Mon, May 27, 2013 at 1:38 PM, Chris Smith smi...@gmail.com wrote: How about x**x? -- You

[sympy] Re: Undoable integral example

2013-05-27 Thread Rathmann
But Sympy returns a result for x**x as well. It is just not in terms that most will find easy to understand: integrate(x**x) Piecewise((x*x**x*gamma(x + 1)/gamma(x + 2), Abs(x) 1), (x*x**x*gamma(x + 1)/gamma(x + 2) + gamma(x + 1)/gamma(x + 2) + gamma(-x - 1)/gamma(-x), Abs(1/x) 1),

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread someone
I don't get that from master. It is wrong, check for example the first branch of Piecewise: In [210]: Piecewise((x*x**x*gamma(x + 1)/gamma(x + 2), Abs(x) 1), (x*x**x*gamma(x + 1)/gamma(x + 2) + gamma(x + 1)/gamma(x + 2) + gamma(-x - 1)/gamma(-x), Abs(1/x) 1), (meijerg(((1,), (x + 2,)), ((x

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread Aaron Meurer
That was before the fix from this issue: https://code.google.com/p/sympy/issues/detail?id=3761, which was a wrong result. So very likely that answer is wrong. Aaron Meurer On Mon, May 27, 2013 at 1:45 PM, Rathmann rathmann...@gmail.com wrote: But Sympy returns a result for x**x as well. It is

Re: [sympy] Re: Undoable integral example

2013-05-27 Thread someone
Maybe someone could try it in Axiom This is from Fricas compiled today: (1) - integrate(x^x, x) x ++%A (1) | %A d%A ++ Type: Union(Expression(Integer),...) -- You received this message because you are

[sympy] Google Summer of Code Students Announced

2013-05-27 Thread Aaron Meurer
Hi everyone. As many of you may have noticed, Google has announced the results for Google Summer of Code. I am proud to announce that we got six slots from Google. The following projects have been accepted: Student (Project): Mentor - Chetna Gupta (Risch algorithm for symbolic integration):

[sympy] Re: Undoable integral example

2013-05-27 Thread Rathmann
Oops, sorry. Apparently I haven't been updating my branch properly. On Monday, May 27, 2013 10:53:22 AM UTC-7, Aaron Meurer wrote: What is a good example of a (preferably simple) integral that SymPy will not likely be able to ever do, because there really aren't any closed forms of it,

[sympy] Linear Algebra Module

2013-05-27 Thread Saurabh Jha
Hi, I had proposed a proposal this year on Linear Algebra Module[1]. This proposal is of high priority in sympy and I want to implement it. Before starting head on to implement it, I wanted to discuss with community about the soundness of the proposed architecture. Any suggestions for

Re: [sympy] Functional equations

2013-05-27 Thread F. B.
Wolfram Mathematica has RSolve: http://reference.wolfram.com/mathematica/ref/RSolve.html It is mainly used to solve recurrence equations, though it is able to accept functional equations too. Wikipedia on recurrence equations: http://en.wikipedia.org/wiki/Recurrence_relation On Monday, May

Re: [sympy] Functional equations

2013-05-27 Thread Aaron Meurer
SymPy also has rsolve(), but it only solves recurrence relations. Aaron Meurer On May 27, 2013, at 4:23 PM, F. B. franz.bona...@gmail.com wrote: Wolfram Mathematica has RSolve: http://reference.wolfram.com/mathematica/ref/RSolve.html It is mainly used to solve recurrence equations, though it

Re: [sympy] Linear Algebra Module

2013-05-27 Thread Aaron Meurer
I think it's fine enough. It's best to just start with some pull requests, and we can work out the details of what should be what way from there. Start out small so that you don't waste much work if you change things, but once you get on track, you can do more. Aaron Meurer On Mon, May 27, 2013