def PosNormal(mean, sigma):
    x = np.random.normal(loc = mean, scale = sigma, size=1)[0]
    return(x if x>=0 else PosNormal(mean,sigma))