Stuff

clc;   
clear; 
format long;
C=0.55*1.2*0.4;
I=1;            
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)=200000;
        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);
        display(' '); 
    end
    keq(i,1)=r; 
        
end

m=0;
n=0;

for i=1:nx  
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;  

for i=1:nx
    for j=1:nm(i,1)
	fza(j,i)=u(j,i)*g(i,1); 
    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