derivative ((1+B^2) * P * R) / (B^2 * P + R)