function Ex0512()
% EX0512    使用全局變量繪制二階系統時域響應
global X
X=0:0.1:20;
z1=0.3;
ex0502(z1);
hold on
z1=0.5;
ex0502(z1);
z1=0.707;
ex0502(z1);

function ex0502(zeta)
%子函數，畫二階系統時域曲線
global X
y=1-1/sqrt(1-zeta^2)*exp(-zeta*X).*sin(sqrt(1-zeta^2)*X+acos(zeta));
plot(X,y);