National Instruments NI MATRIXx Xmath - Manual (Page 109)

Chapter 6 Tutorial
Xmath Model Reduction Module 6-14 ni.com
Generate Figure 6-10:
syscl = feedback(sysol);
sysolr=sys*syscr;
sysclr=feedback(sysolr);
compare(syscl,sysclr,w,{radians,type=5})
f10=plot({keep,legend=["original","reduced"]})
Figure 6-10. Closed-Loop Gain with wtbalance