Thank you for your response. My simulation is for a single wavelength and the source injection is normal to the symmetrical axis of the cylinder as shown in the previous image. I got the sweep results of scattering cross section for 3 different wavelengths as below:

and the analytical solution for extinction efficiency must be similar to below graph:

what should I do to get the same results? Meanwhile, how can I calculate the extinction efficiency, directly?