A high-order transport approximation, the SKN method, is introduced for solving the integral transport equation. The method relies on approximating the integral transport kernels by a sum of diffusionlike kernels. The integral equation is equivalent to a set of coupled second-order differential equations for which blackbody and reflecting boundary conditions are established. These SKN equations are solved for benchmark problems. The benchmark problems include one- and two-dimensional homogeneous cell configurations. The solutions using SKN are compared with discrete ordinates and other high-order transport theory solutions.