clear;
format long
Kp = 0.554477220576924;
Kd = 6.553720021533881;
dKp = 0.0001;
dKd = 0.001;
% Titik ke-1
    x(1) = Kp - dKp
    y(1) = Kd - dKd
    Go(1) = input('SIAP 1 = ')
    K = length(J);
    z(1) = J(K)
    pause;
% Titik ke-2
    x(2) = Kp
    y(1) = Kd - dKd
    Go(2) = input('SIAP 2 = ')
    K = length(J);
    z(2) = J(K)
    pause;
% Titik ke-3
    x(3) = Kp + dKp
    y(1) = Kd - dKd
    Go(3) = input('SIAP 3 = ')
    K = length(J);
    z(3) = J(K)
    pause;
% Titik ke-4
    x(3) = Kp + dKp
    y(2) = Kd
    Go(4) = input('SIAP 4 = ')
    K = length(J);
    z(4) = J(K)
    pause;
% Titik ke-5
    x(2) = Kp 
    y(2) = Kd 
    Go(5) = input('SIAP 5 =')
    K = length(J);
    z(5) = J(K)
    pause;
% Titik ke-6
    x(1) = Kp - dKp
    y(2) = Kd
    Go(6) = input('SIAP 6 = ')
    K = length(J);
    z(6) = J(K)
    pause;
% Titik ke-7
    x(1) = Kp - dKp
    y(3) = Kd + dKd
    Go(7) = input('SIAP 7 = ')
    K = length(J);
    z(7) = J(K)
    pause;
% Titik ke-8
    x(2) = Kp
    y(3) = Kd + dKd
    Go(8) = input('SIAP 8 = ')
    K = length(J);
    z(8) = J(K)
    pause;
 % Titik ke-9
    x(3) = Kp + dKp
    y(3) = Kd + dKd
    Go(9) = input('SIAP 9 = ')
    K = length(J);
    z(9) = J(K)
    pause;
clf
plot(Go,z), grid on, xlabel('k-th POINTS'), ylabel('J')
