dblquads.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% dblquads.m
% 二重积分函数示例
% 确定积分上下限
x=linspace(0,pi,20);
y=linspace(-pi,pi,20);
[xx,yy]=meshgrid(x,y);
zz=sncs(xx,yy);
% 绘制函数图形
mesh(xx,yy,zz);
title('sin(x)cos(y)+1')
% 计算积分及其相对误差
area=dblquad(@sncs,0,pi,-pi,pi)
rel_err=(area-2*pi^2)/(2*pi^2)