n = 5;
p = 2;
P = rand(n,n);
A = P*diag( 2*rand(n,1)-1 )/P;
B = P*rand(n,p);
C = rand(p,n)/P;
D = zeros(p,p);
x0=zeros(n,1);
disp( eig(A) );

nguess = n;
sim('acquisition2.slx',1000);
data = simout.data(1:1000,:);
[A1,B1,C1,result] = scident(simout.data,nguess);
D1 = zeros(p,p);

Te=simout.time(2);