octave:2> nb=51; nv=22; modu=1e6*ones(nb,1);
octave:3> load('conec.txt');
octave:4> load('coord.txt');
octave:5> xx=coo2x(coord);
octave:6> plotsimple(xx,conec,-1,11,-4,2)
octave:7> kglo=rigidez(conec,modu,coord);
octave:8> ld(42)=-1000;
octave:9> ld=zeros(2*nv,1);
octave:10> ld(42)=-1000;
octave:11> iden=eye(2*nv,2*nv);
octave:12> kk=kglo;
octave:13> kk(1:4,:)=iden(1:4,:);
octave:14> uu=kk\ld;
octave:15> norm(uu,inf)
ans =  0.67937
octave:16> hold on;
octave:17> plotsimple(xx+4*uu,conec,-1,11,-4,2)
octave:18> tracao=tracoes(conec,modu,coord,uu);
octave:19> tracao(2)
ans =  2.6322e-05
octave:20> uu(1:2)
ans =

   5.7109e-12
  -6.3888e-12

octave:21> uu(3:4)
ans =

   2.4180e-12
   1.9933e-11

octave:22> tracao([45 50 51])
ans =

  -552.34
   447.66
   781.13

octave:23> 552.34-781,13*0.707
ans = -228.66
ans =  9.1910
octave:24> 552.34-781.13*0.707
ans =  0.081090
octave:25> 447.66+781.13*0.7071-1000
ans = -0.0029770
octave:26> rea1(1,1)=kglo(1,:)*uu;
octave:27> rea1(2,1)=kglo(2,:)*uu;
octave:28> rea2(1,1)=kglo(3,:)*uu;
octave:29> rea2(2,1)=kglo(4,:)*uu;
octave:30> rea1,rea2
rea1 =

   1.0000e+04
   5.0000e+02

rea2 =

  -1.0000e+04
   5.0000e+02

octave:31> [w,iw]=max(tracao)
w =  9500.0
iw =  4
octave:32> [w,iw]=min(tracao)
w = -9500.0
iw =  1
octave:33> diary off