2c2 < ## Copyright (C) 2021 Charles Praplan --- > ## Copyright (C) 2021-2022 Charles Praplan 82a83 > clear nA; 87c88 < A = A(1:nA-1); --- > A = A(1:nB-1); 180a182,195 > > ## Test with 3 trailing zeros both in numerator and in denominator > %!test > %! sos = [1, 1, 0, 1, 0.5, 0; 1, 1, 0, 1, 0.5, 0; 1, 1, 0, 1, 0.5, 0]; > %! [Bh, Ah] = sos2tf (sos); > %! assert (Bh, [1, 3, 3, 1] , 10*eps); > %! assert (Ah, [1, 1.5 0.75 0.125], 10*eps); > > ## Test with 3 leading zeros both in numerator and in denominator > %!test > %! sos = [0, 1, 1, 0, 1, 0.5; 0, 1, 1, 0, 1, 0.5;0, 1, 1, 0, 1, 0.5]; > %! [Bh, Ah] = sos2tf (sos); > %! assert (Bh, [1, 3, 3, 1] , 10*eps); > %! assert (Ah, [1, 1.5 0.75 0.125], 10*eps);