例1:
clear,clc
fill([0,1,1,0],[0,0,1,0],'y');hold on%绘制积分区域
fill([0.55,0.6,0.6,0.55,0.55],[0,0,0.6,0.55,0],'r') %绘制单元条
gtext('y=x');
pause;
gtext('x=1');
pause;
gtext('y=0') %按照矩形区域调用 dblquad 函数,程序如下:
I=dblquad(@(x,y)(x.^2+y.^2).*(y-x<0),0,1,0,1)
或者 I=integral2(@(x,y)(x.^2+y.^2).*(y-x<0),0,1,0,1)
求得 I=0.3333
例2:
解法一(符号解):
clc, clear
syms r z theta
I1=int((r^2+z^2)*r, z, r^2,sqrt(2-r^2