s = tf('s');

g = (1+1.5*s)/(3*s*(s+1)^3);

bode(g), grid