In this work the combined finite difference and orthogonal function approximation methods have been proposed for the numerical solution of the one-dimensional parabolic equation with two integrals conditions. The time variable is approximated using a finite difference scheme. But the orthogonal function approximation is employed for discretizing the space variable by using orthogonal cosine function.