au.polarizationUncertainty(I, Q, U, sigma):
CASA <2>: au.polarizationUncertainty(1.26, 0.048, 0.015, 0.005) percentage polarization = 3.991% +- 0.397% Out[2]: 3.9912034795810021