Cannot sample enough valid points. (more)

\[1 \le y \le 9999\]
\[\begin{array}{l} \mathbf{if}\;\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)} = 0.0:\\ \;\;\;\;1\\ \mathbf{else}:\\ \;\;\;\;\frac{e^{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)}} - 1}{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)}}\\ \end{array}\]
\begin{array}{l}
\mathbf{if}\;\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)} = 0.0:\\
\;\;\;\;1\\

\mathbf{else}:\\
\;\;\;\;\frac{e^{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)}} - 1}{\left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) \cdot \left(\left|y - \sqrt{y \cdot y + 1}\right| - \frac{1}{y + \sqrt{y \cdot y + 1}}\right) + {\left({10}^{-300}\right)}^{\left(10^{4} \cdot \left(y + 1\right)\right)}}\\

\end{array}
double f(double y) {
        double r1000641 = y;
        double r1000642 = r1000641 * r1000641;
        double r1000643 = 1.0;
        double r1000644 = r1000642 + r1000643;
        double r1000645 = sqrt(r1000644);
        double r1000646 = r1000641 - r1000645;
        double r1000647 = fabs(r1000646);
        double r1000648 = r1000641 + r1000645;
        double r1000649 = r1000643 / r1000648;
        double r1000650 = r1000647 - r1000649;
        double r1000651 = r1000650 * r1000650;
        double r1000652 = 10.0;
        double r1000653 = -300.0;
        double r1000654 = pow(r1000652, r1000653);
        double r1000655 = 10000.0;
        double r1000656 = r1000641 + r1000643;
        double r1000657 = r1000655 * r1000656;
        double r1000658 = pow(r1000654, r1000657);
        double r1000659 = r1000651 + r1000658;
        double r1000660 = 0.0;
        double r1000661 = r1000659 == r1000660;
        double r1000662 = exp(r1000659);
        double r1000663 = r1000662 - r1000643;
        double r1000664 = r1000663 / r1000659;
        double r1000665 = r1000661 ? r1000643 : r1000664;
        return r1000665;
}