Using prony to obtain MA coefficients

x is an ARMA(p,q) time series. How can the q MA coefficients be obtained using the prony function.

My tests, using simulated data, show that [b a]=prony(x,q,p) gives the AR coefficients in the vector a,

but the vector b does not return the q MA coeffients.