This paper presents an algorithm combining non-linear programming techniques with methods of solving initial value problems to solve a non-linear boundary value problem describing the temperature field in a parallel-flow multi-channel heat exchanger. The method is illustrated by an example of determining temperature profiles in a four-channel heat exchanger. The calculations are performed with the aid of the computer program Mathcad® 2000 Professional.