clc;   % Limpia la pantalla actual y la deja en blanco.
clear; % Borra todas las variables utilizadas anteriormente.
C=0.55*1.2*0.3; % Coeficiente Sismico Maximo segun Nch433, para zona sismica 3, suelo tipo III y R=4.
I=1;            % Coeficiente de Importancia Estructura, en este caso categoria C.
P=1;            

a=input('Ingrese el area de la losa en [cm^2]: ');
display(' ');
S=C*I*P*a/(100^2); 
cm=input('Ingrese distancia en el Eje X, desde el eje más a la izquierda hasta el Centroide de la losa en [cm]: ');
display(' ');
h=input('Ingrese la altura entre pisos en [cm]: ');
display(' ');
nx=input('Cuantos ejes hay en el sentido X:  '); 
display(' ');

for i=1:nx
r=0;
nm(i,1)=input(['Ingrese el Numero de Muros del Eje [' num2str(i) ']: ']);

d(i,1)=input(['Ingrese distancia del eje mas a la izquierda al eje [' num2str(i) '] en [cm]: ']); 
    for w=1:nm(i,1)
        im(w,1)=input(['Ingrese el Largo del Muro [' num2str(w) ']: en [cm]: ']);
        im(w,2)=input(['Ingrese el Espesor del Muro [' num2str(w) ']: en [cm]: ']);
        im(w,3)=input(['Ingrese el Material del Muro [' num2str(w) '], [1] si es Hormigon, [2] si es Albañileria: ']);
    
        if im(w,3)==1 
        im(w,3)=191001;  % h2', fc'=160
        end
        if im(w,3)==2
        im(w,3)=15000;
        end
        
        k(w,i)=(im(w,1)^3*im(w,2)*im(w,3)/h^3);
        
        r=k(w,i)+r;
        u(w,i)=k(w,i); % Esta matriz guarda las rigideces de cada muro.
        display(' '); 
    end
    keq(i,1)=r; % Esta matriz guarda las rigideces equivalentes por eje.
          
end
m=0;
n=0;
for i=1:nx  % Calculo del Centro de Rigidez
m=keq(i,1)*d(i,1)+m;
n=keq(i,1)+n;
end

cr=m/n;

for i=1:(nx-1)
    if cr > d(i+1,1)
        t(nx+1,i+1)=abs(cr-d(i+1,1))*keq(i,1);
    end
    if cr <= d(i+1,1)
        t(nx+1,i+1)=-abs(cr-d(i+1,1))*keq(i,1);
    end
end
 
z(1,1)=S;

if cr > cm
    z(nx+1,1)=S*abs(cr-cm);
end
if cr <= cm
    z(nx+1,1)=-S*abs(cr-cm);
end
t(nx+1,1)=cr*keq(1,1);
t(nx+1,nx+1)=0;
t(1,nx+1)=0;
 
for i=1:(nx-1)
    z(i+1,1)=0;
    t(i+1,1)=1;
end
 
for i=1:(nx)
    t(1,i)=keq(i,1);
end
 
for i=1:(nx-1)
    t(i+1,nx+1)=-d(i+1,1);
end
 
for i=2:(nx)
    t(i,i)=-1;
end
 
g=t\z;  % G es la matriz con los Delta de desplazamiento y theta
for i=1:nx
    for j=1:nm(i,1)
	fza(j,i)=u(j,i)*g(i,1); % F=k*delta
    
    
    end
end    
display(' ');
display(' Resultados: ');
display('       EJE Nº      |      MURO Nº      |     FUERZA (TonF)     |');

b=0;    
w=0;
for i=1:nx
    for j=1:nm(i,1)
    w=w+1;    
        b(w,1)=i;
        b(w,2)=j;
        b(w,3)=fza(j,i);
    end
end
b % Muestra los resultados finales de las fuerzas que toma cada muro.