请教一个跌落测试可靠库的计算问题

不需要一直跌落到出现故障，只要满足20次不出现故障就停止测试或者在20次以内出现故障就停止测试。
做了10个样品的测试，有9个20次后都没有故障，1个在13次跌落时出现故障。请问如何计算样品在单次跌落中的可靠度？

之前看了下面一篇文章，其计算可靠度的一个前提是：得出每一个样品的出现故障时的次数，与上面问题的前提不一致。
实验数据记录
编号 跌落高度 跌落次数t
1 1米 14
2 1米 15
3 1米 13
4 1米 16
5 1米 10
6 1米 9
7 1米 14
8 1米 12
9 1米 15
10 1米 12
数据处理
以此十组数据作为子样的观测值，用解析法算出均值μ与标准离差估计值σ
均值μ=1/n∑_█(i=1@)^n▒t_i=1/10×（14+15+13+16+10+9+14+12+15+12）=13次
标准离差σ=√((∑▒〖t_1^2-1/n（〗∑▒〖t_i）〗^2)/(n-1))=√((1736-1/10×16900)/(10-1))=2.26次
（1）将上述实验数据按大小次序排列起来，得到如下结果
10件式样按大小排列的结果
序号i 1 2 3 4 5 6 7 8 9 10
次数t 9 10 12 12 13 14 14 15 15 16

（2）根据中位秩计算公式计算出其失效概率估计值F（t_i）=(i-0.3)/(n+0.4)
序号i Xi（ti） Yi[F(ti)] x_i^2 y_i^2 x_iy_i
1 9 0.067 81 0.00449 0.603
2 10 0.163 100 0.02657 1.630
3 12 0.260 144 0.0676 3.120
4 12 0.356 144 0.12674 4.272
5 13 0.452 169 0.2043 5.876
6 14 0.548 196 0.30030 7.672
7 14 0.644 196 0.41474 9.016
8 15 0.740 225 0.5476 11.100
9 15 0.836 225 0.6989 12.540
10 16 0.932 156 0.86862 14.912
∑ 130 5 1736 3.25986 70.741
L_xx=∑▒〖x_1^2-1/n（〗∑▒〖x_i）〗^2=1736-1/10×16900=46
L_yy=∑▒〖y_1^2-1/n（〗∑▒〖y_i）〗^2=3.25986-1/10×25=0.76
L_xy=∑_(i=1)^n▒〖x_iy_i〗-1/n(∑_(i=1)^n▒x_i)(∑_(i=1)^n▒y_i)=70.741-1/10×130×5=5.741
b=L_xy/L_xx=5.741/46=0.1248
a=0.5-0.1248×13=-1.1224

（3）则人工方法拟合的直线为
Y=-1.1224+0.1248X
（4）相关系数为
r=L_xy/√(L_xxL_yy)=5.741/√(46×0.76)=0.971
查表得自由度为8的一行中，起码相关值为0.632，可看出r值大大超过起码值。因此，可以认为x与y这两个变量之间存在着线性关系。则拟合的曲线是由意义的，这样可进一步确认这种手机的抗摔能力的概率分布是可以用正态分布来描述。即X~N（13,5.11）
（5）当X取[1,20]之间的数时，算出手机的可靠度
X 可靠度R
1 0.97478
2 0.96045
3 0.94156
4 0.91522
5 0.8837
6 0.8439
7 0.7988
8 0.7532
9 0.7057
10 0.6628
11 0.6303
12 0.6082
13 0.3989
14 0.3918
15 0.3697
16 0.3372
17 0.2943
18 0.2468
19 0.2012
20 0.1561