6-48.m

MATLAB科学计算与工程分析源代码源程序5

clear
x=-1:2/50:1;
y=-2:4/50:2;
z=-3;
mesh(x,y,z)
hold on
mesh(x,y,-z)
mesh(x,-y,z)
mesh(x,-y,-z)
mesh(-x,y,z)
mesh(-x,y,-z)
mesh(-x,-y,-z)
mesh(-x,-y,z)
%积分区域如图6-3所示。
%下面确定积分限：
view(0,90)
title('沿x轴侧视')
view(90,0)
title('沿y轴侧视')
%由以上两个图形以及椭球面的性质，可以得到
syms x y z
f=x^2+y^2+z^3;
a1=-sqrt(1-(x^2));
a2=sqrt(1-(x^2));
b1=-3*sqrt(1-x^2-(y/2)^2);
b2=3*sqrt(1-x^2-(y/2)^2);
fdz=int(f,z,b1,b2) ;
fdzdy=int(fdz,y,a1,a2) ;
fdzdydx=int(fdzdy,1,2) ;
simplify(fdzdydx)