% sript for testing structure estimating function

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

[o1,o2,o3]=straux1(L1,D1,10,L2,D2,20);