Dropping Higher Order Terms (from an earlier post)

*To*: mathgroup at smc.vnet.net*Subject*: [mg45996] Dropping Higher Order Terms (from an earlier post)*From*: Harold.Noffke at wpafb.af.mil (Harold Noffke)*Date*: Mon, 2 Feb 2004 05:20:33 -0500 (EST)*Sender*: owner-wri-mathgroup at wolfram.com

I have been studying David Park's solution (second one in this post), but always receive a "Take" error of some sort. On the other hand, In[1]/Out[1] immediately below provides the answer I would expect to see. . . . . . . . . . . . . . . . . . . . . In[1]:= (Series[f[x]/g[x], {x, 0, 3}] // Normal)*g[x] f[0] f'[0] f[0] g'[0] Out[1]= g[x] (---- + x (----- - ----------) + g[0] g[0] 2 g[0] 2 2 f'[0] g'[0] f[0] g'[0] f''[0] > x (-(-----------) + ----------- + ------ - 2 3 2 g[0] g[0] g[0] f[0] g''[0] > -----------) + 2 2 g[0] 2 3 3 f'[0] g'[0] f[0] g'[0] g'[0] f''[0] > x (------------ - ----------- - ------------ - 3 4 2 g[0] g[0] 2 g[0] (3) f'[0] g''[0] f[0] g'[0] g''[0] f [0] > ------------ + ----------------- + ------- - 2 3 6 g[0] 2 g[0] g[0] (3) f[0] g [0] > ------------)) 2 6 g[0] . . . . . . . . . . . . . . . . . . . . David -- I do not understand why you took the approach immediately below, as I cannot get it's solution past the point of producing "Take" errors. Begin David Park's solution ... "Let's assume that you don't want to disturb any higher order terms that might be contained in g[x]. Then something like the following would work." f[x_] = Array[a, 10, 0] . Array[x^#1 & , 10, 0]*g[x] (a[0] + x*a[1] + x^2*a[2] + x^3*a[3] + x^4*a[4] + x^5*a[5] + x^6*a[6] + x^7*a[7] + x^8*a[8] + x^9*a[9])* g[x] MapAt[Take[#, 3] &, f[x], 1] (a[0] + x*a[1] + x^2*a[2])*g[x] . . . . . . . . . . . . . . . . . . . . Regards, Harold