Probability of forming a triangle

Estimating_TriangleProb.py

import numpy as np

Nsamples = 10000
p_estimate = 0.
s_estimate = 0.

for j in range(Nsamples):

    ps = np.random.random(2)
    x = min(ps)
    y = max(ps)

    if ( ( x <= 0.5 ) and ( y >= 0.5 ) and ( y - x <= 0.5 ) ):
        p_estimate += 1.

p_estimate /= Nsamples

print( 'P(triangle) = %5.4f' % p_estimate )
print( 'Relative error = %5.4f' % ( np.abs( p_estimate - 0.25 )/0.25 ) )