c_6=0; c_5=9; c_4=6; c_3=4; c_2=5; c_1=9; c_0=0; a = [1, -0.1*c_5, 0.1*c_4, -0.c_3]; b = [0.1*c_2, -0.1*c_1, 0.1*c_0]; %impz(b, a, 10) % opd B stepz(b, a, 25) % opd C