In the teaching of annular fins of rectangular profile, the main obstacle lies in solving the quasi one-dimensional heat conduction equation, the modified Bessel equation of first kind. In the modified Bessel equation, the variable coefficient
multiplying the first order derivative of temperature
is problematic. The principal objective of the present paper on engineering education is to solve the modified Bessel equation of first kind in approximate analytical form. Specifically, we seek to apply the mean value theorem for integrals to the variable coefficient
, viewed as an auxiliary function in the domain of the annular fin extending from the inner radius
to the outer radius
. It is demonstrated in a convincing manner that approximate analytical temperature distributions having exponential functions are easy to obtain without resorting to the exact analytical temperature distribution embodying four modified Bessel functions of first kind. Furthermore, the easiness in calculating heat transfer rates in annular fins of rectangular profile for realistic combinations of the two controlling parameters: the normalized radius ratio and the dimensionless thermo-geometrical parameter is verifiable.