(4-((1/2)(1-b^2)q-b+3/(2q))^2)^(1/2)(1-b^2)>=b; b=0.55701571