musynfit 311
# Compare with a random perturbation
yinfrandp = infrpert*step
ymurandp = murpert*step
gph5 = ctrlplot([step,yinfnom(1,1),yinfrandp(1,1),...
ymunom(1,1),ymurandp(1,1)]);
gph5 = plot(gph5,{legend=["input step";"Kinf nominal";...
"Kinf rand. pert.";"Kmu nominal";"Kmu rand. pert."]})?