В статье описывается алгоритм вычисления матричной экспоненты, основанный на представлении интерполяционного многочлена в виде линейной комбинации ортогональных функций. Численные эксперименты показывают, что этот алгоритм в некоторых случаях приводит к большей точности, чем прямое использование интерполяционных
многочленов в форме Лагранжа и Ньютона.