早前在上数学分析课的时候,研究了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))