gradual adjustment of R

seems to improve stability, converged for phi = 0.01