despite of the "bode_phase_unwrap" switch bode_phase tends to wrap the phase for high order transferfunction e.g. H(s) : 1/(1+s)^8;
Considering the poles and zeros of the transferfunction sperately wrapping can be avoided for any transferfunction.
The modified bode2.mac file provides a implementation of such an unwrap function.
The bode_gain is also modified and can plot a constant transferfunction (without an unique variable).