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