这个式子用matlab怎么写啊

alpha2=0.1;alpha3=0.5;
r3=@(t3)alpha3.p.(n.*T-t3).*gampdf(t3,Lf.*gama,eta);
s1=integral(r1,(n-1).*T,n.*T); % 第n次检测设备状态为X(t)>Lf的概率；
r4=@(t3,t4)alpha2.p.(t3-t4).*gampdf(t4,Lp.*gama,eta).*gampdf(t3,Lf.*gama,eta);
y=@(t3)t3.*gampdf(t3,Lf.*gama,eta);
Q1=integral2(r4,(n-1).*T,y,(n-1).*T,n.*T);
这个二重积分的式子使用matlab怎样写，总是报错xmax必须为浮点标量