| T.R | Title | User | Personal
Name | Date | Lines |
|---|
| 1663.1 | Brute Force solution... | SSAG::LARY | Laughter & hope & a sock in the eye | Mon Sep 14 1992 04:44 | 32 |
| 1663.2 | an attempt at a start to solve this | STAR::ABBASI | Spell checking is a family value | Mon Sep 14 1992 05:34 | 33 |
|
all roots are real if
q^3+r^2 <= 0
where
q= 1/3 B - 1/9 A^2
and
r= 1/6 (BA-3C)- 1/27 A^3
so ask, what is probability that q^3+r^2 <=0 ?
or after expanding and simplifying:
what is probability that
1/27 B^3 - 1/108 B^2 A^2 - 1/6 BCA + 1/4 C^2 + 1/27 CA^3 <=0 ?
since B,A,C are all >= 0, then we must have
1/27 B^3 + 1/4 C^2 + 1/27 CA^3 <= 1/108 B^2 A^2 + 1/6 BCA --- (1)
for all roots to be real.
what is probability of (1) being true?
rewrite (1) as
4 B^3 + 27 C^2 + 4 CA^3 <= B^2 A^2 + 18 BCA -- (2)
so we ask, what is probability of (2) being true?
I leave the rest as an excersise for the reader, I want to go and sleep ;-)
/Nasser
|
| 1663.3 | | SSAG::LARY | Laughter & hope & a sock in the eye | Mon Sep 28 1992 17:10 | 15 |
| 1663.4 | | VMSDEV::HALLYB | Fish have no concept of fire. | Mon Sep 28 1992 19:01 | 2 |
| p.s., a hack BASIC Monte Carlo program, run for mega-trials,
numerically verified Richie's answer.
|
| 1663.5 | Onward... | SSAG::LARY | Laughter & hope & a sock in the eye | Mon Sep 28 1992 23:23 | 20 |
| Thanks for verifying it!
Now, of course, as is customary in MATH notes, we need to generalize the
problem. Let P(N) be the probability that a unary Nth degree polynomial
with uniform random coefficients in [0,1] has all real roots.
(1) Prove that P(N) > 0 when N > 0. This looks like it would be easy if
I knew enough algebra/analysis, but I can't get a handle on it...
(2) Given (1),
P(1) = 1
P(2) = 1/12 (easy to derive)
P(3) = 1/2880 (from note 1663.1)
P(4) = ??
P(N) = ?? (closed form function of N)
Sum(N=1 to inf) P(N) = ????? (numeric value)
|