Let S: [0, 1]→[0, 1] be a chaotic map and let f* be a stationary density of the Frobenius-Perron operator PS: L1→L1 associated with S. We develop a numerical algorithm for approximating f*, using the maximum entropy approach to an under-determined moment problem and the Chebyshev polynomials for the stability consideration. Numerical experiments show considerable improvements to both the original maximum entropy method and the discrete maximum entropy method.