[Date Index]
[Thread Index]
[Author Index]
Re: problem with Expand
*To*: mathgroup at smc.vnet.net
*Subject*: [mg26259] Re: problem with Expand
*From*: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
*Date*: Sun, 10 Dec 2000 00:19:34 -0500 (EST)
*Organization*: Universitaet Leipzig
*References*: <90kqsc$r3g@smc.vnet.net>
*Sender*: owner-wri-mathgroup at wolfram.com
Hi,
a bit of basic algebra
-a - b == -(a+b)
and
Expand[pn1 - pn - k/2*fn*fn]
is
Expand[pn1 - (pn + k/2*fn*fn)]
Expand[pn1]- Expand[pn+k/2*fn*fn]
Don't use a computer algebra before you have
finished the primary school.
Regards
Jens
Peter Pollner wrote:
>
> hi,
>
> I got a strange answer from Expand in the following situation:
> fn1, pn1, fn and pn are polynomials upto order to in x. I want to check
> the identity:
>
> pn1 - pn - k/2*fn*fn = 0 + O(x^3)
>
> the answer depends on the usage of Expand:
> Expand[pn1 - pn - k/2*fn*fn]
> and
> Expand[pn1] - Expand[pn - k/2*fn*fn]
> are not the same (the first expansion contains O(x^2) terms)!
>
> Here is the full input:
>
> pn1 = p0*x - k/12*n(n + 1)(2n + 1)x^2p0^2 - k/2*n(n + 1)p0*f0*x^2 -
> k/2*(n + 1)x^2f0^2
>
> fn1 = f0*x + (n + 1)p0*x - k/24*n(n + 1)(n + 1)(n + 2)x^2 p0^2 -
> k/6*n(n + 1)(n + 2)x^2p0*f0 - k/4*(n + 1)(n + 2)x^2f0^2
>
> pn = p0*x - k/12*(n - 1)n(2(n - 1) + 1)x^2p0^2 - k/2*n(n - 1)p0*f0*x^2 -
> k/2*n x^2f0^2
>
> fn = f0*x + n p0*x - k/24*(n - 1)n^2(n + 1)x^2 p0^2 -
> k/6*(n - 1)n(n + 1)x^2p0*f0 - k/4*n(n + 1)x^2f0^2
>
> Expand[pn1 - pn - k/2*fn*fn]
>
> Expand[pn1] - Expand[pn - k/2*fn*fn]
>
> I am using Mathematica 4.0.1.0, Platform : X
>
> Is it a bug?
>
> Please send me answers to the email address:
> Peter.Pollner at Physik.Uni-Marburg.De
>
> Thanks
> Peter Pollner
Prev by Date:
**RE: Easy way to simplify expr to a+bi?**
Next by Date:
**Graphics3D[Plot3D[]] vs. ParametricPlot3D[]**
Previous by thread:
**problem with Expand**
Next by thread:
**SIMPLIFY: Results seem to depend on the form of the expression**
| |