function tracao = tracoes(conec,\
modu,coord,uu)
%nb nro de barras, nv nro de nos
%modu(1:nb) valor de EA
%coord coordenadas
nb=length(conec);
nv=length(coord);
tracao=zeros(nb,1);
for ib=1:nb
ka=conec(ib,1); kb=conec(ib,2);
Xa=coord(ka,:)';
Xb=coord(kb,:)';
ud(1,1)=uu(2*kb-1)-uu(2*ka-1);
ud(2,1)=uu(2*kb)-uu(2*ka);
d=Xb-Xa;
l0=norm(d);
kkt=modu(ib)/(l0^2);
tracao(ib)=kkt*d'*ud;
end

%%%
end