Alphai(: , 1) = {@(x,y) Alpha(1,1)+Alpha(2,1)*x+Alpha(3,1)*y};
Alphai(: , 2) = {@(x,y) Alpha(1,2)+Alpha(2,2)*x+Alpha(3,2)*y};
Alphai(: , 3) = {@(x,y) Alpha(1,3)+Alpha(2,3)*x+Alpha(3,3)*y};
for ind = 1:3
    temp0 = integral2(Alphai(: , ind), xmin , xmax , ymin , ymax);
end

First input argument must be a function handle

for ind = 1:3
    temp0 = integral2(Alphai{: , ind}, xmin , xmax , ymin , ymax);
end