#include
#define M1 2
#define M2 2
#define M 3 // M = M1+M2-1
#define N1 4
#define N2 3
#define N 6 // N = N1+N2-1
int main() {
double num1[M1] = {1, 2};
double num2[M2] = {1, 2};
double den1[N1] = {1, 20, 0, 0}; //s^2(s+20) = s^3+20*s^2
double den2[N2] = {1, 6, 25};
double num[M];
double den[N];
class CPlot plotbode, plotnyquist;
class CControl sys;
conv(num, num1, num2);
conv(den, den1, den2);
sys.model("tf", num, den);
sys.bode(&plotbode, NULL, NULL, NULL);
sys.nyquist(&plotnyquist, NULL, NULL, NULL);
return 0;
}