
(FPCore (z)
:precision binary64
(let* ((t_0 (- (- 1.0 z) 1.0)) (t_1 (+ t_0 7.0)) (t_2 (+ t_1 0.5)))
(*
(/ PI (sin (* PI z)))
(*
(* (* (sqrt (* PI 2.0)) (pow t_2 (+ t_0 0.5))) (exp (- t_2)))
(+
(+
(+
(+
(+
(+
(+
(+ 0.9999999999998099 (/ 676.5203681218851 (+ t_0 1.0)))
(/ -1259.1392167224028 (+ t_0 2.0)))
(/ 771.3234287776531 (+ t_0 3.0)))
(/ -176.6150291621406 (+ t_0 4.0)))
(/ 12.507343278686905 (+ t_0 5.0)))
(/ -0.13857109526572012 (+ t_0 6.0)))
(/ 9.984369578019572e-6 t_1))
(/ 1.5056327351493116e-7 (+ t_0 8.0)))))))
double code(double z) {
double t_0 = (1.0 - z) - 1.0;
double t_1 = t_0 + 7.0;
double t_2 = t_1 + 0.5;
return (((double) M_PI) / sin((((double) M_PI) * z))) * (((sqrt((((double) M_PI) * 2.0)) * pow(t_2, (t_0 + 0.5))) * exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))));
}
public static double code(double z) {
double t_0 = (1.0 - z) - 1.0;
double t_1 = t_0 + 7.0;
double t_2 = t_1 + 0.5;
return (Math.PI / Math.sin((Math.PI * z))) * (((Math.sqrt((Math.PI * 2.0)) * Math.pow(t_2, (t_0 + 0.5))) * Math.exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))));
}
def code(z): t_0 = (1.0 - z) - 1.0 t_1 = t_0 + 7.0 t_2 = t_1 + 0.5 return (math.pi / math.sin((math.pi * z))) * (((math.sqrt((math.pi * 2.0)) * math.pow(t_2, (t_0 + 0.5))) * math.exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))))
function code(z) t_0 = Float64(Float64(1.0 - z) - 1.0) t_1 = Float64(t_0 + 7.0) t_2 = Float64(t_1 + 0.5) return Float64(Float64(pi / sin(Float64(pi * z))) * Float64(Float64(Float64(sqrt(Float64(pi * 2.0)) * (t_2 ^ Float64(t_0 + 0.5))) * exp(Float64(-t_2))) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(t_0 + 1.0))) + Float64(-1259.1392167224028 / Float64(t_0 + 2.0))) + Float64(771.3234287776531 / Float64(t_0 + 3.0))) + Float64(-176.6150291621406 / Float64(t_0 + 4.0))) + Float64(12.507343278686905 / Float64(t_0 + 5.0))) + Float64(-0.13857109526572012 / Float64(t_0 + 6.0))) + Float64(9.984369578019572e-6 / t_1)) + Float64(1.5056327351493116e-7 / Float64(t_0 + 8.0))))) end
function tmp = code(z) t_0 = (1.0 - z) - 1.0; t_1 = t_0 + 7.0; t_2 = t_1 + 0.5; tmp = (pi / sin((pi * z))) * (((sqrt((pi * 2.0)) * (t_2 ^ (t_0 + 0.5))) * exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0)))); end
code[z_] := Block[{t$95$0 = N[(N[(1.0 - z), $MachinePrecision] - 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + 7.0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 + 0.5), $MachinePrecision]}, N[(N[(Pi / N[Sin[N[(Pi * z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Sqrt[N[(Pi * 2.0), $MachinePrecision]], $MachinePrecision] * N[Power[t$95$2, N[(t$95$0 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Exp[(-t$95$2)], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(771.3234287776531 / N[(t$95$0 + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-176.6150291621406 / N[(t$95$0 + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(12.507343278686905 / N[(t$95$0 + 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-0.13857109526572012 / N[(t$95$0 + 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(t$95$0 + 8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - z\right) - 1\\
t_1 := t\_0 + 7\\
t_2 := t\_1 + 0.5\\
\frac{\pi}{\sin \left(\pi \cdot z\right)} \cdot \left(\left(\left(\sqrt{\pi \cdot 2} \cdot {t\_2}^{\left(t\_0 + 0.5\right)}\right) \cdot e^{-t\_2}\right) \cdot \left(\left(\left(\left(\left(\left(\left(\left(0.9999999999998099 + \frac{676.5203681218851}{t\_0 + 1}\right) + \frac{-1259.1392167224028}{t\_0 + 2}\right) + \frac{771.3234287776531}{t\_0 + 3}\right) + \frac{-176.6150291621406}{t\_0 + 4}\right) + \frac{12.507343278686905}{t\_0 + 5}\right) + \frac{-0.13857109526572012}{t\_0 + 6}\right) + \frac{9.984369578019572 \cdot 10^{-6}}{t\_1}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{t\_0 + 8}\right)\right)
\end{array}
\end{array}
Sampling outcomes in binary64 precision:
Herbie found 20 alternatives:
| Alternative | Accuracy | Speedup |
|---|
(FPCore (z)
:precision binary64
(let* ((t_0 (- (- 1.0 z) 1.0)) (t_1 (+ t_0 7.0)) (t_2 (+ t_1 0.5)))
(*
(/ PI (sin (* PI z)))
(*
(* (* (sqrt (* PI 2.0)) (pow t_2 (+ t_0 0.5))) (exp (- t_2)))
(+
(+
(+
(+
(+
(+
(+
(+ 0.9999999999998099 (/ 676.5203681218851 (+ t_0 1.0)))
(/ -1259.1392167224028 (+ t_0 2.0)))
(/ 771.3234287776531 (+ t_0 3.0)))
(/ -176.6150291621406 (+ t_0 4.0)))
(/ 12.507343278686905 (+ t_0 5.0)))
(/ -0.13857109526572012 (+ t_0 6.0)))
(/ 9.984369578019572e-6 t_1))
(/ 1.5056327351493116e-7 (+ t_0 8.0)))))))
double code(double z) {
double t_0 = (1.0 - z) - 1.0;
double t_1 = t_0 + 7.0;
double t_2 = t_1 + 0.5;
return (((double) M_PI) / sin((((double) M_PI) * z))) * (((sqrt((((double) M_PI) * 2.0)) * pow(t_2, (t_0 + 0.5))) * exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))));
}
public static double code(double z) {
double t_0 = (1.0 - z) - 1.0;
double t_1 = t_0 + 7.0;
double t_2 = t_1 + 0.5;
return (Math.PI / Math.sin((Math.PI * z))) * (((Math.sqrt((Math.PI * 2.0)) * Math.pow(t_2, (t_0 + 0.5))) * Math.exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))));
}
def code(z): t_0 = (1.0 - z) - 1.0 t_1 = t_0 + 7.0 t_2 = t_1 + 0.5 return (math.pi / math.sin((math.pi * z))) * (((math.sqrt((math.pi * 2.0)) * math.pow(t_2, (t_0 + 0.5))) * math.exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0))))
function code(z) t_0 = Float64(Float64(1.0 - z) - 1.0) t_1 = Float64(t_0 + 7.0) t_2 = Float64(t_1 + 0.5) return Float64(Float64(pi / sin(Float64(pi * z))) * Float64(Float64(Float64(sqrt(Float64(pi * 2.0)) * (t_2 ^ Float64(t_0 + 0.5))) * exp(Float64(-t_2))) * Float64(Float64(Float64(Float64(Float64(Float64(Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(t_0 + 1.0))) + Float64(-1259.1392167224028 / Float64(t_0 + 2.0))) + Float64(771.3234287776531 / Float64(t_0 + 3.0))) + Float64(-176.6150291621406 / Float64(t_0 + 4.0))) + Float64(12.507343278686905 / Float64(t_0 + 5.0))) + Float64(-0.13857109526572012 / Float64(t_0 + 6.0))) + Float64(9.984369578019572e-6 / t_1)) + Float64(1.5056327351493116e-7 / Float64(t_0 + 8.0))))) end
function tmp = code(z) t_0 = (1.0 - z) - 1.0; t_1 = t_0 + 7.0; t_2 = t_1 + 0.5; tmp = (pi / sin((pi * z))) * (((sqrt((pi * 2.0)) * (t_2 ^ (t_0 + 0.5))) * exp(-t_2)) * ((((((((0.9999999999998099 + (676.5203681218851 / (t_0 + 1.0))) + (-1259.1392167224028 / (t_0 + 2.0))) + (771.3234287776531 / (t_0 + 3.0))) + (-176.6150291621406 / (t_0 + 4.0))) + (12.507343278686905 / (t_0 + 5.0))) + (-0.13857109526572012 / (t_0 + 6.0))) + (9.984369578019572e-6 / t_1)) + (1.5056327351493116e-7 / (t_0 + 8.0)))); end
code[z_] := Block[{t$95$0 = N[(N[(1.0 - z), $MachinePrecision] - 1.0), $MachinePrecision]}, Block[{t$95$1 = N[(t$95$0 + 7.0), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 + 0.5), $MachinePrecision]}, N[(N[(Pi / N[Sin[N[(Pi * z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Sqrt[N[(Pi * 2.0), $MachinePrecision]], $MachinePrecision] * N[Power[t$95$2, N[(t$95$0 + 0.5), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[Exp[(-t$95$2)], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(N[(N[(N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(t$95$0 + 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(771.3234287776531 / N[(t$95$0 + 3.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-176.6150291621406 / N[(t$95$0 + 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(12.507343278686905 / N[(t$95$0 + 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(-0.13857109526572012 / N[(t$95$0 + 6.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / t$95$1), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(t$95$0 + 8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\begin{array}{l}
\\
\begin{array}{l}
t_0 := \left(1 - z\right) - 1\\
t_1 := t\_0 + 7\\
t_2 := t\_1 + 0.5\\
\frac{\pi}{\sin \left(\pi \cdot z\right)} \cdot \left(\left(\left(\sqrt{\pi \cdot 2} \cdot {t\_2}^{\left(t\_0 + 0.5\right)}\right) \cdot e^{-t\_2}\right) \cdot \left(\left(\left(\left(\left(\left(\left(\left(0.9999999999998099 + \frac{676.5203681218851}{t\_0 + 1}\right) + \frac{-1259.1392167224028}{t\_0 + 2}\right) + \frac{771.3234287776531}{t\_0 + 3}\right) + \frac{-176.6150291621406}{t\_0 + 4}\right) + \frac{12.507343278686905}{t\_0 + 5}\right) + \frac{-0.13857109526572012}{t\_0 + 6}\right) + \frac{9.984369578019572 \cdot 10^{-6}}{t\_1}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{t\_0 + 8}\right)\right)
\end{array}
\end{array}
(FPCore (z)
:precision binary64
(*
(/
(*
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(+
(-
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+ (/ -771.3234287776531 (- 3.0 z)) (/ -1259.1392167224028 (- z 2.0))))
(+
(+ (/ -176.6150291621406 (- 4.0 z)) (/ -12.507343278686905 (- z 5.0)))
(/ 1.5056327351493116e-7 (- 8.0 z)))))
PI)
(sin (* z PI)))
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (exp (- 7.5 z)))))
double code(double z) {
return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * ((double) M_PI)) / sin((z * ((double) M_PI)))) * ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / exp((7.5 - z)));
}
public static double code(double z) {
return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * Math.PI) / Math.sin((z * Math.PI))) * ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.exp((7.5 - z)));
}
def code(z): return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * math.pi) / math.sin((z * math.pi))) * ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.exp((7.5 - z)))
function code(z) return Float64(Float64(Float64(Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) - Float64(Float64(-771.3234287776531 / Float64(3.0 - z)) + Float64(-1259.1392167224028 / Float64(z - 2.0)))) + Float64(Float64(Float64(-176.6150291621406 / Float64(4.0 - z)) + Float64(-12.507343278686905 / Float64(z - 5.0))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))))) * pi) / sin(Float64(z * pi))) * Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / exp(Float64(7.5 - z)))) end
function tmp = code(z) tmp = (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * pi) / sin((z * pi))) * ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / exp((7.5 - z))); end
code[z_] := N[(N[(N[(N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(-771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(z - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-176.6150291621406 / N[(4.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-12.507343278686905 / N[(z - 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * Pi), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \left(\left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) - \left(\frac{-771.3234287776531}{3 - z} + \frac{-1259.1392167224028}{z - 2}\right)\right) + \left(\left(\frac{-176.6150291621406}{4 - z} + \frac{-12.507343278686905}{z - 5}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right)\right)\right) \cdot \pi}{\sin \left(z \cdot \pi\right)} \cdot \frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{e^{7.5 - z}}
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Final simplification99.1%
(FPCore (z)
:precision binary64
(*
(/
(*
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(+
(-
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+ (/ -771.3234287776531 (- 3.0 z)) (/ -1259.1392167224028 (- z 2.0))))
(+
(+ (/ -176.6150291621406 (- 4.0 z)) (/ -12.507343278686905 (- z 5.0)))
(/ 1.5056327351493116e-7 (- 8.0 z)))))
PI)
(exp (- 7.5 z)))
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))))
double code(double z) {
return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * ((double) M_PI)) / exp((7.5 - z))) * ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI))));
}
public static double code(double z) {
return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * Math.PI) / Math.exp((7.5 - z))) * ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI)));
}
def code(z): return (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * math.pi) / math.exp((7.5 - z))) * ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi)))
function code(z) return Float64(Float64(Float64(Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) - Float64(Float64(-771.3234287776531 / Float64(3.0 - z)) + Float64(-1259.1392167224028 / Float64(z - 2.0)))) + Float64(Float64(Float64(-176.6150291621406 / Float64(4.0 - z)) + Float64(-12.507343278686905 / Float64(z - 5.0))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))))) * pi) / exp(Float64(7.5 - z))) * Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi)))) end
function tmp = code(z) tmp = (((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + (1.5056327351493116e-7 / (8.0 - z))))) * pi) / exp((7.5 - z))) * ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))); end
code[z_] := N[(N[(N[(N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(-771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(z - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-176.6150291621406 / N[(4.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-12.507343278686905 / N[(z - 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * Pi), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \left(\left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) - \left(\frac{-771.3234287776531}{3 - z} + \frac{-1259.1392167224028}{z - 2}\right)\right) + \left(\left(\frac{-176.6150291621406}{4 - z} + \frac{-12.507343278686905}{z - 5}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right)\right)\right) \cdot \pi}{e^{7.5 - z}} \cdot \frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)}
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.0%
Final simplification99.0%
(FPCore (z)
:precision binary64
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(*
(+
(/ -0.13857109526572012 (- 6.0 z))
(+
(-
(/ 9.984369578019572e-6 (- 7.0 z))
(+
(+ (/ 12.507343278686905 (- z 5.0)) (/ -176.6150291621406 (- z 4.0)))
(/ 1.5056327351493116e-7 (- z 8.0))))
(+
0.9999999999998099
(+
(/ 676.5203681218851 (- 1.0 z))
(+
(/ -1259.1392167224028 (- 2.0 z))
(/ 771.3234287776531 (- 3.0 z)))))))
(/ PI (exp (- 7.5 z))))))
double code(double z) {
return ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * (((-0.13857109526572012 / (6.0 - z)) + (((9.984369578019572e-6 / (7.0 - z)) - (((12.507343278686905 / (z - 5.0)) + (-176.6150291621406 / (z - 4.0))) + (1.5056327351493116e-7 / (z - 8.0)))) + (0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))))) * (((double) M_PI) / exp((7.5 - z))));
}
public static double code(double z) {
return ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * (((-0.13857109526572012 / (6.0 - z)) + (((9.984369578019572e-6 / (7.0 - z)) - (((12.507343278686905 / (z - 5.0)) + (-176.6150291621406 / (z - 4.0))) + (1.5056327351493116e-7 / (z - 8.0)))) + (0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))))) * (Math.PI / Math.exp((7.5 - z))));
}
def code(z): return ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * (((-0.13857109526572012 / (6.0 - z)) + (((9.984369578019572e-6 / (7.0 - z)) - (((12.507343278686905 / (z - 5.0)) + (-176.6150291621406 / (z - 4.0))) + (1.5056327351493116e-7 / (z - 8.0)))) + (0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))))) * (math.pi / math.exp((7.5 - z))))
function code(z) return Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(Float64(Float64(9.984369578019572e-6 / Float64(7.0 - z)) - Float64(Float64(Float64(12.507343278686905 / Float64(z - 5.0)) + Float64(-176.6150291621406 / Float64(z - 4.0))) + Float64(1.5056327351493116e-7 / Float64(z - 8.0)))) + Float64(0.9999999999998099 + Float64(Float64(676.5203681218851 / Float64(1.0 - z)) + Float64(Float64(-1259.1392167224028 / Float64(2.0 - z)) + Float64(771.3234287776531 / Float64(3.0 - z))))))) * Float64(pi / exp(Float64(7.5 - z))))) end
function tmp = code(z) tmp = ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * (((-0.13857109526572012 / (6.0 - z)) + (((9.984369578019572e-6 / (7.0 - z)) - (((12.507343278686905 / (z - 5.0)) + (-176.6150291621406 / (z - 4.0))) + (1.5056327351493116e-7 / (z - 8.0)))) + (0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))))) * (pi / exp((7.5 - z)))); end
code[z_] := N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(12.507343278686905 / N[(z - 5.0), $MachinePrecision]), $MachinePrecision] + N[(-176.6150291621406 / N[(z - 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(z - 8.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(0.9999999999998099 + N[(N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(-1259.1392167224028 / N[(2.0 - z), $MachinePrecision]), $MachinePrecision] + N[(771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(Pi / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \left(\left(\frac{-0.13857109526572012}{6 - z} + \left(\left(\frac{9.984369578019572 \cdot 10^{-6}}{7 - z} - \left(\left(\frac{12.507343278686905}{z - 5} + \frac{-176.6150291621406}{z - 4}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{z - 8}\right)\right) + \left(0.9999999999998099 + \left(\frac{676.5203681218851}{1 - z} + \left(\frac{-1259.1392167224028}{2 - z} + \frac{771.3234287776531}{3 - z}\right)\right)\right)\right)\right) \cdot \frac{\pi}{e^{7.5 - z}}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Applied egg-rr99.0%
Final simplification99.0%
(FPCore (z)
:precision binary64
(*
(+
(/ 12.507343278686905 (- 5.0 z))
(+
(+
(/ -1259.1392167224028 (- 2.0 z))
(+
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+ (/ 771.3234287776531 (- 3.0 z)) (/ 176.6150291621406 (- z 4.0)))))
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(/ 1.5056327351493116e-7 (- 8.0 z)))))
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(/ PI (exp (- 7.5 z))))))
double code(double z) {
return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * (((double) M_PI) / exp((7.5 - z))));
}
public static double code(double z) {
return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * (Math.PI / Math.exp((7.5 - z))));
}
def code(z): return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * (math.pi / math.exp((7.5 - z))))
function code(z) return Float64(Float64(Float64(12.507343278686905 / Float64(5.0 - z)) + Float64(Float64(Float64(-1259.1392167224028 / Float64(2.0 - z)) + Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) + Float64(Float64(771.3234287776531 / Float64(3.0 - z)) + Float64(176.6150291621406 / Float64(z - 4.0))))) + Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))))) * Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(pi / exp(Float64(7.5 - z))))) end
function tmp = code(z) tmp = ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * (pi / exp((7.5 - z)))); end
code[z_] := N[(N[(N[(12.507343278686905 / N[(5.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-1259.1392167224028 / N[(2.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(176.6150291621406 / N[(z - 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(Pi / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{12.507343278686905}{5 - z} + \left(\left(\frac{-1259.1392167224028}{2 - z} + \left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) + \left(\frac{771.3234287776531}{3 - z} + \frac{176.6150291621406}{z - 4}\right)\right)\right) + \left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right)\right)\right) \cdot \left(\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \frac{\pi}{e^{7.5 - z}}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr98.9%
*-commutativeN/A
associate-/l/N/A
*-commutativeN/A
associate-*l/N/A
times-fracN/A
*-lowering-*.f64N/A
Applied egg-rr98.9%
Final simplification98.9%
(FPCore (z)
:precision binary64
(*
(+
(/ 12.507343278686905 (- 5.0 z))
(+
(+
(/ -1259.1392167224028 (- 2.0 z))
(+
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+ (/ 771.3234287776531 (- 3.0 z)) (/ 176.6150291621406 (- z 4.0)))))
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(/ 1.5056327351493116e-7 (- 8.0 z)))))
(*
PI
(/
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (exp (- 7.5 z)))
(sin (* z PI))))))
double code(double z) {
return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (((double) M_PI) * (((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / exp((7.5 - z))) / sin((z * ((double) M_PI)))));
}
public static double code(double z) {
return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (Math.PI * (((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.exp((7.5 - z))) / Math.sin((z * Math.PI))));
}
def code(z): return ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (math.pi * (((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.exp((7.5 - z))) / math.sin((z * math.pi))))
function code(z) return Float64(Float64(Float64(12.507343278686905 / Float64(5.0 - z)) + Float64(Float64(Float64(-1259.1392167224028 / Float64(2.0 - z)) + Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) + Float64(Float64(771.3234287776531 / Float64(3.0 - z)) + Float64(176.6150291621406 / Float64(z - 4.0))))) + Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))))) * Float64(pi * Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / exp(Float64(7.5 - z))) / sin(Float64(z * pi))))) end
function tmp = code(z) tmp = ((12.507343278686905 / (5.0 - z)) + (((-1259.1392167224028 / (2.0 - z)) + ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) + ((771.3234287776531 / (3.0 - z)) + (176.6150291621406 / (z - 4.0))))) + (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))))) * (pi * (((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / exp((7.5 - z))) / sin((z * pi)))); end
code[z_] := N[(N[(N[(12.507343278686905 / N[(5.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-1259.1392167224028 / N[(2.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(176.6150291621406 / N[(z - 4.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(Pi * N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\frac{12.507343278686905}{5 - z} + \left(\left(\frac{-1259.1392167224028}{2 - z} + \left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) + \left(\frac{771.3234287776531}{3 - z} + \frac{176.6150291621406}{z - 4}\right)\right)\right) + \left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right)\right)\right) \cdot \left(\pi \cdot \frac{\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{e^{7.5 - z}}}{\sin \left(z \cdot \pi\right)}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr98.9%
Final simplification98.9%
(FPCore (z)
:precision binary64
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(/
(*
PI
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(+
(-
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+ (/ -771.3234287776531 (- 3.0 z)) (/ -1259.1392167224028 (- z 2.0))))
(+
(+ (/ -176.6150291621406 (- 4.0 z)) (/ -12.507343278686905 (- z 5.0)))
1.8820409189366395e-8))))
(exp (- 7.5 z)))))
double code(double z) {
return ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * ((((double) M_PI) * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + 1.8820409189366395e-8)))) / exp((7.5 - z)));
}
public static double code(double z) {
return ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * ((Math.PI * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + 1.8820409189366395e-8)))) / Math.exp((7.5 - z)));
}
def code(z): return ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * ((math.pi * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + 1.8820409189366395e-8)))) / math.exp((7.5 - z)))
function code(z) return Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(Float64(pi * Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) - Float64(Float64(-771.3234287776531 / Float64(3.0 - z)) + Float64(-1259.1392167224028 / Float64(z - 2.0)))) + Float64(Float64(Float64(-176.6150291621406 / Float64(4.0 - z)) + Float64(-12.507343278686905 / Float64(z - 5.0))) + 1.8820409189366395e-8)))) / exp(Float64(7.5 - z)))) end
function tmp = code(z) tmp = ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * ((pi * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0)))) + (((-176.6150291621406 / (4.0 - z)) + (-12.507343278686905 / (z - 5.0))) + 1.8820409189366395e-8)))) / exp((7.5 - z))); end
code[z_] := N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(Pi * N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(-771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(z - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-176.6150291621406 / N[(4.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-12.507343278686905 / N[(z - 5.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.8820409189366395e-8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \frac{\pi \cdot \left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \left(\left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) - \left(\frac{-771.3234287776531}{3 - z} + \frac{-1259.1392167224028}{z - 2}\right)\right) + \left(\left(\frac{-176.6150291621406}{4 - z} + \frac{-12.507343278686905}{z - 5}\right) + 1.8820409189366395 \cdot 10^{-8}\right)\right)\right)}{e^{7.5 - z}}
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.0%
Taylor expanded in z around 0
Simplified98.8%
Final simplification98.8%
(FPCore (z)
:precision binary64
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(*
(/ PI (exp (- 7.5 z)))
(+
(/ -0.13857109526572012 (- 6.0 z))
(+
(+
0.9999999999998099
(+
(/ 676.5203681218851 (- 1.0 z))
(+ (/ -1259.1392167224028 (- 2.0 z)) (/ 771.3234287776531 (- 3.0 z)))))
(+
(/ 9.984369578019572e-6 (- 7.0 z))
(+
(+ (/ -176.6150291621406 (- 4.0 z)) (/ 12.507343278686905 (- 5.0 z)))
1.8820409189366395e-8)))))))
double code(double z) {
return ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * ((((double) M_PI) / exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + (((-176.6150291621406 / (4.0 - z)) + (12.507343278686905 / (5.0 - z))) + 1.8820409189366395e-8)))));
}
public static double code(double z) {
return ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * ((Math.PI / Math.exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + (((-176.6150291621406 / (4.0 - z)) + (12.507343278686905 / (5.0 - z))) + 1.8820409189366395e-8)))));
}
def code(z): return ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * ((math.pi / math.exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + (((-176.6150291621406 / (4.0 - z)) + (12.507343278686905 / (5.0 - z))) + 1.8820409189366395e-8)))))
function code(z) return Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(Float64(pi / exp(Float64(7.5 - z))) * Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(Float64(0.9999999999998099 + Float64(Float64(676.5203681218851 / Float64(1.0 - z)) + Float64(Float64(-1259.1392167224028 / Float64(2.0 - z)) + Float64(771.3234287776531 / Float64(3.0 - z))))) + Float64(Float64(9.984369578019572e-6 / Float64(7.0 - z)) + Float64(Float64(Float64(-176.6150291621406 / Float64(4.0 - z)) + Float64(12.507343278686905 / Float64(5.0 - z))) + 1.8820409189366395e-8)))))) end
function tmp = code(z) tmp = ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * ((pi / exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + (((-176.6150291621406 / (4.0 - z)) + (12.507343278686905 / (5.0 - z))) + 1.8820409189366395e-8))))); end
code[z_] := N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(Pi / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(0.9999999999998099 + N[(N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(-1259.1392167224028 / N[(2.0 - z), $MachinePrecision]), $MachinePrecision] + N[(771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(-176.6150291621406 / N[(4.0 - z), $MachinePrecision]), $MachinePrecision] + N[(12.507343278686905 / N[(5.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + 1.8820409189366395e-8), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \left(\frac{\pi}{e^{7.5 - z}} \cdot \left(\frac{-0.13857109526572012}{6 - z} + \left(\left(0.9999999999998099 + \left(\frac{676.5203681218851}{1 - z} + \left(\frac{-1259.1392167224028}{2 - z} + \frac{771.3234287776531}{3 - z}\right)\right)\right) + \left(\frac{9.984369578019572 \cdot 10^{-6}}{7 - z} + \left(\left(\frac{-176.6150291621406}{4 - z} + \frac{12.507343278686905}{5 - z}\right) + 1.8820409189366395 \cdot 10^{-8}\right)\right)\right)\right)\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Applied egg-rr99.0%
Taylor expanded in z around 0
Simplified98.8%
Final simplification98.8%
(FPCore (z)
:precision binary64
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(/
(*
PI
(-
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(-
(-
(/ 1.5056327351493116e-7 (- z 8.0))
(+ (* z -10.53814559148631) -41.65228863479777))
(-
(+ 0.9999999999998099 (/ 676.5203681218851 (- 1.0 z)))
(+
(/ -771.3234287776531 (- 3.0 z))
(/ -1259.1392167224028 (- z 2.0)))))))
(exp (- 7.5 z)))))
double code(double z) {
return ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * ((((double) M_PI) * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) - (((1.5056327351493116e-7 / (z - 8.0)) - ((z * -10.53814559148631) + -41.65228863479777)) - ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0))))))) / exp((7.5 - z)));
}
public static double code(double z) {
return ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * ((Math.PI * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) - (((1.5056327351493116e-7 / (z - 8.0)) - ((z * -10.53814559148631) + -41.65228863479777)) - ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0))))))) / Math.exp((7.5 - z)));
}
def code(z): return ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * ((math.pi * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) - (((1.5056327351493116e-7 / (z - 8.0)) - ((z * -10.53814559148631) + -41.65228863479777)) - ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0))))))) / math.exp((7.5 - z)))
function code(z) return Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(Float64(pi * Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) - Float64(Float64(Float64(1.5056327351493116e-7 / Float64(z - 8.0)) - Float64(Float64(z * -10.53814559148631) + -41.65228863479777)) - Float64(Float64(0.9999999999998099 + Float64(676.5203681218851 / Float64(1.0 - z))) - Float64(Float64(-771.3234287776531 / Float64(3.0 - z)) + Float64(-1259.1392167224028 / Float64(z - 2.0))))))) / exp(Float64(7.5 - z)))) end
function tmp = code(z) tmp = ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * ((pi * (((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) - (((1.5056327351493116e-7 / (z - 8.0)) - ((z * -10.53814559148631) + -41.65228863479777)) - ((0.9999999999998099 + (676.5203681218851 / (1.0 - z))) - ((-771.3234287776531 / (3.0 - z)) + (-1259.1392167224028 / (z - 2.0))))))) / exp((7.5 - z))); end
code[z_] := N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(Pi * N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(1.5056327351493116e-7 / N[(z - 8.0), $MachinePrecision]), $MachinePrecision] - N[(N[(z * -10.53814559148631), $MachinePrecision] + -41.65228863479777), $MachinePrecision]), $MachinePrecision] - N[(N[(0.9999999999998099 + N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - N[(N[(-771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision] + N[(-1259.1392167224028 / N[(z - 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \frac{\pi \cdot \left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) - \left(\left(\frac{1.5056327351493116 \cdot 10^{-7}}{z - 8} - \left(z \cdot -10.53814559148631 + -41.65228863479777\right)\right) - \left(\left(0.9999999999998099 + \frac{676.5203681218851}{1 - z}\right) - \left(\frac{-771.3234287776531}{3 - z} + \frac{-1259.1392167224028}{z - 2}\right)\right)\right)\right)}{e^{7.5 - z}}
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.0%
Taylor expanded in z around 0
sub-negN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval98.1%
Simplified98.1%
Final simplification98.1%
(FPCore (z)
:precision binary64
(*
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (sin (* z PI)))
(*
(/ PI (exp (- 7.5 z)))
(+
(/ -0.13857109526572012 (- 6.0 z))
(+
(+
0.9999999999998099
(+
(/ 676.5203681218851 (- 1.0 z))
(+ (/ -1259.1392167224028 (- 2.0 z)) (/ 771.3234287776531 (- 3.0 z)))))
(+
(/ 9.984369578019572e-6 (- 7.0 z))
(+
(/ 1.5056327351493116e-7 (- 8.0 z))
(+ (* z -10.53814559148631) -41.65228863479777))))))))
double code(double z) {
return ((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / sin((z * ((double) M_PI)))) * ((((double) M_PI) / exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + ((1.5056327351493116e-7 / (8.0 - z)) + ((z * -10.53814559148631) + -41.65228863479777))))));
}
public static double code(double z) {
return ((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.sin((z * Math.PI))) * ((Math.PI / Math.exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + ((1.5056327351493116e-7 / (8.0 - z)) + ((z * -10.53814559148631) + -41.65228863479777))))));
}
def code(z): return ((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.sin((z * math.pi))) * ((math.pi / math.exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + ((1.5056327351493116e-7 / (8.0 - z)) + ((z * -10.53814559148631) + -41.65228863479777))))))
function code(z) return Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / sin(Float64(z * pi))) * Float64(Float64(pi / exp(Float64(7.5 - z))) * Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(Float64(0.9999999999998099 + Float64(Float64(676.5203681218851 / Float64(1.0 - z)) + Float64(Float64(-1259.1392167224028 / Float64(2.0 - z)) + Float64(771.3234287776531 / Float64(3.0 - z))))) + Float64(Float64(9.984369578019572e-6 / Float64(7.0 - z)) + Float64(Float64(1.5056327351493116e-7 / Float64(8.0 - z)) + Float64(Float64(z * -10.53814559148631) + -41.65228863479777))))))) end
function tmp = code(z) tmp = ((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / sin((z * pi))) * ((pi / exp((7.5 - z))) * ((-0.13857109526572012 / (6.0 - z)) + ((0.9999999999998099 + ((676.5203681218851 / (1.0 - z)) + ((-1259.1392167224028 / (2.0 - z)) + (771.3234287776531 / (3.0 - z))))) + ((9.984369578019572e-6 / (7.0 - z)) + ((1.5056327351493116e-7 / (8.0 - z)) + ((z * -10.53814559148631) + -41.65228863479777)))))); end
code[z_] := N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(Pi / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(0.9999999999998099 + N[(N[(676.5203681218851 / N[(1.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(-1259.1392167224028 / N[(2.0 - z), $MachinePrecision]), $MachinePrecision] + N[(771.3234287776531 / N[(3.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision] + N[(N[(z * -10.53814559148631), $MachinePrecision] + -41.65228863479777), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{\sin \left(z \cdot \pi\right)} \cdot \left(\frac{\pi}{e^{7.5 - z}} \cdot \left(\frac{-0.13857109526572012}{6 - z} + \left(\left(0.9999999999998099 + \left(\frac{676.5203681218851}{1 - z} + \left(\frac{-1259.1392167224028}{2 - z} + \frac{771.3234287776531}{3 - z}\right)\right)\right) + \left(\frac{9.984369578019572 \cdot 10^{-6}}{7 - z} + \left(\frac{1.5056327351493116 \cdot 10^{-7}}{8 - z} + \left(z \cdot -10.53814559148631 + -41.65228863479777\right)\right)\right)\right)\right)\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Applied egg-rr99.0%
Taylor expanded in z around 0
sub-negN/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
metadata-eval98.1%
Simplified98.1%
Final simplification98.1%
(FPCore (z)
:precision binary64
(*
(*
PI
(/
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (exp (- 7.5 z)))
(sin (* z PI))))
(+
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(/ 1.5056327351493116e-7 (- 8.0 z)))
(+
(/ 12.507343278686905 (- 5.0 z))
(+
260.9048120626994
(* z (+ 436.3997278161676 (* z 544.9358906000987))))))))
double code(double z) {
return (((double) M_PI) * (((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / exp((7.5 - z))) / sin((z * ((double) M_PI))))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + ((12.507343278686905 / (5.0 - z)) + (260.9048120626994 + (z * (436.3997278161676 + (z * 544.9358906000987))))));
}
public static double code(double z) {
return (Math.PI * (((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.exp((7.5 - z))) / Math.sin((z * Math.PI)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + ((12.507343278686905 / (5.0 - z)) + (260.9048120626994 + (z * (436.3997278161676 + (z * 544.9358906000987))))));
}
def code(z): return (math.pi * (((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.exp((7.5 - z))) / math.sin((z * math.pi)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + ((12.507343278686905 / (5.0 - z)) + (260.9048120626994 + (z * (436.3997278161676 + (z * 544.9358906000987))))))
function code(z) return Float64(Float64(pi * Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / exp(Float64(7.5 - z))) / sin(Float64(z * pi)))) * Float64(Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))) + Float64(Float64(12.507343278686905 / Float64(5.0 - z)) + Float64(260.9048120626994 + Float64(z * Float64(436.3997278161676 + Float64(z * 544.9358906000987))))))) end
function tmp = code(z) tmp = (pi * (((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / exp((7.5 - z))) / sin((z * pi)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + ((12.507343278686905 / (5.0 - z)) + (260.9048120626994 + (z * (436.3997278161676 + (z * 544.9358906000987)))))); end
code[z_] := N[(N[(Pi * N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(12.507343278686905 / N[(5.0 - z), $MachinePrecision]), $MachinePrecision] + N[(260.9048120626994 + N[(z * N[(436.3997278161676 + N[(z * 544.9358906000987), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\pi \cdot \frac{\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{e^{7.5 - z}}}{\sin \left(z \cdot \pi\right)}\right) \cdot \left(\left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right) + \left(\frac{12.507343278686905}{5 - z} + \left(260.9048120626994 + z \cdot \left(436.3997278161676 + z \cdot 544.9358906000987\right)\right)\right)\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6497.9%
Simplified97.9%
Final simplification97.9%
(FPCore (z)
:precision binary64
(*
(*
PI
(/
(/ (* (sqrt (* 2.0 PI)) (pow (- 7.5 z) (- 0.5 z))) (exp (- 7.5 z)))
(sin (* z PI))))
(+
(+
(+ (/ -0.13857109526572012 (- 6.0 z)) (/ 9.984369578019572e-6 (- 7.0 z)))
(/ 1.5056327351493116e-7 (- 8.0 z)))
(+
263.4062807184368
(*
z
(+
436.9000215473151
(* z (+ 545.0359493463282 (* z 606.6767878347069)))))))))
double code(double z) {
return (((double) M_PI) * (((sqrt((2.0 * ((double) M_PI))) * pow((7.5 - z), (0.5 - z))) / exp((7.5 - z))) / sin((z * ((double) M_PI))))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + (263.4062807184368 + (z * (436.9000215473151 + (z * (545.0359493463282 + (z * 606.6767878347069)))))));
}
public static double code(double z) {
return (Math.PI * (((Math.sqrt((2.0 * Math.PI)) * Math.pow((7.5 - z), (0.5 - z))) / Math.exp((7.5 - z))) / Math.sin((z * Math.PI)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + (263.4062807184368 + (z * (436.9000215473151 + (z * (545.0359493463282 + (z * 606.6767878347069)))))));
}
def code(z): return (math.pi * (((math.sqrt((2.0 * math.pi)) * math.pow((7.5 - z), (0.5 - z))) / math.exp((7.5 - z))) / math.sin((z * math.pi)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + (263.4062807184368 + (z * (436.9000215473151 + (z * (545.0359493463282 + (z * 606.6767878347069)))))))
function code(z) return Float64(Float64(pi * Float64(Float64(Float64(sqrt(Float64(2.0 * pi)) * (Float64(7.5 - z) ^ Float64(0.5 - z))) / exp(Float64(7.5 - z))) / sin(Float64(z * pi)))) * Float64(Float64(Float64(Float64(-0.13857109526572012 / Float64(6.0 - z)) + Float64(9.984369578019572e-6 / Float64(7.0 - z))) + Float64(1.5056327351493116e-7 / Float64(8.0 - z))) + Float64(263.4062807184368 + Float64(z * Float64(436.9000215473151 + Float64(z * Float64(545.0359493463282 + Float64(z * 606.6767878347069)))))))) end
function tmp = code(z) tmp = (pi * (((sqrt((2.0 * pi)) * ((7.5 - z) ^ (0.5 - z))) / exp((7.5 - z))) / sin((z * pi)))) * ((((-0.13857109526572012 / (6.0 - z)) + (9.984369578019572e-6 / (7.0 - z))) + (1.5056327351493116e-7 / (8.0 - z))) + (263.4062807184368 + (z * (436.9000215473151 + (z * (545.0359493463282 + (z * 606.6767878347069))))))); end
code[z_] := N[(N[(Pi * N[(N[(N[(N[Sqrt[N[(2.0 * Pi), $MachinePrecision]], $MachinePrecision] * N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[Sin[N[(z * Pi), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[(N[(-0.13857109526572012 / N[(6.0 - z), $MachinePrecision]), $MachinePrecision] + N[(9.984369578019572e-6 / N[(7.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(1.5056327351493116e-7 / N[(8.0 - z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(263.4062807184368 + N[(z * N[(436.9000215473151 + N[(z * N[(545.0359493463282 + N[(z * 606.6767878347069), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(\pi \cdot \frac{\frac{\sqrt{2 \cdot \pi} \cdot {\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{e^{7.5 - z}}}{\sin \left(z \cdot \pi\right)}\right) \cdot \left(\left(\left(\frac{-0.13857109526572012}{6 - z} + \frac{9.984369578019572 \cdot 10^{-6}}{7 - z}\right) + \frac{1.5056327351493116 \cdot 10^{-7}}{8 - z}\right) + \left(263.4062807184368 + z \cdot \left(436.9000215473151 + z \cdot \left(545.0359493463282 + z \cdot 606.6767878347069\right)\right)\right)\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-lowering-*.f64N/A
+-lowering-+.f64N/A
*-commutativeN/A
*-lowering-*.f6497.7%
Simplified97.7%
Final simplification97.7%
(FPCore (z)
:precision binary64
(*
(/
(+
263.3831869810514
(*
z
(+
436.8961725563396
(* z (+ 545.0353078428827 (* (* PI PI) 43.89719783017524))))))
z)
(* (exp (+ z -7.5)) (/ (pow (- 7.5 z) (- 0.5 z)) (pow (* 2.0 PI) -0.5)))))
double code(double z) {
return ((263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((((double) M_PI) * ((double) M_PI)) * 43.89719783017524)))))) / z) * (exp((z + -7.5)) * (pow((7.5 - z), (0.5 - z)) / pow((2.0 * ((double) M_PI)), -0.5)));
}
public static double code(double z) {
return ((263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((Math.PI * Math.PI) * 43.89719783017524)))))) / z) * (Math.exp((z + -7.5)) * (Math.pow((7.5 - z), (0.5 - z)) / Math.pow((2.0 * Math.PI), -0.5)));
}
def code(z): return ((263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((math.pi * math.pi) * 43.89719783017524)))))) / z) * (math.exp((z + -7.5)) * (math.pow((7.5 - z), (0.5 - z)) / math.pow((2.0 * math.pi), -0.5)))
function code(z) return Float64(Float64(Float64(263.3831869810514 + Float64(z * Float64(436.8961725563396 + Float64(z * Float64(545.0353078428827 + Float64(Float64(pi * pi) * 43.89719783017524)))))) / z) * Float64(exp(Float64(z + -7.5)) * Float64((Float64(7.5 - z) ^ Float64(0.5 - z)) / (Float64(2.0 * pi) ^ -0.5)))) end
function tmp = code(z) tmp = ((263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((pi * pi) * 43.89719783017524)))))) / z) * (exp((z + -7.5)) * (((7.5 - z) ^ (0.5 - z)) / ((2.0 * pi) ^ -0.5))); end
code[z_] := N[(N[(N[(263.3831869810514 + N[(z * N[(436.8961725563396 + N[(z * N[(545.0353078428827 + N[(N[(Pi * Pi), $MachinePrecision] * 43.89719783017524), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] * N[(N[Exp[N[(z + -7.5), $MachinePrecision]], $MachinePrecision] * N[(N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision] / N[Power[N[(2.0 * Pi), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{263.3831869810514 + z \cdot \left(436.8961725563396 + z \cdot \left(545.0353078428827 + \left(\pi \cdot \pi\right) \cdot 43.89719783017524\right)\right)}{z} \cdot \left(e^{z + -7.5} \cdot \frac{{\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{{\left(2 \cdot \pi\right)}^{-0.5}}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified96.2%
associate-*l/N/A
clear-numN/A
associate-*l/N/A
div-invN/A
times-fracN/A
*-lowering-*.f64N/A
exp-diffN/A
associate-/r/N/A
rec-expN/A
metadata-evalN/A
prod-expN/A
exp-lowering-exp.f64N/A
+-lowering-+.f64N/A
/-lowering-/.f64N/A
pow-lowering-pow.f64N/A
--lowering--.f64N/A
--lowering--.f64N/A
Applied egg-rr97.0%
Final simplification97.0%
(FPCore (z)
:precision binary64
(*
(+
263.3831869810514
(*
z
(+
436.8961725563396
(* z (+ 545.0353078428827 (* (* PI PI) 43.89719783017524))))))
(/
(/ (pow (- 7.5 z) (- 0.5 z)) (* (exp (- 7.5 z)) (pow (* 2.0 PI) -0.5)))
z)))
double code(double z) {
return (263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((((double) M_PI) * ((double) M_PI)) * 43.89719783017524)))))) * ((pow((7.5 - z), (0.5 - z)) / (exp((7.5 - z)) * pow((2.0 * ((double) M_PI)), -0.5))) / z);
}
public static double code(double z) {
return (263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((Math.PI * Math.PI) * 43.89719783017524)))))) * ((Math.pow((7.5 - z), (0.5 - z)) / (Math.exp((7.5 - z)) * Math.pow((2.0 * Math.PI), -0.5))) / z);
}
def code(z): return (263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((math.pi * math.pi) * 43.89719783017524)))))) * ((math.pow((7.5 - z), (0.5 - z)) / (math.exp((7.5 - z)) * math.pow((2.0 * math.pi), -0.5))) / z)
function code(z) return Float64(Float64(263.3831869810514 + Float64(z * Float64(436.8961725563396 + Float64(z * Float64(545.0353078428827 + Float64(Float64(pi * pi) * 43.89719783017524)))))) * Float64(Float64((Float64(7.5 - z) ^ Float64(0.5 - z)) / Float64(exp(Float64(7.5 - z)) * (Float64(2.0 * pi) ^ -0.5))) / z)) end
function tmp = code(z) tmp = (263.3831869810514 + (z * (436.8961725563396 + (z * (545.0353078428827 + ((pi * pi) * 43.89719783017524)))))) * ((((7.5 - z) ^ (0.5 - z)) / (exp((7.5 - z)) * ((2.0 * pi) ^ -0.5))) / z); end
code[z_] := N[(N[(263.3831869810514 + N[(z * N[(436.8961725563396 + N[(z * N[(545.0353078428827 + N[(N[(Pi * Pi), $MachinePrecision] * 43.89719783017524), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[(N[Power[N[(7.5 - z), $MachinePrecision], N[(0.5 - z), $MachinePrecision]], $MachinePrecision] / N[(N[Exp[N[(7.5 - z), $MachinePrecision]], $MachinePrecision] * N[Power[N[(2.0 * Pi), $MachinePrecision], -0.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(263.3831869810514 + z \cdot \left(436.8961725563396 + z \cdot \left(545.0353078428827 + \left(\pi \cdot \pi\right) \cdot 43.89719783017524\right)\right)\right) \cdot \frac{\frac{{\left(7.5 - z\right)}^{\left(0.5 - z\right)}}{e^{7.5 - z} \cdot {\left(2 \cdot \pi\right)}^{-0.5}}}{z}
\end{array}
Initial program 97.0%
Simplified97.7%
Applied egg-rr99.1%
Taylor expanded in z around 0
/-lowering-/.f64N/A
Simplified96.2%
Applied egg-rr97.0%
(FPCore (z) :precision binary64 (* 263.3831869810514 (/ (* (sqrt 15.0) (/ (sqrt PI) (exp 7.5))) z)))
double code(double z) {
return 263.3831869810514 * ((sqrt(15.0) * (sqrt(((double) M_PI)) / exp(7.5))) / z);
}
public static double code(double z) {
return 263.3831869810514 * ((Math.sqrt(15.0) * (Math.sqrt(Math.PI) / Math.exp(7.5))) / z);
}
def code(z): return 263.3831869810514 * ((math.sqrt(15.0) * (math.sqrt(math.pi) / math.exp(7.5))) / z)
function code(z) return Float64(263.3831869810514 * Float64(Float64(sqrt(15.0) * Float64(sqrt(pi) / exp(7.5))) / z)) end
function tmp = code(z) tmp = 263.3831869810514 * ((sqrt(15.0) * (sqrt(pi) / exp(7.5))) / z); end
code[z_] := N[(263.3831869810514 * N[(N[(N[Sqrt[15.0], $MachinePrecision] * N[(N[Sqrt[Pi], $MachinePrecision] / N[Exp[7.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
263.3831869810514 \cdot \frac{\sqrt{15} \cdot \frac{\sqrt{\pi}}{e^{7.5}}}{z}
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
times-fracN/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
metadata-evalN/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
exp-lowering-exp.f6496.9%
Applied egg-rr96.9%
(FPCore (z) :precision binary64 (* 263.3831869810514 (* (/ (sqrt PI) (exp 7.5)) (/ (sqrt 15.0) z))))
double code(double z) {
return 263.3831869810514 * ((sqrt(((double) M_PI)) / exp(7.5)) * (sqrt(15.0) / z));
}
public static double code(double z) {
return 263.3831869810514 * ((Math.sqrt(Math.PI) / Math.exp(7.5)) * (Math.sqrt(15.0) / z));
}
def code(z): return 263.3831869810514 * ((math.sqrt(math.pi) / math.exp(7.5)) * (math.sqrt(15.0) / z))
function code(z) return Float64(263.3831869810514 * Float64(Float64(sqrt(pi) / exp(7.5)) * Float64(sqrt(15.0) / z))) end
function tmp = code(z) tmp = 263.3831869810514 * ((sqrt(pi) / exp(7.5)) * (sqrt(15.0) / z)); end
code[z_] := N[(263.3831869810514 * N[(N[(N[Sqrt[Pi], $MachinePrecision] / N[Exp[7.5], $MachinePrecision]), $MachinePrecision] * N[(N[Sqrt[15.0], $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
263.3831869810514 \cdot \left(\frac{\sqrt{\pi}}{e^{7.5}} \cdot \frac{\sqrt{15}}{z}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
times-fracN/A
*-commutativeN/A
*-lowering-*.f64N/A
/-lowering-/.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
exp-lowering-exp.f64N/A
/-lowering-/.f64N/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
metadata-eval96.7%
Applied egg-rr96.7%
(FPCore (z) :precision binary64 (* (/ (* 263.3831869810514 (pow (* PI 15.0) 0.5)) z) (exp -7.5)))
double code(double z) {
return ((263.3831869810514 * pow((((double) M_PI) * 15.0), 0.5)) / z) * exp(-7.5);
}
public static double code(double z) {
return ((263.3831869810514 * Math.pow((Math.PI * 15.0), 0.5)) / z) * Math.exp(-7.5);
}
def code(z): return ((263.3831869810514 * math.pow((math.pi * 15.0), 0.5)) / z) * math.exp(-7.5)
function code(z) return Float64(Float64(Float64(263.3831869810514 * (Float64(pi * 15.0) ^ 0.5)) / z) * exp(-7.5)) end
function tmp = code(z) tmp = ((263.3831869810514 * ((pi * 15.0) ^ 0.5)) / z) * exp(-7.5); end
code[z_] := N[(N[(N[(263.3831869810514 * N[Power[N[(Pi * 15.0), $MachinePrecision], 0.5], $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision] * N[Exp[-7.5], $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{263.3831869810514 \cdot {\left(\pi \cdot 15\right)}^{0.5}}{z} \cdot e^{-7.5}
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
rec-expN/A
exp-lowering-exp.f64N/A
metadata-evalN/A
sqrt-unprodN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
PI-lowering-PI.f6496.2%
Applied egg-rr96.2%
*-commutativeN/A
clear-numN/A
un-div-invN/A
sqrt-prodN/A
div-invN/A
rec-expN/A
metadata-evalN/A
frac-timesN/A
*-commutativeN/A
associate-*l/N/A
div-invN/A
rec-expN/A
metadata-evalN/A
*-lowering-*.f64N/A
Applied egg-rr96.1%
associate-*r*N/A
*-lowering-*.f64N/A
associate-*r/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
pow1/2N/A
pow-lowering-pow.f64N/A
*-lowering-*.f64N/A
PI-lowering-PI.f64N/A
exp-lowering-exp.f6496.4%
Applied egg-rr96.4%
(FPCore (z) :precision binary64 (/ (* 263.3831869810514 (sqrt (* PI 15.0))) (/ z (exp -7.5))))
double code(double z) {
return (263.3831869810514 * sqrt((((double) M_PI) * 15.0))) / (z / exp(-7.5));
}
public static double code(double z) {
return (263.3831869810514 * Math.sqrt((Math.PI * 15.0))) / (z / Math.exp(-7.5));
}
def code(z): return (263.3831869810514 * math.sqrt((math.pi * 15.0))) / (z / math.exp(-7.5))
function code(z) return Float64(Float64(263.3831869810514 * sqrt(Float64(pi * 15.0))) / Float64(z / exp(-7.5))) end
function tmp = code(z) tmp = (263.3831869810514 * sqrt((pi * 15.0))) / (z / exp(-7.5)); end
code[z_] := N[(N[(263.3831869810514 * N[Sqrt[N[(Pi * 15.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / N[(z / N[Exp[-7.5], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\frac{263.3831869810514 \cdot \sqrt{\pi \cdot 15}}{\frac{z}{e^{-7.5}}}
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
rec-expN/A
exp-lowering-exp.f64N/A
metadata-evalN/A
sqrt-unprodN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
PI-lowering-PI.f6496.2%
Applied egg-rr96.2%
*-commutativeN/A
*-commutativeN/A
clear-numN/A
un-div-invN/A
sqrt-prodN/A
div-invN/A
rec-expN/A
metadata-evalN/A
frac-timesN/A
*-commutativeN/A
*-commutativeN/A
frac-timesN/A
sqrt-prodN/A
*-commutativeN/A
associate-*l/N/A
/-lowering-/.f64N/A
Applied egg-rr96.4%
Final simplification96.4%
(FPCore (z) :precision binary64 (* (* 263.3831869810514 (sqrt (* PI 15.0))) (/ (exp -7.5) z)))
double code(double z) {
return (263.3831869810514 * sqrt((((double) M_PI) * 15.0))) * (exp(-7.5) / z);
}
public static double code(double z) {
return (263.3831869810514 * Math.sqrt((Math.PI * 15.0))) * (Math.exp(-7.5) / z);
}
def code(z): return (263.3831869810514 * math.sqrt((math.pi * 15.0))) * (math.exp(-7.5) / z)
function code(z) return Float64(Float64(263.3831869810514 * sqrt(Float64(pi * 15.0))) * Float64(exp(-7.5) / z)) end
function tmp = code(z) tmp = (263.3831869810514 * sqrt((pi * 15.0))) * (exp(-7.5) / z); end
code[z_] := N[(N[(263.3831869810514 * N[Sqrt[N[(Pi * 15.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] * N[(N[Exp[-7.5], $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
\left(263.3831869810514 \cdot \sqrt{\pi \cdot 15}\right) \cdot \frac{e^{-7.5}}{z}
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
associate-*r/N/A
div-invN/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-unprodN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
PI-lowering-PI.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
rec-expN/A
exp-lowering-exp.f64N/A
metadata-eval96.3%
Applied egg-rr96.3%
Final simplification96.3%
(FPCore (z) :precision binary64 (* 263.3831869810514 (/ (* (exp -7.5) (sqrt (* PI 15.0))) z)))
double code(double z) {
return 263.3831869810514 * ((exp(-7.5) * sqrt((((double) M_PI) * 15.0))) / z);
}
public static double code(double z) {
return 263.3831869810514 * ((Math.exp(-7.5) * Math.sqrt((Math.PI * 15.0))) / z);
}
def code(z): return 263.3831869810514 * ((math.exp(-7.5) * math.sqrt((math.pi * 15.0))) / z)
function code(z) return Float64(263.3831869810514 * Float64(Float64(exp(-7.5) * sqrt(Float64(pi * 15.0))) / z)) end
function tmp = code(z) tmp = 263.3831869810514 * ((exp(-7.5) * sqrt((pi * 15.0))) / z); end
code[z_] := N[(263.3831869810514 * N[(N[(N[Exp[-7.5], $MachinePrecision] * N[Sqrt[N[(Pi * 15.0), $MachinePrecision]], $MachinePrecision]), $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
263.3831869810514 \cdot \frac{e^{-7.5} \cdot \sqrt{\pi \cdot 15}}{z}
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
rec-expN/A
exp-lowering-exp.f64N/A
metadata-evalN/A
sqrt-unprodN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
PI-lowering-PI.f6496.2%
Applied egg-rr96.2%
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f64N/A
sqrt-lowering-sqrt.f64N/A
*-commutativeN/A
*-lowering-*.f64N/A
PI-lowering-PI.f6496.2%
Applied egg-rr96.2%
(FPCore (z) :precision binary64 (* 263.3831869810514 (* (sqrt (* PI 15.0)) (/ (exp -7.5) z))))
double code(double z) {
return 263.3831869810514 * (sqrt((((double) M_PI) * 15.0)) * (exp(-7.5) / z));
}
public static double code(double z) {
return 263.3831869810514 * (Math.sqrt((Math.PI * 15.0)) * (Math.exp(-7.5) / z));
}
def code(z): return 263.3831869810514 * (math.sqrt((math.pi * 15.0)) * (math.exp(-7.5) / z))
function code(z) return Float64(263.3831869810514 * Float64(sqrt(Float64(pi * 15.0)) * Float64(exp(-7.5) / z))) end
function tmp = code(z) tmp = 263.3831869810514 * (sqrt((pi * 15.0)) * (exp(-7.5) / z)); end
code[z_] := N[(263.3831869810514 * N[(N[Sqrt[N[(Pi * 15.0), $MachinePrecision]], $MachinePrecision] * N[(N[Exp[-7.5], $MachinePrecision] / z), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]
\begin{array}{l}
\\
263.3831869810514 \cdot \left(\sqrt{\pi \cdot 15} \cdot \frac{e^{-7.5}}{z}\right)
\end{array}
Initial program 97.0%
Simplified97.7%
Taylor expanded in z around 0
*-lowering-*.f64N/A
associate-*l/N/A
/-lowering-/.f64N/A
*-lowering-*.f64N/A
*-lowering-*.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
sqrt-lowering-sqrt.f64N/A
PI-lowering-PI.f64N/A
*-lowering-*.f64N/A
exp-lowering-exp.f6495.7%
Simplified95.7%
clear-numN/A
associate-/r/N/A
*-lowering-*.f64N/A
*-commutativeN/A
associate-/r*N/A
/-lowering-/.f64N/A
rec-expN/A
exp-lowering-exp.f64N/A
metadata-evalN/A
sqrt-unprodN/A
sqrt-unprodN/A
sqrt-lowering-sqrt.f64N/A
*-lowering-*.f64N/A
metadata-evalN/A
PI-lowering-PI.f6496.2%
Applied egg-rr96.2%
Final simplification96.2%
herbie shell --seed 2024164
(FPCore (z)
:name "Jmat.Real.gamma, branch z less than 0.5"
:precision binary64
:pre (<= z 0.5)
(* (/ PI (sin (* PI z))) (* (* (* (sqrt (* PI 2.0)) (pow (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5) (+ (- (- 1.0 z) 1.0) 0.5))) (exp (- (+ (+ (- (- 1.0 z) 1.0) 7.0) 0.5)))) (+ (+ (+ (+ (+ (+ (+ (+ 0.9999999999998099 (/ 676.5203681218851 (+ (- (- 1.0 z) 1.0) 1.0))) (/ -1259.1392167224028 (+ (- (- 1.0 z) 1.0) 2.0))) (/ 771.3234287776531 (+ (- (- 1.0 z) 1.0) 3.0))) (/ -176.6150291621406 (+ (- (- 1.0 z) 1.0) 4.0))) (/ 12.507343278686905 (+ (- (- 1.0 z) 1.0) 5.0))) (/ -0.13857109526572012 (+ (- (- 1.0 z) 1.0) 6.0))) (/ 9.984369578019572e-6 (+ (- (- 1.0 z) 1.0) 7.0))) (/ 1.5056327351493116e-7 (+ (- (- 1.0 z) 1.0) 8.0))))))