By superposition the field $$\vec {E}$$ inside the ball is given by
$$\vec {E} = \vec {E_{0}} - \dfrac {\vec {P}}{3\epsilon_{0}}$$
On the other hand, if the sphere is not too small, the macroscopic equation $$\vec {P} = (\epsilon - 1) \epsilon_{\theta} \vec {E}$$ must hold. Thus,
$$\vec {E} \left (1 + \dfrac {1}{3} (\epsilon - 1)\right ) = \vec {E_{0}}$$ or, $$\vec {E} = \dfrac {3\vec {E_{0}}}{\epsilon + 2}$$
Also $$\vec {P} = 3\epsilon_{0} \dfrac {\epsilon - 1}{\epsilon + 2} \vec {E}_{0}$$.
We write $$\vec {E} = \vec {E_{0}} - \dfrac {\vec {P}}{2\epsilon_{0}}$$
using here the result of the foregoing problem.
Also $$\vec {P} = (\epsilon - 1) \epsilon_{0} \vec {E}$$
So, $$\vec {E} \left (\dfrac {\epsilon + 1}{2}\right ) = \vec {E}_{0}$$, or, $$\vec {E} = \dfrac {2\vec {E}_{0}}{\epsilon + 1}$$ and $$\vec {P} = 2\epsilon_{0} \dfrac {\epsilon - 1}{\epsilon + 1} \vec {E}_{0}$$.