We should probably proceed by cases, as $f$ will have a vertical asymptote and a hole for $c=0,$ a single vertical asymptote and no holes for $c=1,$ neither vertical asymptotes nor holes when $c>1,$ and otherwise has two vertical asymptotes and no holes. In all cases, $f$ has $y=1$ as its horizontal asymptote, meaning that $\lim_{|x|\to\infty}f(x)=1.$
In the $c>1$ case, differentiating yields $$f'(x)=\frac{x^2-c}{(x^2+2x+c)^2},$$ which is negative in the interval $\left(-\sqrt{c},\sqrt{c}\right),$ zero at $x=\pm\sqrt{c},$ and positive otherwise. Consequently, $f$ achieves a global maximum at $x=-\sqrt{c}$ and a global minimum at $x=\sqrt{c}.$ The minimum is the one we care about, though, since we'll require it to be greater than $-1.$ So, let's consider $$-1 But this is trivially true, since $c>1$ implies that both $2\sqrt c+1$ and $2\sqrt c+2$ are positive, so that $f(c)>0>-1,$ as desired.
I leave the $c=0$ and $c=1$ cases to you.
If $c<1$ and $c\ne 0,$ we know that $x^2+2x+c=(x-j)(x-k)$ for some distinct real numbers $j$ and $k.$ We also know that $x^2+x+c$ takes on a non-zero value at each of $j,k.$ Without loss of generality, suppose that $j Now, at $x=j,$ the denominator changes signs, while the numerator does not, meaning that we have either $$\lim_{x\nearrow j}f(x)=-\infty\text{ and }\lim_{x\searrow j}f(x)=+\infty\tag{1}$$ or $$\lim_{x\nearrow j}f(x)=+\infty\text{ and }\lim_{x\searrow j}f(x)=-\infty.\tag{2}$$ Likewise, we have either $$\lim_{x\nearrow k}f(x)=-\infty\text{ and }\lim_{x\searrow k}f(x)=+\infty\tag{3}$$ or $$\lim_{x\nearrow k}f(x)=+\infty\text{ and }\lim_{x\searrow k}f(x)=-\infty.\tag{4}$$ If $(1)$ holds, then since $f$ is continuous on $(-\infty,j)$ and since $\lim_{x\to-\infty}f(x)=1$ and $\lim_{x\nearrow j}f(x)=-\infty,$ then $(-\infty,1)$ is a subset of the range of $f,$ and so $\left[-1,-\frac13\right]$ is, too. If $(4)$ holds, then by continuity of $f$ on $(k,\infty),$ we similarly have that $\left[-1,-\frac13\right]$ is a subset of the range of $f.$ However, if $(2)$ and $(3)$ both hold, it isn't clear from what we've discussed so far whether $\left[-1,-\frac13\right]$ is a subset of the range of $f$ or not, so let's consider further when these different possibilities may occur.
The only way for $(1)$ to hold is if the numerator is negative at $x=j,$ meaning that $j$ must lie between the two zeroes of $x^2+x+c.$ In other words, by the quadratic formula, we have the following (equivalent) inequalities: $$\frac{-1-\sqrt{1^2-4c}}2<\frac{-2-\sqrt{2^2-4c}}2<\frac{-1+\sqrt{1^2-4c}}2\\-1-\sqrt{1-4c}<-2-\sqrt{4-4c}<-1+\sqrt{1-4c}\\-\sqrt{1-4c}<-1-\sqrt{4-4c}<\sqrt{1-4c}\\\left(-1-\sqrt{4-4c}\right)^2<1-4c\\1+2\sqrt{4-4c}+4-4c<1-4c\\2\sqrt{4-4c}+4<0.$$ Since this is impossible, then $(1)$ cannot hold, and so $(2)$ must hold.
Similarly, $(4)$ only holds if the numerator is negative at $x=k,$ which happens if and only if $$\frac{-1-\sqrt{1^2-4c}}2<\frac{-2+\sqrt{2^2-4c}}2<\frac{-1+\sqrt{1^2-4c}}2\\-1-\sqrt{1-4c}<-2+\sqrt{4-4c}<-1+\sqrt{1-4c}\\-\sqrt{1-4c}<-1+\sqrt{4-4c}<\sqrt{1-4c}\\\left(-1+\sqrt{4-4c}\right)^2<1-4c\\1-2\sqrt{4-4c}+4-4c<1-4c\\-2\sqrt{4-4c}+4<0\\2<\sqrt{4-4c}\\4<4-4c\\c<0.$$ Thus, if $c<0,$ then $\left[-1,-\frac13\right]$ is a subset of the range of $f.$
It remains only to consider $0 in which case $(2)$ and $(3)$ will hold, so that $f$ achieves a local maximum in the interval $(j,k).$ Observe that $$1-c<1\\\sqrt{1-c}<1\\2\sqrt{1-c}<2\\\sqrt{4-4c}<2\\-2+\sqrt{4-4c}<0\\\frac{-2+\sqrt{4-4c}}2<0\\k<0,$$ meaning that $f$ achieves its local maximum at $x=-\sqrt{c},$ and a local minimum at $x=\sqrt{c}$ by the first derivative test. Since $f$ is increasing and continuous on $(-\infty,j),$ and since $\lim_{x\to-\infty}f(x)=1,$ then $f$ is positive on $(-\infty,j).$ Further, since $f\bigl(\sqrt c\bigr)$ is positive by prior work, then since $f$'s minimum value in $(k,\infty)$ occurs at $x=\sqrt c,$ then we have that $f$ is positive on $(k,\infty).$ Thus, we need only ensure that $f\left(-\sqrt c\right)<-\frac13.$ Since $0 then $c<\sqrt c,$ and so the following are equivalent: $$f\left(-\sqrt c\right)<-\frac13\\-3f\left(-\sqrt c\right)>1\\-3\cdot\frac{c-\sqrt c+c}{c-2\sqrt c+c}>1\\-3\cdot\frac{2c-\sqrt c}{2c-2\sqrt c}>1\\-3\left(2c-\sqrt c\right)<2c-2\sqrt c\\-6c+3\sqrt c<2c-2\sqrt c\\5\sqrt c<8c\\\frac58<\sqrt c\\\frac{25}{64}
Putting it all together, we find that the range of $f$ does not contain the interval in question if and only if $c>\frac{25}{64}.$