MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

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?