Re: Coercion into series
- To: mathgroup at smc.vnet.net
- Subject: [mg71469] Re: Coercion into series
- From: bghiggins at ucdavis.edu
- Date: Mon, 20 Nov 2006 02:43:50 -0500 (EST)
- References: <ejosoj$n41$1@smc.vnet.net>
Dimitris, You can write your equation as -3 - a^4 + b*c + (b + c)*y=0 Now a series expansion of the above equation about y=0 gives the above equation for O[x]^1 Thus it is clear that Mathematica has found the solution for the special case b+c=0 Why Solve does this is not clear to me. Though it is consistent as Solve[b*c + (b + c)*y^4 - (a^4 + 3) == 0 + O[y]^5, {b, c}] gives the same solution. Whereas, Solve[b*c + (b + c)*y^4 - (a^4 + 3) == 0 + O[y]^4, {b, c}] gives {{b -> (3 + a^4)/c}} Note that the series expansion of the solution Normal[Series[b /. Solve[3 + a^4 + b*c - (b + c)*y^4 == 0, b], {y, 0, 5}]] is {-(3/c) - a^4/c + (1 - 3/c^2 - a^4/c^2)*y^4} which of course agrees with Solve[b*c + (b + c)*y^4 - (a^4 + 3) == 0 + O[y]^5, {b, c}], if b+c=0 Cheers, Brian dimitris wrote: > Any time that an object like O[x] appears in a sum of terms, > Mathematica will convert the whole sum into a power series. > > E.g. > > a x + Exp[x] + O[x]^3 > SeriesData[x, 0, {1, 1 + a, 1/2}, 0, 3, 1] > > Cos[z] + O[z, Pi/2]^4 > SeriesData[z, Pi/2, {-1, 0, 1/6}, 1, 4, 1] > > and so on... > > Consider the following equation > > eq=(y + b)*(y + c) == a^4 + y^2 + 3; > > Obviously > > Solve[eq, {b, c}] > Solve::svars: Equations may not give solutions for all "solve" > variables. > {{b -> (3 + a^4 - c*y)/(c + y)}} > > However how the appearance of O[y]^3 below "makes" Solve give solution > for {b,c} > > Solve[(y + b)*(y + c) == (a^4 + y^2 + 3) + O[y]^3, {b, c}] > {{b -> -Sqrt[-3 - a^4], c -> Sqrt[-3 - a^4]}, {b -> Sqrt[-3 - a^4], c > -> -Sqrt[-3 - a^4]}} > > Thanks a lot. > > Dimitris