integrate exp(-sqrt((z-y)/b))*exp(-sqrt(y/b))/sqrt(y*(z-y)) for y from 0 to z