Normalization of the wavefunction obtained from perturbation theory based on a matrix method

We present the derivation of the normalization constant for the perturbation matrix method recently proposed. The method is tested on the problem of a binary waveguide array for which an exact and an approximate solution are known. In our analysis, we show that to third order the normalized matrix method approximate solution gives results coinciding with the exact known solution.