Add unitary tests for cacsd
[scilab.git] / scilab / modules / cacsd / tests / unit_tests / ricc.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2012 - SE - Sylvestre Ledru
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 A = [8,8,9;1,8,6;2,5,9];
8 B = [0,6;7,8;4,0];
9 C = [8,5,5;9,8,1;5,0,7];
10 C=C*C';
11 R=[2,9; 5,7];
12 R=R*R'+eye();
13 B=B*inv(R)*B';
14 X=ricc(A,B,C,'cont');
15 X_ref=[44.84688,2.304207,66.647924;2.304207,27.392655,10.625534;66.647924,10.625534,118.86653];
16 assert_checkalmostequal(X, X_ref, 1e-6);