Reference: Chapter 7 of Scattering of Electromagnetic Waves: Numerical Simulations
ptscext.m simulates the extinction coefficient due to a collection of randomly distributed point scatterers. The point scatterers are divided into clusters.
The figures below are generated with the following parameters:
(a) L=50, wavel=0.02, tai=45, fp=0.008905, phi=0, seed=123456, nc=1, npc=250, lc=L=50
(b) L=50, wavel=0.02, tai=45, fp=0.008905, phi=0, seed=123456, nc=25, npc=10, lc=1 and 0.2