%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Programa de analisis modal de edificios %
%Calcula auto valores y auto vectores    %
%Para la vibracion libre de la estructura%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clc
%----------------------------------------------------
%Modulo de lectura de datos
fid= fopen ('datos.dat','r');
[npis,K,M,hpis]=lecdat(fid)
fclose (fid);
%----------------------------------------------------

%----------------------------------------------------
%Modulo de analisis modal
[FI,OM]=modos(npis,K,M);
%----------------------------------------------------

%----------------------------------------------------
%Modulo de normalizacion de la matriz modal
[beta]=normal(npis,M,FI);
%----------------------------------------------------

%----------------------------------------------------
%Modulo de Determinacion de aceleraciones de diseño
%y de fuerzas de piso por modos
[F,omega,V,UI]=aceleracion(npis,M,FI,OM,hpis);
%----------------------------------------------------

%----------------------------------------------------
%Modulo de calculo de Fuerzas de piso
[VV,FF]=srss(npis,F,V);
[VVV,FFF]=cqc(npis,F,OM,omega,V);
%----------------------------------------------------

%----------------------------------------------------
%Modulo para graficar los resultados
[qq]=graficos(npis,hpis,UI,VV,VVV);

%----------------------------------------------------