Help in building the magnitude and group delay characteristic of filters

Heelo! Can anyone help with the functions to build the magnitude and the group dealay characteristic for third order normalized(wp=1rad/s) filters :a) Butterworth pass-band attenuation=3dB b)Chebyshev pass-band attenuation=1dB c) Ellyptic ,pass-band attenuation=1dB,stop-band attenuation=30dB d)Bessel I want to build them in the same graphic.