Complex Least Squares FIR filter design using Levinson' s algorithm
From this term in D1:

exp(-j*om1(231:430) *2)


it looks like you're asking it for a filter with a group delay of 2 samples even though the number of taps is 61.  Try setting this to:

exp(-j*om1(231:430) * 61/2)


i.e. set a desired group delay of about half the filter length.