Also mostly a duplicate.  ``v`` isn't used except for in the line below where it's calculated. Can be merged with the code 50 lines up, only needs an ``if b is not defined: b = (1-xbar) * fac ``.