I have already calculated the polarisation that is
$$ \mathbf{P} = \frac{\rho_f r}{2} \left( 1 - \frac{\epsilon_0}{\epsilon} \right) \hat{r} . $$
I tried to use the following formulas to calculate the density bound charges. For the surface bound charge I got:
$$ \sigma_{b1} = \mathbf{P} \cdot...