clear;
clc;
disp('     L                ___');
disp(' -----------           |') ;
disp(' |          |  A       |');
disp(' ----    ---           |   '); 
disp('    |   |              |  D');
disp('    |   |              |');
disp('    |   |              |');
disp('    |   |              |');
disp('    |___|             _|_');
disp('      B                ');
disp(' ');

fcprima=input('INGRESE "fc prima" en megapascales: ');
ES=2100000;
EC=15100*sqrt(fcprima);
FY=2800;
n=ES/EC;

L=input('INGRESE L (largo del ala): ');
A=input('INGRESE A (espesor del ala): ');
B=input('INGRESE B (espesor del alma): ');
D=input('INGRESE D (alto de la viga): ');
b=input('INGRESE CUANTAS BARRAS DE ACERO: ');
d=input('INGRESE DIAMETRO DE LAS BARRAS EN MM: ');
r=input('INGRESE DISTANCIA DESDE TOPE DEL ALA HASTA EL CENTRO DE LAS BARRAS: ');


as=(3.1415*d^2/400); % en cm^2
ac=L*A+B*(D-A)+(n-1)*as;


cm=[L*A*A/2+B*(D-A)*(A+(D-A)/2)+as*(n-1)*(r)]/[L*A+B*(D-A)+(n-1)*as];
cm

i=[L*A^3/12+L*A*(cm-A/2)^2+B*(D-A)^3/12+B*(D-A)*(A+(D-A)/2-cm)^2+(n-1)*as*(d-cm)^2];
i;