§6 Mathmatica求定积分以及相关应用问题练习解答.docVIP

§6 Mathematica求定积分以及相关应用问题练习解答 1. 用Mathematica求解下列定积分： (1); (2); (3); (4) ; (5) 2. 计算下列积分的数值积分 (1); (2) 3. 设，求 4. 分别用矩形法、梯形法、抛物线法计算定积分. 5. 求由两条曲线与围成的平面区域的面积. 6. 求半径为的圆的周长. 7. 求星形线, 的全长. 8. 求圆绕轴旋转一周的旋转体（环体）的体积. 练习5.6答案 1. (1) 解 In[1]:= Out[1]= 解 In[2]:= Out[2]=SinIntegrate[(] 解 In[3]:= Out[3]= (4) 解 In[4]:= Out[4]= (5) 解 In[5]:=Integrate[Log[x],{x,a,b}] Out[5]=a-b-aLog[a]+bLog[b] 2. (1) 解 In[1]:=NIntegrate[Sqrt[1+Sin[x]^3],{x,0,1}] Out[1]=1.08268 (2) 解 In[2]:=NIntegrate[Sin[x]/x,{x,0,Pi}] Out[2]=1.85194 3. 解 In[1]:=f[x_]:=If[x0,1/(1+E^x),1/(1+x)] NIntegrate[f[x-1],{x,0,2}] Out[1]=1.31326 4. 解 (1)矩形法 In[2]:=Clear[y,x,s1,n,b,a]; n=20;a=1;b=5; y[x_]:=; s1=(b-a)/n*Sum[y[a+i(b-a)/n],{i,0,n-1}]//N; s2=(b-a)/n*Sum[y[a+i(b-a)/n],{i,1,n}]//N; Print[“s1=”,s1” s2=”,s2] Out[2]=s1=8.72358 s2=9.03513 梯形法 In[3]:=Clear[y,x,a,b,ss3,s3]; y[x_]:=; n=20;a=1;b=5; ss3=Sum[y[a+i*(b-a)/n],{i,1,n-1}]; s3=(y[a]/2+y[b]/2+ss3)*(b-a)/n //N; Print[“s3=”,s3] Out[3]=s3=8.87936 抛物线法 In[4]:=Clear[y,,x,a,b,s3]; y[x_]:=; n=20;a=1;b=5;m=10; ss1=Sum[(1+(-1)^i)*y[a+i*(b-a)/n],{i,1,n-1}];(*ss1=2y2+2y4+…+2yn-2*) ss2=Sum[(1-(-1)^i)*y[a+i*(b-a)/n],{i,1,n-1}];(*ss2=2y1+2y3+…+2yn-1*) s4=N[(y[a]+y[b]+ss1+2ss2)*(b-a)/3/n,20]; Print[“s4=”,s4] Out[4]=s4=8.8791919143863116989 5. 解 首先画出函数的图形，如图1所示 In[1]:=Plot[{x^2,-Sqrt[x],Sqrt[x]},{x,0,2}] 图1 Out[1]=-Graphics- 然后求出两条曲线的交点： In[2]:=Solve[{y-x^2==0,y^2-x==0},{x,y}] Out[2]={{x(0,y(0},{x(1,y(1}} 再以为积分变量求面积： In[3]:=s=Integrate[Sqrt[y]-y^2,{y,0,1}] Out[3]= 6. 解 取圆心在原点半径为的圆的参数方程为 , 取=1,首先画出圆的图形，如图2所示 In[1]:=ParametricPlot[{Cos[t],Sin[t]},{t,0,2Pi},AspectRatio→Automatic] 图2 Out[1]=-Graphics- 再利用定积分计算曲线的弧长 In[2]:=dx=D[r*Cos[t],t] Out[2]=-rSin[t] In[3]:=dy=D[r*Sin[t],t] Out[3]=rCos[t] In[4]:=s=Integrate[Sqrt[dy^2+dx^2],{t,0,2Pi}] Out[4]=2(r 7. 解 取,首先画出星形线的图形，如图3所示 In[1]:=ParametricPlot[{(Cos[t]^3),(Sin[t]^3)},{t,0,2Pi},AspectRatio(Automatic] 图3 Out[1]= -Graphics