An analytical, two‐dimensional, stationary model of heat flow in edge‐emitting laser is revisited. In the work, we show how to use this approach to be able to calculate the temperature of the entire device including the most susceptible for thermal runaway region, namely the vicinity of the joint of mirrors and an active layer. Numerical tools based on our considerations are implemented in Matlab code and available at the software developer's web site. A procedure of investigation of surface recombination at facets with the help of our model combined with mirror temperature measurements is proposed. To get an insight in reasonableness of our theoretical research, we calculated temperature profiles and surface recombination velocities for example devices and successfully confronted them with results found in the literature.