Exception for Harley's example

\[0 < c_p \land 0 < c_n\]

Math

\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}} \]

FPCore

(FPCore (c_p c_n t s)
 :precision binary64
 (/
  (*
   (pow (/ 1.0 (+ 1.0 (exp (- s)))) c_p)
   (pow (- 1.0 (/ 1.0 (+ 1.0 (exp (- s))))) c_n))
  (*
   (pow (/ 1.0 (+ 1.0 (exp (- t)))) c_p)
   (pow (- 1.0 (/ 1.0 (+ 1.0 (exp (- t))))) c_n))))

C

double code(double c_p, double c_n, double t, double s) {
	return (pow((1.0 / (1.0 + exp(-s))), c_p) * pow((1.0 - (1.0 / (1.0 + exp(-s)))), c_n)) / (pow((1.0 / (1.0 + exp(-t))), c_p) * pow((1.0 - (1.0 / (1.0 + exp(-t)))), c_n));
}

Fortran

real(8) function code(c_p, c_n, t, s)
    real(8), intent (in) :: c_p
    real(8), intent (in) :: c_n
    real(8), intent (in) :: t
    real(8), intent (in) :: s
    code = (((1.0d0 / (1.0d0 + exp(-s))) ** c_p) * ((1.0d0 - (1.0d0 / (1.0d0 + exp(-s)))) ** c_n)) / (((1.0d0 / (1.0d0 + exp(-t))) ** c_p) * ((1.0d0 - (1.0d0 / (1.0d0 + exp(-t)))) ** c_n))
end function

Java

public static double code(double c_p, double c_n, double t, double s) {
	return (Math.pow((1.0 / (1.0 + Math.exp(-s))), c_p) * Math.pow((1.0 - (1.0 / (1.0 + Math.exp(-s)))), c_n)) / (Math.pow((1.0 / (1.0 + Math.exp(-t))), c_p) * Math.pow((1.0 - (1.0 / (1.0 + Math.exp(-t)))), c_n));
}

Python

def code(c_p, c_n, t, s):
	return (math.pow((1.0 / (1.0 + math.exp(-s))), c_p) * math.pow((1.0 - (1.0 / (1.0 + math.exp(-s)))), c_n)) / (math.pow((1.0 / (1.0 + math.exp(-t))), c_p) * math.pow((1.0 - (1.0 / (1.0 + math.exp(-t)))), c_n))

Julia

function code(c_p, c_n, t, s)
	return Float64(Float64((Float64(1.0 / Float64(1.0 + exp(Float64(-s)))) ^ c_p) * (Float64(1.0 - Float64(1.0 / Float64(1.0 + exp(Float64(-s))))) ^ c_n)) / Float64((Float64(1.0 / Float64(1.0 + exp(Float64(-t)))) ^ c_p) * (Float64(1.0 - Float64(1.0 / Float64(1.0 + exp(Float64(-t))))) ^ c_n)))
end

MATLAB

function tmp = code(c_p, c_n, t, s)
	tmp = (((1.0 / (1.0 + exp(-s))) ^ c_p) * ((1.0 - (1.0 / (1.0 + exp(-s)))) ^ c_n)) / (((1.0 / (1.0 + exp(-t))) ^ c_p) * ((1.0 - (1.0 / (1.0 + exp(-t)))) ^ c_n));
end

Wolfram

code[c$95$p_, c$95$n_, t_, s_] := N[(N[(N[Power[N[(1.0 / N[(1.0 + N[Exp[(-s)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], c$95$p], $MachinePrecision] * N[Power[N[(1.0 - N[(1.0 / N[(1.0 + N[Exp[(-s)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], c$95$n], $MachinePrecision]), $MachinePrecision] / N[(N[Power[N[(1.0 / N[(1.0 + N[Exp[(-t)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision], c$95$p], $MachinePrecision] * N[Power[N[(1.0 - N[(1.0 / N[(1.0 + N[Exp[(-t)], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], c$95$n], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]