Abstract: C-Bézier and H-Bézier basis functions are applied to the space of trial function and test function in the finite element method and the errors are analyzed. By means of numerical examples, the error accuracy of numerical solution is improved by one to three orders of magnitude compared with Lagrange basis function, it is shown that these two kinds of basis functions have better approximations in simulating specific heat conduction problems.