L = 0.03;
R = 2;
Km = 2;
J_th = 15;
B_th = 0.5;
J_phi = 8;
B_phi = 0.5;
T = 0.05;

num1 = Km;
den1 = [(J_th*L) (J_th*R+B_th*L) (B_th*R+Km^2) 0];

tf1 = tf(num1,den1);
dig1 = c2d(tf1,T);