import numpy.random as rd
import numpy as np
import numpy.linalg as alg

def A(n) :
    return np.array([rd.binomial(1,.5,n) for _ in range(n)])
    
def estSDP(M) :
    test = np.array_equal(M,M.T)
    L = alg.eigvals(M)
    for k in L :
        test = test and (k > 0)
    return test
    
def prop(n) :
    s = 0
    for _ in range(10000) :
        if estSDP(A(n)) :
            s += 1
    return s/10000