一元函数积分学与空间图形的画法.docVIP

项目二 一元函数积分学与空间图形的画法 实验1 一元函数积分 基础实验 实验目的 掌握用Mathematica计算不定积分与定积分的方法. 通过作图和观察, 深入理解 定积分的概念和思想方法. 初步了解定积分的近似计算方法. 理解变上限积分的概念. 提高应用 定积分解决各种问题的能力. 基本命令 1.计算不定积分与定积分的命令Integrate: 求不定积分时, 其基本格式为 Integrate[f[x],x] 如输入 Integrate[x^2+a,x] 则输出 其中a是常数. 注意积分常数C被省略. 求定积分时, 其基本格式为 Integrate[f[x], x,a,b ] 其中a是积分下限, b是积分上限. 如输入 Integrate[Sin[x], x,0,Pi/2 ] 则输出 1 注:Mathematica有很多的命令可以用相应的运算符号来代替. 例如,命令Integrate可用积分号 代替, 命令Sum可以用连加号代替, 命令Product可用连乘号代替. 因此只要调出这些运 算符号, 就可以代替通过键盘输入命令. 调用这些命令,只要打开左上角的File菜单,点击Palettes 中的BasicCalculations, 再点击Calculus就可以得到不定积分号、定积分号、求和号、求偏导数 号等等. 为了行文方便, 下面仍然使用键盘输入命令, 但也可以试用这些数学符号直接计算. 2.数值积分命令NIntegrate 用于求定积分的近似值. 其基本格式为 NIntegrate[f[x], x,a,b ] 如输入 NIntegrate[Sin[x^2], x,0,1 ] 则输出 0.310268 3.循环语句For 循环语句的基本形式是 For[循环变量的起始值, 测试条件, 增量, 运算对象] 运行此命令时, 将多次对后面的对象进行运算, 直到循环变量不满足测试条件时为止. 这里必须 用三个逗号分开这四个部分. 如果运算对象由多个命令组成, 命令之间用分号隔开. 例如, 输入 t 0; For[j 1,j 10,j++,t t+j]; t 则循环变量j从取值1开始, 到10结束. 每次增加1. 执行结果, 输出变量t的最终值1+2+… +10 55. 注:For语句中的 j++ 实际表示j j+1 实验内容 用定义计算定积分 当在上连续时, 因此可将 与 作为的近似值. 为了下面计算的方便, 在例1.1中定义这两个近似值为和n的函 数. 例1.1 计算的近似值. 输入 s1[f_, a_,b_ ,n_]: N[ b-a /n*Sum[f[a+k* b-a /n], k,0,n-1 ]]; s2[f_, a_,b_ ,n_]: N[ b-a /n*Sum[f[a+k* b-a /n], k,1,n ]]; 再输入 Clear[f];f[x_] x^2; js1 Table[ 2^n,s1[f, 0,1 ,2^n],s2[f, 0,1 ,2^n] , n,1,10 ]; TableForm[js1,TableHeadings- None, n, s1, s2 ] 输出为 n s1 s2 2 0.125 0.625 4 0.21875 0.46875 8 0.273438 0.398438 16 0.302734 0.365234 32 0.317871 0.349121 64 0.325562 0.341187 128 0.329437 0.33725 256 0.331383 0.335289 512 0.332357 0.334311 1024 0.332845 0.333822 这是的一系列近似值. 且有 不定积分计算 例1.2 求 输入 Integrate[x^2* 1-x^3 ^5,x] 则输出 例1.3 求 输入 Integrate[x^2*ArcTan[x],x] 则输出 定积分计算 例1.4 求 输入 Integrate[Abs[x-2], x,0,4 ] 则输出 4 例1.5 求 输入 Integrate[Sqrt[4-x^2], x,1,2 ] 则输出 例1.6 求 输入 Integrate[Exp[-x^2], x,0,1 ] 则输出 其中Erf是误差函数, 它不是初等函数. 改为求数值积分, 输入 NIntegrate[Exp[-x^2], x,0,1 ] 则有结果 0.746824. 变上限积分 例1.7 求 输入 D[Integrate[w[x], x,0,Cos[x]^2 ],x] 则输出 -2 Cos[x] Sin[x]w[Cos[x]2] 注意这里使用了复合函数求导公式. 例1.8 画出变上限函数及