Based on the density-matrix approach and iterative treatment, a detailed procedure for the calculation of the linear and nonlinear intersubband refractive index changes (RICs) in wurtzite GaN-based coupling quantum wells (CQWs) is given. The simple analytical formulas for electronic eigenstates...