If two parameters p1=10*exp(-a^2/(2b)) and p0=2*exp(-c^2/d), the random variable X is sample from Bernoulli(p1/(p0+p1)). How to write this code in
p1=10*exp(-a^2/(2b))
p0=2*exp(-c^2/d)