Re: A way around the limitations of Re[] and Im[]
- To: mathgroup at smc.vnet.net
- Subject: [mg51118] Re: A way around the limitations of Re[] and Im[]
- From: carlos at colorado.edu (Carlos Felippa)
- Date: Tue, 5 Oct 2004 04:37:04 -0400 (EDT)
- References: <cjr8lg$opp$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
carlos at colorado.edu (Carlos Felippa) wrote in message news:<cjr8lg$opp$1 at smc.vnet.net>...
> As you know, Re[expr] and Im[expr] are left unevaluated when expr is
> not numeric. I had good luck in simple cases with the following
> substitution trick:
>
> Real[c_]:=c/.I->0; Imag[c_]:=(c-Real[c])/.I->1;
>
> Do you see any problem with these definitions?
Not yet there. Found limitations because n*I is internally stored as
Complex[0,n]. This version has done better:
RealPart[c_]:=ComplexExpand[c]/.{Complex[0,_]->0};
Examples:
Print[ RealPart[(a+I*b)^2] //InputForm];
Print[ RealPart[(a+I*b)^n] //InputForm];
Print[ RealPart[Exp[I*x]/I] //InputForm];
Print[ RealPart[Cosh[x+y*I]] //InputForm];
Results:
a^2 - b^2
Abs[a]^n*Cos[n*Arg[a]]
Sin[x]
Cos[y]*Cosh[x]
- Follow-Ups:
- Re: Re: A way around the limitations of Re[] and Im[]
- From: Murray Eisenberg <murray@math.umass.edu>
- Re: Re: A way around the limitations of Re[] and Im[]