Exception for Octave 3.8, jcobi/2

\[\left(\alpha > -1 \land \beta > -1\right) \land i > 0\]

Math

\[\frac{\frac{\frac{\left(\alpha + \beta\right) \cdot \left(\beta - \alpha\right)}{\left(\alpha + \beta\right) + 2 \cdot i}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) + 2} + 1}{2} \]

FPCore

(FPCore (alpha beta i)
 :precision binary64
 (/
  (+
   (/
    (/ (* (+ alpha beta) (- beta alpha)) (+ (+ alpha beta) (* 2.0 i)))
    (+ (+ (+ alpha beta) (* 2.0 i)) 2.0))
   1.0)
  2.0))

C

double code(double alpha, double beta, double i) {
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
}

Fortran

real(8) function code(alpha, beta, i)
    real(8), intent (in) :: alpha
    real(8), intent (in) :: beta
    real(8), intent (in) :: i
    code = (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0d0 * i))) / (((alpha + beta) + (2.0d0 * i)) + 2.0d0)) + 1.0d0) / 2.0d0
end function

Java

public static double code(double alpha, double beta, double i) {
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
}

Python

def code(alpha, beta, i):
	return (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0

Julia

function code(alpha, beta, i)
	return Float64(Float64(Float64(Float64(Float64(Float64(alpha + beta) * Float64(beta - alpha)) / Float64(Float64(alpha + beta) + Float64(2.0 * i))) / Float64(Float64(Float64(alpha + beta) + Float64(2.0 * i)) + 2.0)) + 1.0) / 2.0)
end

MATLAB

function tmp = code(alpha, beta, i)
	tmp = (((((alpha + beta) * (beta - alpha)) / ((alpha + beta) + (2.0 * i))) / (((alpha + beta) + (2.0 * i)) + 2.0)) + 1.0) / 2.0;
end

Wolfram

code[alpha_, beta_, i_] := N[(N[(N[(N[(N[(N[(alpha + beta), $MachinePrecision] * N[(beta - alpha), $MachinePrecision]), $MachinePrecision] / N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision] + 2.0), $MachinePrecision]), $MachinePrecision] + 1.0), $MachinePrecision] / 2.0), $MachinePrecision]