%Detects eigenvectors of a sparse stiffness matrix that are localized at ends. %If end effects propagate along 18 parameters N=30; %number of base pairs n=6*N+6*(N-1); %number of parameters A = constructStiffHom(N); B = A(1:18,1:18); C = A(19:30,1:18); [x,y]=eig(B); y = diag(y); nb=[]; for i=1:size(x,2), nb=[nb, norm(C*x(:,i))]; end; th=0.8 findem1 = find(nb