close all; clear all; s = tf('s'); m = 1/4000; k = 1; c = 0.001; %c = 0.0001; %c = -0.001; P = 1/(m*s^2+c*s+k); %P = s/(m*s^2+c*s+k); figure(1); bode(P); mod0909_1;