function f3 = fun3_prob1(x3,y3,alpha3,beta3)
f3 = exp(-y3/beta3)*(y3/beta3)^(floor(y3)-1)/(floor(alpha3)*gamma(beta3));
end