load('conec.txt');
load('coord.txt');
nb=length(conec)
nv=length(coord)
modu=ones(nb,1);
Kglo=rigidez(conec,modu,coord);
KK=Kglo;
iden=1000*eye(2*nv,2*nv);
KK(1:4,:)=iden(1:4,:);
KK(2*nv-3:2*nv,:)=iden(2*nv-3:2*nv,:);
[s l]=eig(KK);
d=diag(l);
[sd id]=sort(d);
xx=coo2x(coord);
for j=1:20
y=s(:,id(j));
plotsimple(xx+3*y,conec,-1,nv/2+1,-3,3);
pause
end