From: SourceForge.net <no...@so...> - 2010-05-23 20:44:43
|
Bugs item #3005835, was opened at 2010-05-23 03:16 Message generated for change (Comment added) made by crategus You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3005835&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 - Taylor Group: None >Status: Closed >Resolution: Works For Me Priority: 5 Private: No Submitted By: bh () Assigned to: Nobody/Anonymous (nobody) Summary: taylor Initial Comment: taylor(f(x+h),h,0,2); gives f(x)+(at('diff(f(x+h),h,1),h=0))*h+((at('diff(f(x+h),h,2),h=0))*h^2)/2+... How do I reduce it to the expected result? ev(%,diff,nouns) does not help. ---------------------------------------------------------------------- >Comment By: Dieter Kaiser (crategus) Date: 2010-05-23 22:44 Message: The reported example does not show a bug. The suggested solution to load the package pdiff seems to give the desired results. Closing this bug report as works for me. Dieter Kaiser ---------------------------------------------------------------------- Comment By: bh () Date: 2010-05-23 17:50 Message: Thanks. The pdiff package is nice. It would seem worthy enough to appear in the manual! ---------------------------------------------------------------------- Comment By: Aleksas Domarkas (alex108) Date: 2010-05-23 10:11 Message: Some variants: (%i1) load(pdiff)$ (%i2) taylor(f(x+h),x,0,3); (%o2)/T/ f(h)+f[(1)](h)*x+(f[(2)](h)*x^2)/2+(f[(3)](h)*x^3)/6+... or (%i3) powerseries (f(x+h), x, 0); (%o3) sum((x^i1*(at('diff(f(x+h),x,i1),x=0)))/i1!,i1,0,inf) (%i4) niceindices (%); (%o4) sum((x^i*(at('diff(f(x+h),x,i),x=0)))/i!,i,0,inf) ---------------------------------------------------------------------- You can respond by visiting: https://sourceforge.net/tracker/?func=detail&atid=104933&aid=3005835&group_id=4933 |