Provide functions to compute max(q) and/or range(cutoff)

[hofnerb]

• [ ] For assumption = "unimodal" it must hold:

cutoff >= 1/2 + min(theta^2, 1 / (2*B) + 3/4 * theta^2)

• [ ] For assumption = "r-concave": Export an interface to maxQ