优化积分方法在Matlab中的实践



早前在上数学分析课的时候，研究了matlab上的不同积分方法的速度以及误差率。

分别是：

一般积分法

示性函数法

蒙特卡洛法

并在这一过程中发掘出了Matlab新带的Mupad展现的强大画图功能（之前Matlab的画图效果就不吐嘲了）

一并记录如下。

以下是习题12.3部分例题基于matlab2014b的解答编写

(鉴于matlab2014b加入了带有便捷功能的mupad模块，额外试用mupad进行处理.)
(默认数值显示方式 format rat）

求下面曲面所围区域的体积

例1，z=x^{2}+y^{2},abs{x}=1,abs{y}=1,z=0.

(1)绘制积分区域

在命令窗口打开mupad窗口

>>mupad

在mupad命令窗口绘制积分区域

>>
plot(plot::Implicit3d(x^2 + y^2 =z,x=-2..2, y=-2..2, z=0..2),
     plot::Implicit3d(abs(x)=1, x=-2..2, y=-2..2, z=0..2),
     plot::Implicit3d(abs(y)=1, x=-2..2, y=-2..2, z=0..2),
     plot::Implicit3d(z=0, x=-2..2, y=-2..2, z=0..2))