I wrote the script with the help of the previous questions, but the graph is different from the paper.
I want to plot the conversion efficiency graph when the rod size is l=220nm and w=120nm.
The paper I refer to is https://www.researchgate.net/publication/348888528_Bandwidth_Unlimited_Polarization-Maintaining_Metasurfaces.
This is a picture of a script and a plot I wrote.


Below is a picture of the plot in the paper
