pol_order function

Polarization Order Parameter