Stuff

clc;
clear;
q=input('ingrese Q basal:  ');
h=input('Ingrese altura "h" del muro:   ');
display(' ');
e=input('Ingrese el Espesor (e) del muro:   ');
display(' ');
cg=input('Ingrese el Centro Geométrico de la planta:   ');
display(' ');
ejes=input('Ingrese el Número de Ejes:   ');
display(' ');
var1=1;
for i=1:ejes
    Eje=i;
    fprintf('Eje:%2.0f ',Eje);
    nmuros=input('Ingrese el numero de muros para el eje:  ');
    display(' ');
    muros(i,1)=nmuros;
    d(1,i)=input('ingrese distancia al primer eje, siendo cero para el primero:  ');
    display(' ');
    for j=1:nmuros
        MATRIX(var1,1)=i;
        Muro=j;
        MATRIX(var1,2)=j;
        fprintf('Muro:%2.0f ',Muro);
        murito=input('Ingrese el Largo del muro del eje:  ');
        MATRIX(var1,3)=murito;
    display(' ');
    fprintf('Muro: %2.0f ',Muro);
        f=input('Ingrese el tipo de Material: 1 si es hormigón armado, 2 si es albañilería:  ');
        if f==1
            E=20000;
        else
            E=15000;
        end
  rigidez(1,j)=(12*E*(e*murito^3/12))/h^3;
  MATRIX(var1,4)=d(1,i);
  MATRIX(var1,5)=rigidez(1,j);
  var1=var1+1;
        
end
  
    
Keje=sum(rigidez,2);
v(1,i)=Keje; 
rigidez=zeros(1,nmuros);

        
end
display('  ');
display('Vector de Rigideces POR EJE');
display(v);
cresistencia=v*d'/sum(v,2);
fprintf('Centro Resistente: %1.0f',cresistencia);
display('  ');
i=1;
for i=1:(ejes-1)
 if cresistencia > d(1,i)
 A(ejes+1,i+1)=-abs(cresistencia-d(1,i+1))*v(1,i);
 end
 if cresistencia < d(1,i)
 A(ejes+1,i+1)=abs(cresistencia-d(1,i+1))*v(1,i);
 end
end
 
if cresistencia > cg
B(ejes+1,1)=q*abs(cresistencia-cg);
end
if cresistencia <=cg
B(ejes+1,1)=-q*abs(cresistencia-cg);
end
 
B(1,1)=q;
A(ejes+1,1)=cresistencia*v(1,1);
A(1,ejes+1)=0;
 
for j=1:(ejes-1)
B(j+1,1)=0; 
A(j+1,1)=1;
end
 
for g=1:(ejes)
A(1,g)=v(1,g);
end
 
for h=1:(ejes-1)
A(h+1,ejes+1)=-d(1,h+1);
end
 
for h=2:(ejes)
A(h,h)=-1;
end
C=A\B;  
res=1;
C=C*10^6;
for i=1:ejes
  for k=1:muros(i,1)
     R(res,1)=i;
     R(res,2)=k;
     R(res,3)=C(i,1)*MATRIX(res,5);
     res=res+1;
  end
end    
display('    EJE     Nº MURO     FUERZA');
display(R);