Issue 1163: make log(sin(x+1)+1).limit(x, 0) simpler
http://code.google.com/p/sympy/issues/detail?id=1163
New issue report by ondrej.certik:
It returns the correct result, but a little bit more complex than it should
be:
In [1]: log(sin(x+1)+1).limit(x,0).evalf()
Out[1]: 0.610564700497503
In [2]: log(sin(x+1)+1).subs(x, 0).evalf()
Out[2]: 0.610564700497503
In [3]: log(sin(x+1)+1).limit(x, 0)
Out[3]:
49 15⋅sin(1) 6⋅sin(1) 163⋅sin(1)
- ── - ────────────────────── + ────────── - ──────────────────────────── -
──
20 2 1 + sin(1) 2
1 + 2⋅sin(1) + sin (1) ⎛ 2 ⎞
15
15⋅⎝1 + 2⋅sin(1) + sin (1)⎠
62⋅sin(1)
───────────────────────────────────────────────── -
──────────────────────────
⎛ 2 3 4 ⎞ ⎛ 2
⎞ ⎛
⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ ⎝1 + 2⋅sin(1) + sin
(1)⎠⋅⎝
sin(1) 6
─────────────────────────────────────────────── +
────────────────────────────
2 3 4 ⎞ ⎛
1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ (1 + sin(1))⋅⎝1 +
4⋅sin(1) +
⋅sin(1) 20⋅sin(1)
───────────────────────────────── + ───────────────────────────────────── +
lo
2 3 4 ⎞ ⎛ 2 ⎞
6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ (1 + sin(1))⋅⎝1 + 2⋅sin(1) + sin (1)⎠
6 31
g(1 + sin(1)) + ────────── -
─────────────────────────────────────────────────
1 + sin(1) ⎛ 2 3 4
30⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin
(1
15 1
── - ────────────────────────── -
────────────────────────────────────────────
⎞ ⎛ 2 ⎞ ⎛ 2 ⎞ ⎛
)⎠ 2⋅⎝1 + 2⋅sin(1) + sin (1)⎠ 6⋅⎝1 + 2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1)
+ 6
6
─────────────────────────────── +
────────────────────────────────────────────
2 3 4 ⎞ ⎛ 2
⋅sin (1) + 4⋅sin (1) + sin (1)⎠ 5⋅(1 + sin(1))⋅⎝1 + 4⋅sin(1) + 6⋅sin (1)
+ 4
20
16
─────────────────── + ─────────────────────────────────────── -
──────────────
3 4 ⎞ ⎛ 2 ⎞
⋅sin (1) + sin (1)⎠ 3⋅(1 + sin(1))⋅⎝1 + 2⋅sin(1) + sin (1)⎠ ⎛
60⋅⎝1 +
2⋅sin(
2 2
3 163⋅sin (1) 31⋅sin
(1)
────────────── - ──────────────────────────── -
──────────────────────────────
2 2 ⎛ 2
2 ⎞ ⎛ 2 ⎞ 5⋅⎝1 + 4⋅sin(1) + 6⋅sin (1)
+
1) + sin (1)⎠ 10⋅⎝1 + 2⋅sin(1) + sin (1)⎠
2
15⋅sin (1)
──────────────────── - ────────────────────────── +
──────────────────────────
3 4 ⎞ ⎛ 2 ⎞ ⎛
4⋅sin (1) + sin (1)⎠ 2⋅⎝1 + 2⋅sin(1) + sin (1)⎠ (1 + sin(1))⋅⎝1 +
4⋅sin(1)
2 2
12⋅sin (1) 20⋅sin (1)
─────────────────────────────────── + ─────────────────────────────────────
-
2 3 4 ⎞ ⎛ 2 ⎞
+ 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ (1 + sin(1))⋅⎝1 + 2⋅sin(1) + sin (1)⎠
2
5⋅sin (1)
───────────────────────────────────────────────────────────────────────────
-
⎛ 2 ⎞ ⎛ 2 3 4 ⎞
2⋅⎝1 + 2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠
3 3
163⋅sin (1) 62⋅sin (1)
──────────────────────────── -
───────────────────────────────────────────────
2 ⎛ 2
3 4
⎛ 2 ⎞ 15⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) +
sin
15⋅⎝1 + 2⋅sin(1) + sin (1)⎠
3
12⋅sin (1)
──── + ───────────────────────────────────────────────────────────── -
───────
⎞ ⎛ 2 3 4 ⎞ ⎛
(1)⎠ (1 + sin(1))⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ 3⋅⎝1
+
3
10⋅sin (1)
──────────────────────────────────────────────────────────────────── +
───────
2 ⎞ ⎛ 2 3 4 ⎞
2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ 3⋅(1
+
3 4
20⋅sin (1) 163⋅sin (1)
──────────────────────────────── - ──────────────────────────── -
────────────
⎛ 2 ⎞ 2 ⎛
sin(1))⋅⎝1 + 2⋅sin(1) + sin (1)⎠ ⎛ 2 ⎞ 30⋅⎝1 +
4⋅si
60⋅⎝1 + 2⋅sin(1) + sin (1)⎠
4 4
31⋅sin (1) 6⋅sin
(1)
─────────────────────────────────────── +
────────────────────────────────────
2 3 4 ⎞ ⎛
2
n(1) + 6⋅sin (1) + 4⋅sin (1) + sin (1)⎠ (1 + sin(1))⋅⎝1 + 4⋅sin(1) +
6⋅sin (
4
5⋅sin (1)
───────────────────────── -
──────────────────────────────────────────────────
3 4 ⎞ ⎛ 2 ⎞ ⎛
2
1) + 4⋅sin (1) + sin (1)⎠ 2⋅⎝1 + 2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1) +
6⋅sin (
5
sin (1)
───────────────────────── -
──────────────────────────────────────────────────
3 4 ⎞ ⎛ 2 ⎞ ⎛ 2
1) + 4⋅sin (1) + sin (1)⎠ ⎝1 + 2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1) + 6⋅sin
(1)
5
6⋅sin (1)
─────────────────────── +
────────────────────────────────────────────────────
3 4 ⎞ ⎛ 2 3
+ 4⋅sin (1) + sin (1)⎠ 5⋅(1 + sin(1))⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin
(1)
6
sin (1)
─────────── -
────────────────────────────────────────────────────────────────
4 ⎞ ⎛ 2 ⎞ ⎛ 2 3
+ sin (1)⎠ 6⋅⎝1 + 2⋅sin(1) + sin (1)⎠⋅⎝1 + 4⋅sin(1) + 6⋅sin (1) + 4⋅sin
(1)
───────────
4 ⎞
+ sin (1)⎠
Issue attributes:
Status: Accepted
Owner: ----
Labels: Type-Defect Priority-Medium
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to [EMAIL PROTECTED]
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en
-~----------~----~----~----~------~----~------~--~---
