2c2 < ## Copyright (C) 2021 Charles Praplan --- > ## Copyright (C) 2021-2022 Charles Praplan 64c64 < ## 0 1.0000 1.0000 0 1.0000 0.9791 --- > ## 1.0000 1.0000 0 1.0000 0.9791 0 78c78 < DoNotCombineReal = 1; --- > DoNotCombineReal = 0; 168c168 < SOS=[0, 0, 1, 0, 0, 1]; --- > SOS=[0, 0, 1, 0, 0, 1]; # leading zeros will be removed 170a171,183 > # Removing leading zeros if present in numerator and denominator > for count = 1:size(SOS,1) > B=SOS(count, 1:3); > A=SOS(count, 4:6); > while B(1)==0 && A(1)==0 > A(1) = []; # faster than A = circshift(A,-1); > A(3) = 0; # " > B(1) = []; > B(3) = 0; > endwhile > SOS(count,:) = [B,A]; > endfor > 188c201 < %! assert (sos, [0, 0, 1, 0, 0, 1], 100*eps); --- > %! assert (sos, [1, 0, 0, 1, 0, 0], 100*eps); 192c205 < %! assert (sos, [0, 0, 1, 0, 0, 1], 100*eps); --- > %! assert (sos, [1, 0, 0, 1, 0, 0], 100*eps); 196c209 < %! assert (sos, [0, 0, 2, 0, 0, 1], 100*eps); --- > %! assert (sos, [2, 0, 0, 1, 0, 0], 100*eps); 200c213 < %! assert (sos, [0, 0, 1, 0, 0, 1], 100*eps); --- > %! assert (sos, [1, 0, 0, 1, 0, 0], 100*eps); 205c218 < %! assert (sos, [0, 0, 1, 0, 1, 0], 100*eps); --- > %! assert (sos, [0, 1, 0, 1, 0, 0], 100*eps); 209c222 < %! assert (sos, [0, 1, 0, 0, 0, 1], 100*eps); --- > %! assert (sos, [1, 0, 0, 0, 1, 0], 100*eps);