% Define random disturbance
rng(0); % set seed for reproducibility
M = rss(3,1,1); % generate 3x1 random matrix with damping=1 and natural frequency=1