Re: using FindRoot to find multiple answers in a domain?
- To: mathgroup at smc.vnet.net
- Subject: [mg69261] Re: using FindRoot to find multiple answers in a domain?
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Mon, 4 Sep 2006 02:41:25 -0400 (EDT)
- Organization: The University of Western Australia
- References: <firstname.lastname@example.org> <email@example.com>
- Sender: owner-wri-mathgroup at wolfram.com
In article <ed3rou$rv5$1 at smc.vnet.net>, "akil" <akomur at wanadoo.nl>
> IntervalRoots package and RootSearch do not work for the kind of formulas I
Are you sure?
> For example in the domain [0,Pi/2] we want to solve (I only used FindRoot
> because this was the fastest):
> I put the package I test in on
> http://home.wanadoo.nl/akomur/testRootSearch.nb so that you can see the
> formulas I use.
> Copying them In here ruins the formulas.
Copy as Plain text. Your first formula is
nlv = Max[-81.24275115593154 Cot[beta] - 48.489352280270914,
Min[-46.844746272761526 Cot[beta] - 71.4213555357176,
-46.480384751324436 Cot[beta] - 71.66426321667566]]
Note that PiecewiseExpand can be applied to such formulae. If you
simplify the result
you deduce that nlv is just
-81.24275115593154 Cot[beta] - 48.489352280270914 if Cot[beta] <= 2/3
-46.844746272761526 Cot[beta] - 71.4213555357176 if Cot[beta] >= 2/3
You can check that these values are consistent when Cot[beta] == 2/3.
Moreover, the catch-all condition
-46.480384751324436 Cot[beta] - 71.66426321667566
is redundant for real beta since it only applies when Cot[beta] == 2/3.
A brief examination of nuvb shows that it involves ratios of (powers of)
expressions similar to nlv. I expect that simplification of these
expressions prior to constructing nuvb would lead to a final form that
is not that hard to handle.
I expect that if you step back and tell us more about the original
problem, then a more elegant approach would present itself.
Paul Abbott Phone: 61 8 6488 2734
School of Physics, M013 Fax: +61 8 6488 1014
The University of Western Australia (CRICOS Provider No 00126G)
Prev by Date:
Re: Co-Displaying Combinatorica and Graphics Objects
Next by Date:
Re: comboboxes and paste to notebook
Previous by thread:
Re: scaled complementary error function in Mathematica?
Next by thread:
Re: Trigonometric simplification