function lambda = spectrum(K, B)
    [v,lambda] = eig(full(K), full(B));
%     [v,lambda] = eig(full(K));
    lambda = diag(lambda);
end