课程要所得单位阶跃响应曲线中的时域性能指标，上升时间、峰值时间、最大超调量、调整时间（2%）这些，自己按照书上写了不知道怎么改
代码如下：
t=0:0.01:2.0;
yss=1;dta=0.02;
%
nG=[100];
tao=sqrt(8)/10;dG=[1 20*tao 100];G1=tf(nG,dG);
[y1,T]=step(G1,t);
%
r=1;while y1(r)<yss;r=r+1;end
trl=(r-1)*0.001;
%
[ymax,tp]=max(y1);tp1=(tp-1)*0.001;
%
mpl=(ymax-yss)/yss;
%
s=1001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;end
tsl=(s-1)*0.001;
运行结果及报错：
索引超出数组元素的数目(201)。
Untitled
索引超出数组元素的数目(201)。
出错 Untitled (第 15 行)
s=1001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;end
求指示我应该怎么改
