/*******************************
 * File name: p9-13d_uncomp.ch
 ******************************/

#include <control.h>

int main() {
    // system gain
    double k = 9273;

    // open-loop poles
    array double complex op[3] = {complex(0, 0),
                                  complex(-20, 0),
                                  complex(-40, 0)};
    array double num[1] = {1};
    array double den[1] = {1};
    double tf = 1;

    CControl sys1, sys2, *sys3;
    CPlot plot;

    sys1.model("zpk", NULL, op, k);
    sys2.model("tf", num, den);
    sys3 = sys1.feedback(&sys2);
    sys3->grid(1);
    sys3->step(&plot, NULL, NULL, NULL, tf);
    
    return 0;
}