function Kglo = rigidez(conec,\
modu,coord)
%nb nro de barras, nv nro de nos
%modu(1:nb) valor de EA
%coord coordenadas
nb=length(conec);
nv=length(coord);
Kglo=zeros(2*nv,2*nv);
for ib=1:nb
Kloc=zeros(4,4);
na=conec(ib,1); nb=conec(ib,2);
Xa=coord(na,:)';
Xb=coord(nb,:)';
d=Xb-Xa;
l0=norm(d);
kk=modu(ib)/(l0^3);
aux=d*d';
Kloc(1:2,1:2)=kk*aux;
Kloc(1:2,3:4)=-kk*aux;
Kloc(3:4,1:2)=-kk*aux;
Kloc(3:4,3:4)=kk*aux;
%Kloc agora contem a matriz de rigidez 
% da barra
loc2glo=[2*na-1,2*na,2*nb-1,2*nb];
for j=1:4
for k=1:4
jglo=loc2glo(j);
kglo=loc2glo(k);
Kglo(jglo,kglo)=Kglo(jglo,kglo)+Kloc(j,k);
end
end
end

%%%
end