Python code for simulating highly focused laser beams using the Richards–Wolf diffraction integral. It computes the 3D vectorial electromagnetic field in the focal region and derives the corresponding 3D Stokes parameters for full polarization analysis.