Path: csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!gegeweb.org!eternal-september.org!feeder.eternal-september.org!.POSTED!not-for-mail From: Joe Riel Newsgroups: comp.soft-sys.math.maple Subject: Re: Complex Roots Date: Fri, 01 Jul 2011 10:39:21 -0700 Organization: A noiseless patient Spider Lines: 31 Message-ID: <87sjqp294m.fsf@san.rr.com> References: Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii Injection-Info: mx04.eternal-september.org; posting-host="7daQ3AF9ALJlnU9jGWSG5Q"; logging-data="7168"; mail-complaints-to="abuse@eternal-september.org"; posting-account="U2FsdGVkX19vPbckLRcVfoieFWk1xFzQ" User-Agent: Gnus/5.13 (Gnus v5.13) Emacs/23.3 (gnu/linux) Cancel-Lock: sha1:cEhYA1RhdfRPnEBI5xEh0cw6ciY= sha1:SAvkUYCeYhx/WOZIfn64EPt1WB4= Xref: x330-a1.tempe.blueboxinc.net comp.soft-sys.math.maple:198 Joel writes: > I am trying to find all the complex roots of the equation T given > below: > > w := 60*(2*3.14159); > R12 := .5; > L12 := 3/w; > Z12 := R12+s*L12+I*w*L12; > kp := 0.5e-3; > kv := 0.5e-3; > Q1 := 384; > Q2 := 375; > > apw1 := s+kp*(-Q1+w*L12/abs(Z12)^2); > apw2 := s+kp*(-Q2+w*L12/abs(Z12)^2); > aqv1 := 1+kv*(Q1+w*L12/abs(Z12)^2); > aqv2 := 1+kv*(Q2+w*L12/abs(Z12)^2); > T := apw1*apw2*aqv1*aqv2; > > > I have used the solve, fsolve (with complex argument), and RootOf > function but all that is returned is one real root or nothing at all. It only has two roots, both real. That is apparent from the form of the expressions. One way to produce both of them is map(fsolve, {op(T)}); -- Joe Riel