Linear profile multilayer approximation for computing reflection or transmission coefficients of plane wave by a layered medium, which is formal extension of standard staircase stratification method, is studied theoretically and numerically. The problem is reduced to evaluate

matrix product whose elements include the Bessel functions of order

and

. The polynomial approximation for the Bessel function is constructed by applying modified Chebyshev method to facilitate the numerical calculation.