While repenner's answer is 100% correct, I think his method is very slightly wrong. He calculates the fraction of the sky blocked by small spherical shell then integrates that to show it is approaches 100% blocked as the "max radius, S" (or upper limit of the integral) becomes ever larger. The error is one of "over counting." I. e. some of the blocking made by shell n will also be blocked by shell n-1 stars, (Two stars perfectly alined with earth radius) etc. His integral converges to infinity a little more slowly, if this "double blocking" were not included, counted at least twice.
I believe if you would check again, I have already addressed that, because I addressed the fraction of the sky which is DARK.
Obviously, you do get an overcounting problem if you try and add fractions of the sky covered by each shell because $$\lim \limits_{S\to\infty} FS > 1$$. How can you have more than 1 full sky? But, as with other probability arguments, focusing on the opposite outcome -- the fraction of the shell not blocked by stars -- is the way to go.
Thus each dR shell has a fraction of sky with stars, F dR, and a portion of sky without stars, 1 - F dR.
So we want to multiply the "dark" part of all the shells. If $$ dR = \frac{S}{n} $$ then we want $$\lim \limits_{n\to\infty} \prod \limits_{k=1}^{n} ( 1 - F \frac{S}{n} )$$ which is almost like an integral: $$\int \limits_{0}^{S} g(R) dR = \lim \limits_{n\to\infty} \sum \limits_{k=1}^{n} g \left( \frac{k - \frac{1}{2}}{n} S \right) \times \frac{S}{n} $$. But we can made this analogy precise by replacing $$\prod g \left( \frac{k - \frac{1}{2}}{n} S \right) $$ with $$\sum \ln g \left( \frac{k - \frac{1}{2}}{n} S \right) $$. So what is $$\ln \left( 1 - F dR \right) $$ ?
The natural logarithm of this latter quantity is, -F dR. This is because $$\lim \limits_{x \to 0} \ln (1 + x) = x + O(x^2)$$.
Thus the logarithm of free sky of k shells is -k F dR and for all shells out to distance S is -F S.
Say $$F = 10^{-7} \, \textrm{pc}^{-1}$$ (ridiculously high, but useful for illustration). Compare the two models for fraction of sky covered by stars:
$$ \begin{array}{l|ll} S & FS & 1 - e^{-FS} \\ \hline
\\ 100,000 \textrm{pc} & 0.01 & 0.009950166
\\ 200,000 \textrm{pc} & 0.02 & 0.01980133
\\ 500,000 \textrm{pc} & 0.05 & 0.04877058
\\ 1 \textrm{Mpc} & 0.1 & 0.09516258
\\ 2 \textrm{Mpc} & 0.2 & 0.1812692
\\ 5 \textrm{Mpc} & 0.5 & 0.3934693
\\ 10 \textrm{Mpc} & 1 & 0.6321206
\\ 20 \textrm{Mpc} & 2 & 0.8646647
\\ 50 \textrm{Mpc} & 5 & 0.9932621
\\ 100 \textrm{Mpc} & 10 & 0.9999546 \end{array}$$
// Edited to change f -> g for the general function in illustration of integration because my previous post used f in a problem-specific manner