Design Stable IIR Filter in MATLAB
After working through some possible solutions, we figured out it would be easier to design argitrary magnitude and group delay separately in MATLAB:

fdesign.arbmag
fdesign.arbgrpdelay