clear;m=1;c=1;k=10;ki=10;kd=10;yset=0.1;kp=[10,-9.09,-10];
for n=1:3
    dydt=@(t,y) [y(2);(-(c+kd)*y(2)-(k+kp(n))*y(1)+kp(n)*yset+y(3))/m;ki*(yset-y(1))];
    [t,y]=ode45(dydt,[0,50],[0;0;0]);
    subplot(1,3,n);plot(t,y(:,1));grid('on');
    title(['kp=',num2str(kp(n))]);xlabel('Time[s]');ylabel('Height[m]'); 
end