Here's one way to solve the problem by using PROC OPTMODEL in SAS/OR:

proc optmodel;
var a >= 0, b >= 1 < = 3;
max prob = cdf("beta", 0.8, a, b) - cdf("beta", 0.6, a, b);
solve;
print a b prob;
quit;