The stationary energy current through crystalline lattice in contact with external heat reservoirs is analyzed as the model of the lattice thermal conduction. The reservoir is chozen as the pair of the white noise random force and the frictional one. This reservoir yields exact solutions for harmonic regular lattice of one, two and three dimension: the results is naturally without temperature gradient, but the contact resistance limits the heat current to a finite value. The same formulation is applied to anharmonic linear chains, and a procedure of numerical integration is presented. Executions on the computer are shown for three cases of regular anharmonic lattice. The results are in accord with the usual notion that the stronger the anharmonicity the larger the thermal resistance.