xa=1xb=2
ya=1
yb=2
x=(xa-0.5):0.01:(xb+0.5)
y=(ya-0.5):0.01:(yb+0.5)
[X,Y]=meshgrid(x,y)
Z=X.*Y
M=(Y<ya|Y>yb)|(X<Y|X>xb)
Z(M)=0
surf(X,Y,Z)
hold on
title('待求體積的曲頂柱體')
x2=[1,1.5,2]
x3=repmat(x2,3,1)
y2=x3
z2=[0 2 4]'
z3=repmat(z2,1,3)
surf(x3,y2,z3)
material metal
colormap(white)
ssym=vpa(int(int('x*y','x','y',2),'y',1,2))

|