Given {a,b,c} in [-L,L]{a,b,c}[L,L] What is the probability that the roots of a x^2+b x + c = 0ax2+bx+c=0 be real?

1 Answer
Nov 19, 2016

(41+6ln(2))/72~~0.627241+6ln(2)720.6272

Explanation:

We will operate on the assumptions L>0L>0 and that the probability is equal that (a, b, c) = (x_1, x_2, x_3)(a,b,c)=(x1,x2,x3) for all (x_1, x_2, x_3) in [-L, L]^3(x1,x2,x3)[L,L]3

Initial observations:

  • ax^2+bx+c = 0ax2+bx+c=0 has real roots if and only if the discriminant of the quadratic is nonnegative, that is, b^2-4ac >= 0b24ac0, or b^2>=4acb24ac.

  • We can partition the probability space into the two cases in which "sgn"(a) = "sgn"(c)sgn(a)=sgn(c) or "sgn"(a)!="sgn"(c)sgn(a)sgn(c), and
    P("sgn"(a)="sgn"(c)) = P("sgn"(a)!="sgn"(c)) = 1/2P(sgn(a)=sgn(c))=P(sgn(a)sgn(c))=12

  • If "sgn"(a) != "sgn"(c)sgn(a)sgn(c), then 4ac <= 04ac0, meaning b^2>=4acb24ac.

  • By symmetry:
    P(b^2>=4ac|a, c<=0) = P(b^2>=4ac|a, c>=0)P(b24aca,c0)=P(b24aca,c0)

  • As b^2=(-b)^2b2=(b)2 for all bb, we can restrict bb to [0, L][0,L] without changing the probability of it falling within a certain range. Thus
    P(b^2>=4ac|a, c>=0) = P(b^2>=4ac|a, b, c>=0)P(b24aca,c0)=P(b24aca,b,c0)

With these, we can reformulate the problem as follows:

P(ax^2+bx+c" has real roots") = P(b^2>=4ac)P(ax2+bx+c has real roots)=P(b24ac)

=1/2P(b^2>=4ac|"sgn"(a)!="sgn"(c))=12P(b24acsgn(a)sgn(c))

+1/2P(b^2>=4ac|"sgn"(a)="sgn"(c))+12P(b24acsgn(a)=sgn(c))

=1/2(1)+1/2P(b^2>=4ac|a, c<0 or a, c >0)=12(1)+12P(b24aca,c<0ora,c>0)

=1/2+1/2P(b^2>=4ac|a, c>0)=12+12P(b24aca,c>0)

=1/2+1/2P(b^2>=4ac|a, b, c>0)=12+12P(b24aca,b,c>0)

=1/2+1/2P(b^2>=4ac|(a, b, c) in [0, L]^3)=12+12P(b24ac(a,b,c)[0,L]3)

We will now calculate P(b^2>=4ac|(a, b, c) in [0, L]^3)P(b24ac(a,b,c)[0,L]3)


If we consider a coordinate system in 3-space, then [0, L]^3[0,L]3 is equivalent to a cube with side length LL resting in the first octant and having sides colinear with the axes. Let SS be the solid bounded by this cube and above by the surface y^2 = 4xzy2=4xz, and V_SVS be the volume of this solid. Then the probability that b^2>=4acb24ac given an arbitrary (a, b, c) in [0, L]^3(a,b,c)[0,L]3 is equal to the probability that (a, b, c)(a,b,c) falls within SS given that (a, b, c)(a,b,c) falls within the cube, i.e. V_S/V_"cube" = V_S/L^3VSVcube=VSL3

To find V_SVS, we will first consider the area A_bAb of the slice of SS found by fixing y=by=b for some b in [0, L]b[0,L], and then we will add up all of these areas by integrating A_bAb as bb goes from 00 to LL.

If we fix y=by=b for some bb, then A_bAb is the area bounded by the lines x=0, x=L, z=0, z=Lx=0,x=L,z=0,z=L and the curve 4xz=b^24xz=b2, which we can rewrite as z=b^2/(4x)z=b24x

![desmos.com](useruploads.socratic.org)

Notice that the curve intersects the square at the points where x=Lx=L or z=Lz=L, that is, at (x, z) in{(L, b^2/(4L)), (b^2/(4L), L)}(x,z){(L,b24L),(b24L,L)}

With that, we can now set up our integrals. The upper bound from x=0x=0 to x=b^2/(4L)x=b24L is the line z=Lz=L. The upper bound from x=b^2/(4L)x=b24L to x=Lx=L is the curve z=b^2/(4x)z=b24x. Thus

A_b = int_0^(b^2/(4L))Ldx + int_(b^2/(4L))^Lb^2/(4x)dxAb=b24L0Ldx+Lb24Lb24xdx

and

V_S = int_0^LA_bdbVS=L0Abdb

Omitting the full process of integration to save space, we find the result

V_S = (5+6ln(2))/36L^3VS=5+6ln(2)36L3

Thus

P(b^2>=4ac|a, b, c in [0, L]) = V_S/L^3= (5+6ln(2))/36P(b24aca,b,c[0,L])=VSL3=5+6ln(2)36

Substituting this into our initial equation, we get our final result:

P(ax^2+bx+c" has real roots")P(ax2+bx+c has real roots)

=1/2+1/2P(b^2>=4ac|(a, b, c) in [0, L]^3)=12+12P(b24ac(a,b,c)[0,L]3)

=1/2+1/2((5+6ln(2))/36)=12+12(5+6ln(2)36)

=(41+6ln(2))/72=41+6ln(2)72

~~0.62720.6272