How to solve the following four integration in MATLAB

Hello all, I am trying to solve the following expression involving four integrations in MATLAB but note getting it correctly.
This is how I had tried to code:
integrand = @(y, h, z1, z2) (1 / gamma(m_0)) * …
gammainc( ((u_1 ./ …
(a_1 * zeta_1 * y .* g_1_abs_sq .* h .* P_h1_2 – …
u_1 * a_2 * zeta_1 * y .* g_1_abs_sq .* h .*P_h1_2 – …
u_1 * a_1 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_1 * zeta_1 *P_h1_2 * z2 .* g_1_abs_sq .* h – …
u_1 * a_2 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_2 * zeta_1 * P_h1_2 * z2 .* g_1_abs_sq .* h)) / …
(Omega_0 / m_0)),m_0) .* …
(1/gamma(m_0))*((m_0/Omega_0)^m_0)*y^(m_0-1)*exp(-m_0*y/Omega_0).* 1/(2*pi*A_1).*lambda_1*exp(-lambda_1*z1).*lambda_2*exp(-lambda_2*z2);
outer_integral = @(z2) arrayfun(@(z2_val) integral(@(y,h,z1) integrand(y,h,z1,z2_val ), 0, y_max(z2_val)), 0:1000);

% Perform the integration
op = integral(outer_integral, 0, 1000);
I am not getting why I am getting such errors:
Not enough input arguments.
Error in Analytical_2>@(z1,z2)(P_h1_2*(a_1-u_1*a_2))./(u_1*a_1*z1+u_1*a_1*P_h1_2*z2+u_1*a_2*z1+u_1*a_2*P_h1_2*z2) (line 96)
(u_1 * a_1 * z1 + u_1 * a_1 * P_h1_2 * z2 + …
Error in Analytical_2>@(z2_val)integral(@(y,h,z1)integrand(y,h,z1,z2_val),0,y_max(z2_val)) (line 132)
outer_integral = @(z2) arrayfun(@(z2_val) integral(@(y,h,z1) integrand(y,h,z1,z2_val ), 0, y_max(z2_val)), 0:1000);
Error in Analytical_2>@(z2)arrayfun(@(z2_val)integral(@(y,h,z1)integrand(y,h,z1,z2_val),0,y_max(z2_val)),0:1000) (line 132)
outer_integral = @(z2) arrayfun(@(z2_val) integral(@(y,h,z1) integrand(y,h,z1,z2_val ), 0, y_max(z2_val)), 0:1000);
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 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);
Error in integral (line 87)
Q = integralCalc(fun,a,b,opstruct);
Error in Analytical_2 (line 136)
op = integral(outer_integral, 0, 1000);

Please note that I had intentionally not given other part of code which is used to produce the values of elements that are used in this integrand. Any help in this regard will be highly appreciated.Hello all, I am trying to solve the following expression involving four integrations in MATLAB but note getting it correctly.
This is how I had tried to code:
integrand = @(y, h, z1, z2) (1 / gamma(m_0)) * …
gammainc( ((u_1 ./ …
(a_1 * zeta_1 * y .* g_1_abs_sq .* h .* P_h1_2 – …
u_1 * a_2 * zeta_1 * y .* g_1_abs_sq .* h .*P_h1_2 – …
u_1 * a_1 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_1 * zeta_1 *P_h1_2 * z2 .* g_1_abs_sq .* h – …
u_1 * a_2 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_2 * zeta_1 * P_h1_2 * z2 .* g_1_abs_sq .* h)) / …
(Omega_0 / m_0)),m_0) .* …
(1/gamma(m_0))*((m_0/Omega_0)^m_0)*y^(m_0-1)*exp(-m_0*y/Omega_0).* 1/(2*pi*A_1).*lambda_1*exp(-lambda_1*z1).*lambda_2*exp(-lambda_2*z2);
outer_integral = @(z2) arrayfun(@(z2_val) integral(@(y,h,z1) integrand(y,h,z1,z2_val ), 0, y_max(z2_val)), 0:1000);

Please note that I had intentionally not given other part of code which is used to produce the values of elements that are used in this integrand. Any help in this regard will be highly appreciated. Hello all, I am trying to solve the following expression involving four integrations in MATLAB but note getting it correctly.
This is how I had tried to code:
integrand = @(y, h, z1, z2) (1 / gamma(m_0)) * …
gammainc( ((u_1 ./ …
(a_1 * zeta_1 * y .* g_1_abs_sq .* h .* P_h1_2 – …
u_1 * a_2 * zeta_1 * y .* g_1_abs_sq .* h .*P_h1_2 – …
u_1 * a_1 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_1 * zeta_1 *P_h1_2 * z2 .* g_1_abs_sq .* h – …
u_1 * a_2 * zeta_1 * g_1_abs_sq .* h .* z1 – …
u_1 * a_2 * zeta_1 * P_h1_2 * z2 .* g_1_abs_sq .* h)) / …
(Omega_0 / m_0)),m_0) .* …
(1/gamma(m_0))*((m_0/Omega_0)^m_0)*y^(m_0-1)*exp(-m_0*y/Omega_0).* 1/(2*pi*A_1).*lambda_1*exp(-lambda_1*z1).*lambda_2*exp(-lambda_2*z2);
outer_integral = @(z2) arrayfun(@(z2_val) integral(@(y,h,z1) integrand(y,h,z1,z2_val ), 0, y_max(z2_val)), 0:1000);

Cart

Cart