I know that $\operatorname{IV-model \space free}=2 \int_{0}^{+\infty}\frac{c_0(T,Ke^{r(T-t)})-c_0(t,Ke^{r(T-t)})}{K^2}\operatorname{d}K$ is calculated using an iterative procedure, i.e. setting a volatility value and comparing theoretical price obtained by software (with that volatility value) and the market price: if theoretical price is less than market price we have to increase the volatility (and vice-versa). This would be the reason:
"You can not invert the BS-equation because implied volatility is repeated both in 1) the PDF of Normal standard like standard deviation of log-returns distribution and 2) the argument of Normal standard in $d_1$ and $d_2$ terms."
Deriving BS-close form I have $$\varphi(S_t)=e^{-r(T-t)}\int_{z_0}^{+\infty}((S_te^{(r-\frac{\sigma^2}{2})(T-t)+\sigma\sqrt{T-t}z})-K)^+\Phi(z)\operatorname{d}z$$ with $z_0=\frac{\operatorname{ln}(\frac{S_t}{K})+(r+\frac{\sigma^2}{2})(T-t)}{\sigma\sqrt{T-t}}$ and $\Phi(z)\doteq \frac{1}{\sqrt{2\pi}}e^{-\frac{z^2}{2}}$. So, knowing that $d_1:=-z_0+\sigma\sqrt{T-t}$ and $d_2:=d_1-\sigma\sqrt{T-t}$ (so the point 2) is trivially true), what is the $\sigma=\operatorname{IV}$ to which it refers the point 1)?
Thanks in advance for any help!