各位好，本人在用matlab验证寿命分布符合指数型设备的可靠性，根据公式此类数据的可靠度R服从负对数分布（公式3.27 3.28），具体涉及的计算公式如图所示。我编写的matlab程序计算结果始终与标准答案有较大差距，不知道是哪里出了问题，请大家帮忙看看。谢谢！
clear;
clc;
z = 2;
z0 = 2;
tau = 300;
tau0 = 200;
t = 10;
eta = tau/t
eta0 = tau0/t
%根据公式3.27以及3.28计算rL  标准答案0.8749
f2 = @(r) power((eta0+eta),(z0+z)).*power((-log(r)),(z0-1)).*power(r,(eta0+eta-1))./gamma(z0+z)
f = @(r) integral(@(r) f2(r),r,1)-0.9
rL = fzero(f,1)