The thermal response of a tunable laser is analyzed by using a mode density method based on a Fourier-Laplace analysis. This method introduces a mode density function for mode distribution of the Fourier-Laplace transform and gives temperature time-dependency in an integral form instead of an infinite weighted summation. When symmetric structures are assumed, the mode density method gives the transient thermal response in a simple form: error functions (spherical-symmetry case) and exponential integral functions (cylindrical-symmetry case). The cylindrical-symmetry analysis was extended to the noncylindrical-symmetry model and the thermal response of the tunable laser was calculated by the mode density method. The result shows good agreement with a Fourier-Laplace analysis (deviation 2%) and experimental results. As a rough estimation, the thermal response of the laser is in proportion to the logarithm of time in some range that depends on the chip and tuning-section size of the laser.