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 r21756 = y;
        double r21757 = r21756 * r21756;
        double r21758 = 1.0;
        double r21759 = r21757 + r21758;
        double r21760 = sqrt(r21759);
        double r21761 = r21756 - r21760;
        double r21762 = fabs(r21761);
        double r21763 = r21756 + r21760;
        double r21764 = r21758 / r21763;
        double r21765 = r21762 - r21764;
        double r21766 = r21765 * r21765;
        double r21767 = 10.0;
        double r21768 = -300.0;
        double r21769 = pow(r21767, r21768);
        double r21770 = 10000.0;
        double r21771 = r21756 + r21758;
        double r21772 = r21770 * r21771;
        double r21773 = pow(r21769, r21772);
        double r21774 = r21766 + r21773;
        double r21775 = 0.0;
        double r21776 = r21774 == r21775;
        double r21777 = exp(r21774);
        double r21778 = r21777 - r21758;
        double r21779 = r21778 / r21774;
        double r21780 = r21776 ? r21758 : r21779;
        return r21780;
}