23/07/20 07:49:05.42 FCPwP+wT0.net
rm(list=ls())
A=c(-2,0,0)
B=c(2,1,1)
f=\(P){
dist(rbind(P,A))+dist(rbind(P,B))
}
sim=\(){
th=runif(1,-pi,pi)
x=runif(1)*cos(th)
y=runif(1)*sin(th)
z=runif(1,0,2)
f(c(x,y,z))
}
d=replicate(1e6,sim())
max(d)
g=\(){
x=runif(1)
y=runif(1)
if(x^2+y^2>1){return(0)
}else{
z=runif(1,0,2)
return(f(c(x,y,z)))
}
}
r=replicate(1e6,g())
BEST::plotPost(r)