Re: Re: Weird result in Mathematica 6
- To: mathgroup at smc.vnet.net
- Subject: [mg76566] Re: [mg76431] Re: [mg76393] Weird result in Mathematica 6
- From: DrMajorBob <drmajorbob at bigfoot.com>
- Date: Thu, 24 May 2007 05:59:05 -0400 (EDT)
- References: <26727995.1179743868970.JavaMail.root@m35> <200705220647.CAA19795@smc.vnet.net> <9745176.1179921476358.JavaMail.root@m35>
- Reply-to: drmajorbob at bigfoot.com
Interesting. But your results are entirely different from mine, for the same input. > which explains what is wrong (error messages can tell you a lot, > sometimes). And frequently, they don't. Here's the error message at THIS machine: FindRoot::lstol: The line search decreased the step size to within \ tolerance specified by AccuracyGoal and PrecisionGoal but was unable \ to find a sufficient decrease in the merit function. You may need \ more than MachinePrecision digits of working precision to meet these \ tolerances. >> Bobby On Tue, 22 May 2007 06:28:08 -0500, Andrzej Kozlowski <akoz at mimuw.edu.pl> wrote: > *This message was transferred with a trial version of CommuniGate(tm) > Pro* > I don't see any connection between these two issues. Moreover, I get: > > FindRoot[h == g, {x, 0}] > FindRoot::njnum:The Jacobian is not a matrix of numbers at {x} = {0.}. >> > {x -> 0.} > > which explains what is wrong (error messages can tell you a lot, > sometimes). Trying a slightly different starting search point: > > FindRoot[h == g, {x, 0.1}] > {x->2.} > > {g, h} /. % > {0., 0.} > > > Andrzej Kozlowski > > > On 22 May 2007, at 15:47, DrMajorBob wrote: > >> Even worse, FindRoot returns a wrong answer: >> >> g = Piecewise[{{0, x < 0}, {2 x - x^2, 0 <= x < 4}, {16 x - x^2, >> x³4}}]; >> h = x - 2; >> FindRoot[h == g, {x, 0}] >> {g, h} /. % >> >> {x->-2.84217*10^-15} >> {0, -2.} >> >> Bobby >> >> On Mon, 21 May 2007 05:01:21 -0500, Sebastian Meznaric >> <meznaric at gmail.com> wrote: >> >>> I was playing around with Mathematica 6 a bit and ran this command to >>> solve for the inverse of the Moebius transformation >>> >>> FullSimplify[ >>> Reduce[(z - a)/(1 - a\[Conjugate] z) == w && a a\[Conjugate] < 1 && >>> w w\[Conjugate] < 1, z]] >>> >>> This is what I got as a result: >>> -1 < w < 1 && -1 < a < 1 && z == (a + w)/(1 + w Conjugate[a]) >>> >>> Why is Mathematica assuming a and w are real? The Moebius >>> transformation is invertible in the unit disc regardless of whether a >>> and w are real or not. Any thoughts? >>> >>> >>> >> >> >> >> --DrMajorBob at bigfoot.com >> > > -- DrMajorBob at bigfoot.com
- References:
- Re: Weird result in Mathematica 6
- From: DrMajorBob <drmajorbob@bigfoot.com>
- Re: Weird result in Mathematica 6