From: SourceForge.net <no...@so...> - 2012-09-19 16:22:15
|
Bugs item #3569562, was opened at 2012-09-19 09:22 Message generated for change (Tracker Item Submitted) made by You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core - Simplification Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: taylor(elliptic_kc(m), m, 0, 1) fails Initial Comment: Maxima cannot compute the Taylor expansion of elliptic_kc(m) for m = 0. Expansions about other points is ok. I think the problem is that diff(elliptic_kc(m),m) at m = 0 produces a form that is indeterminate at 0. Maxima is unable to evaluate the limit as m approaches 0. ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 |
From: SourceForge.net <no...@so...> - 2012-09-20 07:40:01
|
Bugs item #3569562, was opened at 2012-09-19 09:22 Message generated for change (Comment added) made by aleksasd You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 Please note that this message will contain a full copy of the comment thread, including the initial issue submission, for this request, not just the latest update. Category: Lisp Core - Simplification Group: None Status: Open Resolution: None Priority: 5 Private: No Submitted By: https://www.google.com/accounts () Assigned to: Nobody/Anonymous (nobody) Summary: taylor(elliptic_kc(m), m, 0, 1) fails Initial Comment: Maxima cannot compute the Taylor expansion of elliptic_kc(m) for m = 0. Expansions about other points is ok. I think the problem is that diff(elliptic_kc(m),m) at m = 0 produces a form that is indeterminate at 0. Maxima is unable to evaluate the limit as m approaches 0. ---------------------------------------------------------------------- Comment By: Aleksas (aleksasd) Date: 2012-09-20 00:40 Message: From maxima help: Function: elliptic_kc (m) The complete elliptic integral of the first kind, defined as integrate(1/sqrt(1 - m*sin(x)^2), x, 0, %pi/2) Example 1 Compute taylor(elliptic_kc(m), m, 0, 1) (%i1) taylor(integrate(1/sqrt(1-m*sin(x)^2),x,0,%pi/2),m,0,1); (%o1)/T/ %pi/2+((at(integrate(sin(x)^2,x,0,%pi/2),m=0))*m)/2+... (%i2) ev(%,integrate); (%o2)/R/ (%pi*m+4*%pi)/8 (%i3) taylor(%,m,0,1); (%o3)/T/ %pi/2+(%pi*m)/8+... Example 2 Compute taylor(elliptic_kc(m), m, 0, 3) (%i4) taylor(integrate(1/sqrt(1-m*sin(x)^2),x,0,%pi/2),m,0,3)$ (%i5) ev(%,nouns)$ (%i6) taylor(%,m,0,3); (%o6)/T/ %pi/2+(%pi*m)/8+(9*%pi*m^2)/128+(25*%pi*m^3)/512+... best Aleksas D ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3569562&group_id=4933 |