lamda_b1 = 1;
lamda_r1 = 0.5;
lamda_r2 = 0.5;
lamda_br = 0.5;
lamda_rr = 0.3;
a_1 = 0.05;
a_2 = 0.95;
k_1 = (0.1).^2;
k_2 = (0.08).^2;
row_s1 = 5;
row_s = (10./row_s1)^10;
row_r = row_s./2;
a = @(x) (x./(lamda_r2.*row_r));
b = @(x) (x./(a_2-(a_1.*x)));
c = (1./(lamda_br.*row_s));
d = (1./(lamda_b1.*row_s));
f = @(x) exp(-a(x)-(b(x).*(c(x)+d(x))));
g = @(x) 1./(1+((k_2.*lamda_rr.*row_r.*x)./((a_2-(a_1.*x)).*lamda_br.*row_s)));
h = @(x) 1./(1+((k_1.*lamda_r1.*row_r.*x)./((a_2-(a_1.*x)).*lamda_b1.*row_s)));
i = @(x) (1./(1+x));
fun = @(x) f(x).*g(x).*h(x).*i(x);
q = integral(fun,0,Inf)
error----------------------------------
Subscript indices must either be real positive integers or logicals.
Error in @(x)exp(-a(x)-(b(x).*(c(x)+d(x))))
Error in @(x)f(x).*g(x).*h(x).*i(x)
Error in integralCalc/iterateScalarValued (line 314) fx = FUN(t);
Error in integralCalc/vadapt (line 132) [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);
Error in integralCalc (line 83) [q,errbnd] = vadapt(@AToInfInvTransform,interval);
Error in integral (line 88) Q = integralCalc(fun,a,b,opstruct);
Error in test_1 (line 27) q = integral(fun,0,Inf)
