added C++ code of steps 3.a and 3.b;
corrected and sorted references;
added Theorem about  estimation of $\psi(n)$;