function PreobC_D;
w=tf([0 0 1],[50 15 1]);
Ts=4;
tf=c2d(w,Ts);