function X=lauchli(small,n) % function X=lauchli(small,n) X=[ones(1,n) ; small*eye(n,n)];