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(10000 \cdot \left(y + 1\right)\right)} = 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(10000 \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(10000 \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(10000 \cdot \left(y + 1\right)\right)} = 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(10000 \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(10000 \cdot \left(y + 1\right)\right)}}\\

\end{array}
double f(double y) {
        double r1331030 = y;
        double r1331031 = r1331030 * r1331030;
        double r1331032 = 1.0;
        double r1331033 = r1331031 + r1331032;
        double r1331034 = sqrt(r1331033);
        double r1331035 = r1331030 - r1331034;
        double r1331036 = fabs(r1331035);
        double r1331037 = r1331030 + r1331034;
        double r1331038 = r1331032 / r1331037;
        double r1331039 = r1331036 - r1331038;
        double r1331040 = r1331039 * r1331039;
        double r1331041 = 10.0;
        double r1331042 = -300.0;
        double r1331043 = pow(r1331041, r1331042);
        double r1331044 = 10000.0;
        double r1331045 = r1331030 + r1331032;
        double r1331046 = r1331044 * r1331045;
        double r1331047 = pow(r1331043, r1331046);
        double r1331048 = r1331040 + r1331047;
        double r1331049 = 0.0;
        double r1331050 = r1331048 == r1331049;
        double r1331051 = exp(r1331048);
        double r1331052 = r1331051 - r1331032;
        double r1331053 = r1331052 / r1331048;
        double r1331054 = r1331050 ? r1331032 : r1331053;
        return r1331054;
}