[Date Index]
[Thread Index]
[Author Index]
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
Prev by Date:
**RE: Notation and Symbolize Packages**
Next by Date:
**Re: clarification. Re: one liner for a function?**
Previous by thread:
**Request from Design Science**
Next by thread:
**Re: clarification. Re: one liner for a function?**
| |