Re: ALL roots of non-polynomial equation
- To: mathgroup at smc.vnet.net
- Subject: [mg36038] Re: ALL roots of non-polynomial equation
- From: Selwyn Hollis <slhollis at earthlink.net>
- Date: Tue, 13 Aug 2002 05:23:01 -0400 (EDT)
- References: <aj7omm$2up$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Looks like we've been reinventing the wheel, folks. The package
NumericalMath`IntervalRoots` has exactly what we need.
<<NumericalMath`IntervalRoots`
In: IntervalNewton[Sin[3*x2] - 0.2*x, x, Interval[{0, 2}],
10^(-12), MaxRecursion -> 20]
Out: Interval[{0, 1.337913459004121`*^-18}, {0.06666864225044185`,
0.06666864225044247`}, {0.9903214888289147`,
0.9903214888289159`}, {1.4814395017188728`,
1.4814395017188753`}, {1.7387361933371692`, 1.7387361933371717`}]
In: First /@ List @@ %
Out: {0, 0.0666686, 0.990321, 1.48144, 1.73874}
---------
Selwyn Hollis