Dear Adel, Thank you very much for your kind help. As you point out, wraparound can be the key to this problem. I learned a lot from your codes, and I am trying to convert it to my system.
BTW, could you please explain why you used 'trapz(transmission,Ky)' ? In my mind, 'transmission.append(Sm.transmission(0,1))' is enough for the transmission spectrum. So, what's the meaning of the integral using 'trapz' ? (In fact there is a code error for this instruction in my laptop even though I noticed the indentation) Thanks again for your help! Sincerely, Jiaqi
