Creates the matrix [ ones(1,N); MU*eye(N) ] The value MU defaults to sqrt(eps).