%dblquad_exam2.m
%計算二維高斯函數在圓形區域sqrt(x.^2+y.^2)<1上的二重積分
%在圓形區域外填充0的歸一化高斯函數
f=@(x,y)(1/sqrt(pi)*exp(-x.^2)*1/sqrt(pi)*exp(-y.^2)).*(sqrt(x.^2+y.^2)<=1);   
dblquad(f,-2,2,-2,2,1e-6,@quadl)