Say we have an Electric Field, produced by a charge distribution, given as- \begin{equation} \mathbf{E}=c(1-e^{-\alpha r}) \frac{\hat{\mathbf{r}}}{r^2}, \end{equation} $c$ and $\alpha$ being constants.
Now if we are asked to find net charge within a radius $r=\frac{1}{\alpha}$, how should we approach?
My thought is that since the charge distribution is unbounded, we can't have a Gaussian surface that contains all the charge, and so we can't just trivially solve it using Gauss's law \begin{equation} \int{\mathbf{E}}\cdot{d\mathbf{S}}=\frac{Q}{\epsilon_0} \end{equation}
Even if we go forward with the old method using Gauss's Law, the answer comes out wrong.
What should we do? How should we approach?