solve y'(x) = 1/4*(gamma^2 - b^2)*y(x)^2 + b*y(x) - 1 and y(0)=0