lambda = 10e-3 % constante de temps

TF = tf (1, [lambda 1])
% diagramme de bode
figure; bode (TF)
% trouvons la bande passante pour -3 dB -> 100 rad/s = 15.85 Hz

% réponse indicielle
figure; step(TF)

% réponse sinusoidale pour 15,85 Hz 
[u,t] = gensig ('sin', 1/15.85, 0.5, 0.001);
figure; lsim (TF, u, t)
grid on