std1<-sqrt(0.0036)
std2<-sqrt(0.0576)
rho<--1
#
x1<-runif(5000)
x2<-1-x1
var_p<-x1^2*std1^2 + x2^2*std2^2 +2*x1*x2*rho*std1*std2
final<-cbind(x1,x2,var_p)
final2<-final[order(final[,3]),]
head(final2)
            x1        x2        var_p
[1,] 0.7999653 0.2000347 1.081783e-10
[2,] 0.8000401 0.1999599 1.449664e-10
[3,] 0.8000787 0.1999213 5.568790e-10
[4,] 0.8001358 0.1998642 1.658932e-09
[5,] 0.8002285 0.1997715 4.698080e-09
[6,] 0.8003739 0.1996261 1.258483e-08
