In this article a numerical solution of non-isothermal helical flow for the thermodependent Herschel-Bulkley fluids has been presented. The consistency term (mu(o)) in comparison to other properties is more temperature dependent; therefore an exponential function in the form of mu(o) = mu exp(-bT) is considered. The governing equations are solved for two different physical situations; constant temperature walls and for a case of constant heat flux using the implicit finite difference method. The results from numerical solution are compared with results of first law of thermodynamics and a good agreement is noticed.