function y=Ex0502(zeta)
% M函數文件二階系統時域響應曲線
% EX0502     Step response of quadratic system.
% zeta  阻尼系數
% y     時域響應
%

% copyright 2003-08-01

x=0:0.1:20;
y=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt(1-zeta^2)*x+acos(zeta))
plot(x,y)