Mathematica高数4 Mathematica在积分学中的应用 数学软件与数学实验 教学课件.pptVIP

定积分的思想是“分割、近似、求和、取极限”。 若函数f(x)在区间[a，b]上可积，则f(x)在[a，b]上的定积分为 Clear[f]; f[x_]:=x^2+1; a = 0; b = 2; m = 0; g = Plot[f[x], {x, 0, 2}, PlotStyle-RGBColor[1, 0, 0], DisplayFunction-Identity]; For[n = 0, n100, n = n + 4; t1 = {}; t2 = {}; For[i =0, in, i++, x1 = a+(b-a)i/n; x2 = x1 +(b-a)/n; t1 = Append[t1, Graphics[{RGBColor[0, 1, 0], Rectangle[{x1, f[x1]}, {x2, f[x2]}]}]]; t2 = Append[t2, Graphics[{RGBColor[0, 0, 1], Rectangle[{x1, f[x1]}, {x2, 0}]}]] ]; Show[t1, t2, g, PlotLabel-nintervals] ] 注：Graphics函数的Rectangle选项为作填充矩形图。t1为绿色矩形块，t2为蓝色矩形块。 练习： 求解常微分方程 * * Mathematica For Windows Mathematica I n 高等数学 三、 Mathematica在积分学中的应用 1．Integrate[f[x],x] 功能：求函数关于x的不定积分． 2. Integrate[f[x],{x,x0,x1}] 功能:求函数在区间上的定积分． 3. Integrate[f[x,y],{x,x0,x1},{y,y0,y1}] 功能:求函数的二重积分(注: ， 可以是x的函数)． 4. Integrate[f[x,y,z],{x,x0,x1},{y,y0,y1},{z,z0,z1}] 功能:求函数的三重积分． 基本语句： 通过几何与数值相结合的方法演示定积分的概念和定积分的几何意义。 【例】动画演示小矩形面积的和逼近曲边梯形面积 动画演示小矩形面积的和逼近曲边梯形面积 输出部分图形如图 计算 ，其中 求积分举例-1 Clear[f] f[x_,y_]:=x*y Integrate[f[x,y],{x,0,1}, {y,2*x,x^2+1}] %//N 0.0833333 求积分举例-2 Clear[x,y,z] Integrate[x,{x,0,1},{y,0,(1-x)/2},{z,0,1-x-2y}] 计算 ，其中Ω为三个坐标面与平 面 所围成的区域． 求积分举例-3 ， 计算 ，其中 为柱面 和平面 所围成． f[x_,y_]:=z*Sqrt[x^2+y^2] Integrate[f[x,y],{x,0,2},{y,0,Sqrt[2x-x^2]},{z,0,a}] %//N 另: 采用柱坐标，则有 f=z*Sqrt[x^2+y^2]; x=r*Cos[t]; y=r*Sin[t]; Integrate[f*r,{t,0,Pi/2},{r,0,2Cos[t]},{z,0,a}] I=0.888889 a^2 四、 Mathematica在级数中的应用 函数展成幂级数 Series[Exp[x],{x,0,7}] 2 3 4 5 6 7 x x x x x x 8 1 + x + -- + -- + -- + --- + --- + ---- + O[x] 2 6 24 120 720 5040 Normal[%] 2 3 4 5 6 7 x x x x x x 1 + x + -- + -- + -- + --- + --- + ---- 2 6 24 120 720 5040 Sinx的幂极数