function p = p3 ( m )
  a = randn ( 1, m );
  b = randn ( 1, m );
  c = randn ( 1, m );
  det = b.^2 - 4 * a .* c;
  i = ( 0.0 <= det );
  p = sum ( i ) / m;
  return
end