Problem with the evaluation of the product of triple integral

Consider a set of positive values , and .

How can I compute q passing the integrand function as a function handle ? With the followind code an error occurs...

% Inputs:
t       = 15;
sigma   = 0.25;
y       = [1.0 1.2]';

h  = @(x1,x2,x3)       x1.*(t - x2).^x3;
g  = @(x1,x2,x3)       exp(-0.5.*((y - h(x1,x2,x3))./sigma).^2 );  
q  = prod(integral3(@(x1,x2,x3) g(x1,x2,x3), 0, 1, 0, t, 0.3, 1));

Answer 1

[Yongjian Feng]

y is an array. You can't use it directly in integral3. Do the product yourself as a for loop.

t       = 15;
sigma   = 0.25;
y       = [1.0 1.2]';

accu = 1;
for i = 1:length(y)
    h  = @(x1,x2,x3)       x1.*(t - x2).^x3;
    g  = @(x1,x2,x3)       exp(-0.5.*((y(i) - h(x1,x2,x3))./sigma).^2 );  
    q  = integral3(@(x1,x2,x3) g(x1,x2,x3), 0, 1, 0, t, 0.3, 1);
    accu = accu * q;
end
accu

[DIMITRIS GEORGIADIS]

Thank you @Yongjian Feng. Is there a way to compute q using prod() instead of for-loop?

[Yongjian Feng]

Use arrayfun

y = [1.0 1.2];
result = prod(arrayfun(@(x) int3(x), y));


function q = int3(y)
    t       = 15;
    sigma   = 0.25;
    h  = @(x1,x2,x3)       x1.*(t - x2).^x3;
    g  = @(x1,x2,x3)       exp(-0.5.*((y - h(x1,x2,x3))./sigma).^2 );  
    q  = integral3(@(x1,x2,x3) g(x1,x2,x3), 0, 1, 0, t, 0.3, 1);
end

[DIMITRIS GEORGIADIS]

Thank you again @Yongjian Feng!