MathGroup Archive 2001

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

Search the Archive

Re: Simplify


<Matthias.Bode at oppenheim.de> wrote in message
news:9vmrje$hh7$1 at smc.vnet.net...
> Dear Colleagues,
>
> why should I not expect MATHEMATICA to (Full)Simplify
>
> (a^b)^(1/b) or a^b^(1/b) to a?
>

Matthias,
the reason is because they are not equal for general complex-valued a, b, c.
The following notebook illustrates how to make a rule to enforce the desired
simplification (which is correct only for a, b, c real), and also gives a
case where equality fails.  z^w is defined as Exp[w Log[z]], and Log has a
branch cut; this complication is the reason for the failure.

In[1]:=
test = (a^b)^(1/b)
Out[1]=
\!\(\((a\^b)\)\^\(1\/b\)\)
In[2]:=
FullSimplify[%]
Out[2]=
\!\(\((a\^b)\)\^\(1\/b\)\)
In[3]:=
test2 = (a^b)^c
Out[3]=
\!\(\((a\^b)\)\^c\)
In[4]:=
FullSimplify[%]
Out[4]=
\!\(\((a\^b)\)\^c\)
In[5]:=
rule = Power[Power[a_, b_], c_]\[RuleDelayed] Power[a, b c]
Out[5]=
\!\(\((a_\^b_)\)\^c_ \[RuleDelayed] a\^\(b\ c\)\)
In[6]:=
test /. rule
Out[6]=
a
In[7]:=
test2 /. rule
Out[7]=
\!\(a\^\(b\ c\)\)
In[8]:=
a = 2 + 3 I;b = 3 +I; c = -2-4I;
In[9]:=
N[(a^b)^c]
Out[9]=
\!\(4.240023659792295`*^-7 - 7.75911100859049`*^-7\ \[ImaginaryI]\)
In[10]:=
N[a^(b c)]
Out[10]=
34864.2\[InvisibleSpace]-63800.3 \[ImaginaryI]

Alan



  • Prev by Date: RE: Simplify
  • Next by Date: Re: fit to page using inline postscript
  • Previous by thread: RE: Simplify
  • Next by thread: no math font in exporting plots to .eps files