function f (s, n, createValue) {
    const c = divide(1, multiply(d(createValue(0), n), subtract(1, pow(2, subtract(1, s)))))
    let S = createValue(0)
    for (let k = createValue(1); smallerEq(k, n); k = add(k, 1)) {
      S = add(S, divide(multiply((-1) ** (k - 1), d(k, n)), pow(k, s)))
    }
    return multiply(c, S)
  }
Calculate the positive Riemann Zeta function
@param {number} s   a real or complex number with s.re > 1
@param {number} n   a positive integer
@param {(number) => number | BigNumber | Complex} createValue
@return {number}    Riemann Zeta of s