clear close all ex1_dat ex1_geo ex1_fem ex1_KCM ex1_stat ex1_eigen ex1_lti Ts = 100e-6 sysd = c2d(sys,Ts);