We apply the cubic spline projection method to compute a stationary density of the Foias operator for a random map to improve the convergence rate of the linear spline and quadratic spline projection methods. We analyze the convergence of the numerical scheme as well as provide numerical experiments that support our claims.