In this paper a fifth-order numerical scheme is developed and implemented for the solution of homogeneous heat equation ut = a uxx with nonlocal boundary condition as well as for inhomogeneous heat equation ut = a uxx + s ( x,t ) with nonlocal boundary condition. The results obtained show that the numerical method based on the proposed technique is fifth-order accurate as well as L -acceptable. In the development of this method second-order spatial derivative are approximated by fifth-order finite-difference approximations which give a system of first order, linear, ordinary differential equations whose solution satisfies a recurrence relation which leads to the development of algorithm. The algorithm is tested on various heat equations and no oscillations are observed in the experiments. This method is based on partial fraction technique which is useful in parrel processing and it does not require complex arithmetic.