matlab2013中如何求解被积函数为分段函数的二重积分？

∫0 10∫0 20 f(x,y,1)*f(x,y,2)*f(x,y,....10) dxdy，(0,10),(0,20)分别为积分上下限。
被积函数为的i从1到10 f(x,y,i)的求积,
且function a=f(x,y,i)
if(sqrt(x^2+y^2)*i<=2) a=1;
if(sqrt(x^2+y^2)*i>2) a=0;
本人是新手，请高手指点迷津，非常着急，万分感激！

参考代码：
function zd1766370457469573620
dblquad(@f10,0,10,0,20)

function b=f10(x,y)
b = 1;
for i=1:10
    b = b .* f(x,y,i);
end

function a=f(x,y,i)
a = sqrt(x.^2+y.^2)*i <= 2;求得结果：
ans =
    0.0314

%% 第一题 clear clc syms x y i = int(int(exp(-x^2/2)*sin(x^2+y),x,-2,2),y,-1,1) i = -(5^(1/2)*i^(7/2)*pi^(1/2)*sin(1)*((1 - 2*i)^(1/2)*erf((2 + 4*i)^(1/2)) - (1 + 2*i)^(1/2)*erf((2 - 4*i)^(1/2))))/(5*(i/2)^(1/2)) %% 第二题 t = linspace(-pi,pi,100); x = 2*sin(t); y = cos(t); plot(x,y)