Fractional approximations for linear first‐order differential equations with polynomial coefficients—application to E1(x)