global nb nv barras restric massa x0
global extforces totforces tensoes grav
global cvec

barras=load('barras.txt');
nb=length(barras);
nos=load('nos.txt');
nv=length(nos);
xx=zeros(2*nv,1);
restric=zeros(2*nv,1);
extforces=xx;
totforces=xx;
tensoes=zeros(nb,1);
grav=-1;

% construcao de xx, vetor posicao inicial
% construcao de restric, informando restricoes

for i=1:nv
 xx(2*i-1)=nos(i,1);
 xx(2*i)=nos(i,2);
 massa(i)=nos(i,3);
 kode=nos(i,4);
 if (kode >= 10) 
   restric(2*i)=1;
 end
 if (kode == 1 | kode == 11) 
    restric(2*i-1)=1;
 end
end

x0=xx;

% estimativa do tamanho medio das barras
tm=norm(barras(:,3),1)/nb;
cvec = tm;