% sript for testing structure estimating function

A=rand(2);
[L1,D1]=ldform(A);
B=rand(2);
[L2,D2]=ldform(A+B);

% structure estimation
belief = 2;           % default belief
nbest = 1;           % nbest: how many regressors are returned
nrep = 5;         % nrep: number of random repetions of structure estimation
lambda   = 0.9;
k=2;

[o1,o2,o3]=straux1(L1,D1,10,L2,D2,20, belief, nbest, nrep, lambda, k);
