global a b c d h

a = input('a: ');
b = input('b: ');
k = input('k: ');
t = a:0.01:b;
sintk = sin(t.^k);


lens = strlength(split('Байбородов Серафим Максимович'));
a = lens(1);
b = lens(2);
c = lens(3);
d = 3;
h = 8;

w = 0:0.05:2*pi;
z = w*1i;

subplot(2, 2, 2);
plot(w, abs(F(z)));
xlabel('\omega');
ylabel('|F(z)|');
grid on;

subplot(2, 2, 3);
plot(w, unwrap(angle(F(z))));
xlabel('\omega');
ylabel('arg(F(z))');
grid on;

subplot(2, 2, 4);
plot(real(F(z)), imag(F(z)));
xlabel('Re(F(z))');
ylabel('Im(F(z))');
grid on;

subplot(2, 2, 1);
comet(t, sintk);
xlabel('t');
ylabel('sin(t^' + string(k) + ")");
grid on;

function Fz = F(z)
    global a b c d h
    %Fz = exp(z);
    Fz = (d*z + h) ./ (a*z.^2 + b*z + c);
end