* Bug #13531 fixed - Fixed sylv help page to include the real Sylester equation.
[scilab.git] / scilab / modules / cacsd / tests / unit_tests / sylv.tst
index 6c7eecf..c89fcd0 100644 (file)
 A = rand(4, 4); C = rand(4, 3); B = rand(3, 3);
 
 // Solve continuous-time Sylvester equation: A*X+X*B = C
-X = sylv(A, B, C, 'c');
+X = sylv(A, B, C, "c");
 
 expected_norm = 1.D-16;
 
 if norm(A*X+X*B-C)-expected_norm > 1.D-15 then pause,end
 
-// Solve discrete-time Sylvester equation: A*X*B+X = C
-X = sylv(A, B, C, 'd');
+// Solve modified discrete-time Sylvester equation: A*X*B+X = C
+X = sylv(A, B, C, "d");
 
 if norm(A*X*B+X-C)-expected_norm > 1.D-15 then pause,end
+
+// Solve real discrete-time Sylvester equation: A*X*B+X = C
+X = -sylv(-A, B, C, "d");
+
+if norm(A*X*B-X-C)-expected_norm > 1.D-15 then pause,end