I am not quite certain what you want to do.
To express your elliptic filter as second-order-sections, the Signal Processing Toolbox has a built-in function to do just that.
I would just do:
[z,p,k] = ellip(4,0.2,40,[.41 .47]);
[sos,g] = zp2sos(z,p,k);
figure
freqz(sos, 2^14)
The freqz plot is the Bode plot of the filter.
Note that th filtfilt function can use the ‘sos’ and ‘g’ outputs to define your filter to filter your signal.