Update TothStationaryDensity.R