Update TothStationaryDensity.R

Latest changes, working code with plots for book density ratios against u