'ma only, need division for ar, use LagPolynomial'

def spdmapoly(self, w, twosided=False):

    if (w is None):
        w = np.linspace(0, np.pi, nfreq)
    return ((0.5 / np.pi) * self.mapoly(np.exp((w * 1j))))
