clear;
M = input('Maximum Loop = ')
C1 = (1/220);
C2 = (1/2);
for I = 1:M
    k = I
    x(I) = input('Nilai (Kp, Ki,Kd) = ')
    K = length(J);
    y1(I) = C1*RMSEf(K);
    y2(I) = C2*RMSE(K);
    y(I) = J(K)
    pause;
end
clf
plot(x,y,'k-',x,y1,'k:',x,y2,'k:'), grid on, xlabel('K'), ylabel('RMSE,RMSEf and J')
