function fx = fun4_prob1(x,y,alpha,beta)
fx = min((fun1_prob1(y)*fun2_prob1(x,y,alpha,beta))/(fun1_prob1(x)*fun3_prob1(x,y,alpha,beta)),1);
end