close all; clc;

kp = 0.0017198;
ki = 0.04662;
kd = 0;

s = tf('s');

PID = kp + ki/s + kd*s;
G = 1108.7/(1+0.073597*s);

SYS = feedback(series(PID, G), 1);

opt = stepDataOptions;
opt.StepAmplitude = 500;

step(SYS, opt);
t = [0.005:0.01:0.295];
hold on;
plot(t, output03);