x = 0:0.01:100;
s = %i * x;

K = 1.4;
Ti = 3.9;
T = 0;
T0 = 0.85;
n = 4;

y = (1 ./ (Ti * s) + 1) * K;
y = y .* exp(-T * s) ./ ((1 + T0 * s) .^ n);

x = real(y);
y = imag(y);
plot2d(x, y);