clear
close all
format compact

lame_dat
lame_geo2a
lame_fem

lame_KCM
lame_static
lame_eigen

% ------ Dynamic response -----------------------------------------
freqvec = logspace(0,3,100)'; 	w=2*pi*freqvec;
[sys] = lame_lti(mdof,N_motor1,N_motor2,n_1,n_M,freq,Egv,w);

Ts = 100e-6
sysd = c2d(sys,Ts);
Tend = 0.3