We deal with one-dimensional source inverse heat conduction equation. An approach based on least squares support vector machines (LS-SVM) is proposed for semi-analytic approximate solutions. Furthermore, a parameters tuning formulism is shown and stability of the method is presented. The method yields high accuracy and stability solutions in practical examples.