[fact1, b21, b11, b01, c11, c01, zzeros1, zpoles1] = syredi(1, 4, [1, 2, 0, 0], 0.02, 0.001);
assert_checkequal(fact1, fact1_ref);
assert_checkequal(b21, b21_ref);
-assert_checkequal(b11, b11_ref);
+assert_checkalmostequal(b11, b11_ref);
assert_checkequal(b01, b01_ref);
assert_checkalmostequal(c11, c11_ref);
assert_checkalmostequal(c01, c01_ref, [], %eps);
-assert_checkequal(zzeros1, zzeros1_ref);
+assert_checkalmostequal(zzeros1, zzeros1_ref);
assert_checkalmostequal(zpoles1, zpoles1_ref);
[fact2, b22, b12, b02, c12, c02, zzeros2, zpoles2] = syredi(3, 4, [1, 2, 3, 3.01], 0.1, 0.001);
assert_checkalmostequal(fact2, fact2_ref);
[fact1, b21, b11, b01, c11, c01, zzeros1, zpoles1] = syredi(1, 4, [1, 2, 0, 0], 0.02, 0.001);
assert_checkequal(fact1, fact1_ref);
assert_checkequal(b21, b21_ref);
-assert_checkequal(b11, b11_ref);
+assert_checkalmostequal(b11, b11_ref);
assert_checkequal(b01, b01_ref);
assert_checkalmostequal(c11, c11_ref);
assert_checkalmostequal(c01, c01_ref, [], %eps);
-assert_checkequal(zzeros1, zzeros1_ref);
+assert_checkalmostequal(zzeros1, zzeros1_ref);
assert_checkalmostequal(zpoles1, zpoles1_ref);
[fact2, b22, b12, b02, c12, c02, zzeros2, zpoles2] = syredi(3, 4, [1, 2, 3, 3.01], 0.1, 0.001);