function output = p3_stepResponse(time)

global omega_n zetaa omega_d psii

output = 1 - exp(-zetaa * omega_n * time)./sqrt(1-zetaa^2) .* sin(omega_d * time +psii);