Re: Bug in Solve?
- To: mathgroup at smc.vnet.net
- Subject: [mg36849] Re: [mg36834] Bug in Solve?
- From: Daniel Lichtblau <danl at wolfram.com>
- Date: Sun, 29 Sep 2002 02:55:11 -0400 (EDT)
- References: <200209280834.EAA25013@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
David wrote:
>
> Hi,
>
> inside a program I need to solve this linear equation in terms of p1.
> However something odds happens. Sometimes the solution is computed and
> sometimes the result is empty [I mean no output...]. Is this a bug of the
> solve command or am I doing something wrong? The problem is robust to:
> changing name to the variables and other makeups..
>
> Thank you very much,
>
> David
That's the weirdest bug I've seen in weeks. As it happens, it's mine. At
least the inconsistent behavior, that is. I'll fix it, and maybe also
try to address the issue of how to handle approximate numbers in testing
subexpressions for zero.
I've excised your code and put in place a substantially smaller example
that I believe is responsible. The table will tend to give erratic
results.
zz = (-1.*x^7*(-1. + p - 7.*x^6 + p*x^6 + 6.*x^7)*
(7.000000000000002 - 7.000000000000002*x + 14.000000000000004*x^6 -
14.000000000000004*x^7 + 7.000000000000002*x^12 -
6.999999999999998*x^13))/
(p^2*(1. + 0.9*x^6)^2*(1. + x^6)^5);
Table [Developer`ClearCache[]; Developer`ZeroQ[zz], {10}]
One workaround would be to use exact input, say by preprocessing with
Rationalize.
Daniel Lichtblau
Wolfram Research
- References:
- Bug in Solve?
- From: "David" <none@news.acns.nwu.edu>
- Bug in Solve?