matthiasm@8: function phase=princarg(phasein)
matthiasm@8: %phase=princarg(phasein) maps phasein into the [-pi:pi] range
matthiasm@8:     phase=mod(phasein+pi,-2*pi)+pi;