{erf(sqrt(alpha/2))*sqrt(Pi*alpha/(2))*exp(alpha/2), alpha/(1-alpha)} for alpha from 0 to 2