function f2 = fun2_prob1(x2,y2,alpha2,beta2)
f2 = exp(-x2/beta2)*(y2/beta2)^(floor(alpha2)-1)/(floor(alpha2)*gamma(beta2));
end