1、119 项目四 无穷级数与微分方程 实验 1 无穷级数 实验目的 观察无穷级数部分和的变化趋势,进一步理解级数的审敛法以及幂级数部分和对函数的 逼近. 掌握用 Mathematica 求无穷级数的和, 求幂级数的收敛域, 展开函数为幂级数以及展 开周期函数为傅里叶级数的方法. 基本命令 1. 求无穷和的命令 Sum 该命令可用来求无穷和. 例如,输入 Sum1/n2,n,l,Infinity 则输出无穷级数的和为 命令 Sum 与数学中的求和号 相当 6/2 2. 将函数展开为幂级数的命令 Series 该命令的基本格式为 Seriesfx,x,x0,n 它将 展开成关于 的幂级数. 幂级数的
2、最高次幂为 余项用 表)(xf 0x ,)(0nx10)(nx 示. 例如,输入 Seriesyx,x,0,5 则输出带皮亚诺余项的麦克劳林级数 65432 01202106100 xoyxyxyxyy 3. 去掉余项的命令 Normal 在将 展开成幂级数后, 有时为了近似计算或作图, 需要把余项去掉. 只要使用)(xf Normal 命令. 例如 ,输入 SeriesExpx,x,0,6 Normal% 则输出 7 65432 xO!x!x1 4. 强制求值的命令 Evaluate 如果函数是用 Normal 命令定义的, 则当对它进行作图或数值计算时, 可能会出现问题. 例如,输入 fx
3、=NormalSeriesExpx,x,0,3 Plotfx,x,-3,3 则只能输出去掉余项后的展开式 6x213 而得不到函数的图形. 这时要使用强制求值命令 Evaluate, 改成输入 PlotEvaluatefx,x,-3,3 则输出上述函数的图形. 5. 作散点图的命令 ListPlot 120 ListPlot 为平面内作散点图的命令, 其对象是数集,例如,输入 ListPlotTablej2,j,16,PlotStyle-PointSize0,012 则输出坐标为 的散点图.16,3,212 6. 符号“/；”用于定义某种规则,“/；”后面是条件. 例如,输入 Clearg,g
4、f; gx_:=x/;0PointSize0.02; 则输出(2)中级数部分和的的变化趋势图. 例 1.2 画出级数 的部分和分布图. 1)(n 输入命令 Clearsn,g;sn=0;n=1;g=;m=3; 121 While1/n10-m,sn=sn+(-1)(n-1)/n; g=Appendg,GraphicsRGBColorAbsSinn,0,1/n, Linesn,0,sn,1;n+; Showg,PlotRange-0.2,1.3,Axes-True; 则输出所给级数部分和的图形,从图中可观察到它收敛于 0.693 附近的一个数. 例 1.3 设 求 .,!10na1na 输入 C
5、leara; an_=10n/(n!); vals=Tablean,n,1,25; ListPlotvals,PlotStyle-PointSize0.012 则输出 的散点图,从图中可观察 的变化趋势. 输入nana Suman,n,l,Infinity 则输出所求级数的和. 求幂级数的收敛域 例 1.4 求 的收敛域与和函数. 021)3(4nnx 输入 Cleara; an_=4(2n)*(x-3)n/(n+1); stepone=an+1/an/Simplify 则输出 n2)x3(16 再输入 steptwo=Limitstepone,n-Infinity 则输出 )x3(16 这里
6、对 an+1和 an都没有加绝对值 . 因此上式的绝对值小于 1 时, 幂级数收敛; 大于 1 时发散. 为了求出收敛区间的端点, 输入 ydd=Solvesteptwo=1,x zdd=Solvesteptwo=-1,x 则输出 1647x1649x与 由此可知,当 时,级数收敛,当 或 时,级数发散. 71649x 为了判断端点的敛散性, 输入 Simplifyan/.x-(49/16) 则输出右端点处幂级数的一般项为 1n 122 因此,在端点 处,级数发散. 再输入1649x Simplifyan/.x-(47/16) 则输出左端点处幂级数的一般项为 1n)( 因此,在端点 处, 级数
7、收敛. 647x 也可以在收敛域内求得这个级数的和函数. 输入 Sum4(2n)*(x-3)n/(n+1),n,0,Infinity 则输出 )x3(16Log 函数的幂级数展开 例 1.5 求 的 6 阶麦克劳林展开式.cos 输入 SeriesCosx,x,0,6 则输出 7 642xo0x1 注:这是带皮亚诺余项的麦克劳林展开式. 例 1.6 求 在 处的 6 阶泰勒展开式.ln1 输入 SeriesLogx,x,1,6 则输出 .xo6)1(5)x(4)1(3)x(2)1(x 7 例 1.7 求 的 5 阶泰勒展开式.arctn 输入 serl=SeriesArcTanx,x,0,5;
8、 Poly=Normalserl 则输出 的近似多项式xrt 53 通过作图把 和它的近似多项式进行比较. 输入arctn PlotEvaluateArcTanx,Poly,x,-3/2,3/2, PlotStyle-Dashing0.01,GrayLevel0,AspectRatio-l 则输出所作图形, 图中虚线为函数 ,实线为它的近似多项式.xarctn 傅里叶级数 例 1.8 设 是以 为周期的周期函数,它在 的表达式是)(xg2, xxg0,1)( 123 将 展开成傅里叶级数.)(xg 输入 Clearg; gx_:=-1/;-PiRGBColor0,1,0,RGBColor1,0
9、.3,0.5,DisplayFunction-Identity; (*tun是以 n 为参数的作图命令*) tu2=Tabletun,n,1,30,5; (*tu2 是用 Table 命令作出的 6 个图形的集合*) toshow=Partitiontu2,2; (*Partition 是对集合 tu2 作分割, 2 为分割的参数*) ShowGraphicsArraytoshow (*GraphicsArray 是把图形排列的命令*) 则输出 6 个排列着的图形,每两个图形排成一行. 可以看到 越大, 的傅里叶级数的前n)(xg 项和与 越接近.n)(xg 实验习题 1.求下列级数的和: (
10、1) (2) (3) (4);21k ;)12(k ;)2(1k .)1(k 2. 求幂级数 的收敛域与和函数. 01)5(nnx 3. 求函数 的 6 阶麦克劳林多项式.l 4. 求 的 6 阶麦克劳林多项式.xarcsi 5. 设 ,求 的 5 阶和 10 阶麦克劳林多项式,把两个近似多项式和函数的1)(2f)(xf 图形作在一个坐标系内. 124 6. 设 在一个周期内的表达式为 , 将它展开为傅里叶级)(xf 211)(2xxf 数(取 6 项), 并作图. 7. 设 在一个周期内的表达式为 , 将它展开为傅里叶级数)(f ,20)(xf (取 8 项), 并作图. 8. 求级数 的和
11、的近似值. 1sink 实验 2 微分方程 实验目的 理解常微分方程解的概念以及积分曲线和方向场的概念,掌握利用 Mathematica 求微分方程及方程组解的常用命令和方法. 基本命令 1. 求微分方程的解的命令 DSolve 对于可以用积分方法求解的微分方程和微分方程组,可用 Dsolve 命令来求其通解或特解. 例如,求方程 的通解, 输入023y DSolvey x+3y x+2yx=0,yx,x 则输出含有两个任意常数 C1和 C2的通解: 2Ce1xx2 注:在上述命令中,一阶导数符号 是通过键盘上的单引号 输入的 ,二阶导数符号 要 输入两个单引号,而不能输入一个双引号. 又如,
12、求解微分方程的初值问题: ,10,6,034 xxyyy 输入 Dsolveyx+4 yx+3yx=0,y0=6, y0=10,yx,x (*大括号把方程和初始条件放在一起 *) 则输出 x2x3e148(ey 2. 求微分方程的数值解的命令 NDSolve 对于不可以用积分方法求解的微分方程初值问题,可以用 NDSolve 命令来求其特解.例 如 要求方程 5.0,32xyy 的近似解 , 输入)5.10(x NDSolveyx=yx2+x3,y0=0.5,yx,x,0,1.5 (*命令中的x,0,1.5表示相应的区间*) 则输出 y-InterpolatingFunction0.,1.5,
13、 注:因为 NDSolve 命令得到的输出是解 的近似值. 首先在区间0,1.5 内插入一)(xy 系 125 列点 , 计算出在这些点上函数的近似值 , 再通过插值方法得到nx,21 ny,21 在区间上的近似解.)(y 3. 一阶微分方程的方向场 一般地,我们可把一阶微分方程写为 ),(yxf 的形式,其中 是已知函数. 上述微分方程表明:未知函数 在点 处的斜率等于函数),(yxf yx 在点 处的函数值. 因此,可在 平面上的每一点, 作出过该点的以 为斜率f O),(yf 的一条很短的直线(即是未知函数 的切线). 这样得到的一个图形就是微分方程y ),(xf 的方向场. 为了便于观
14、察, 实际上只要在 平面上取适当多的点,作出在这些点的函数的x 切线. 顺着斜率的走向画出符合初始条件的解,就可以得到方程 的近似的积分曲),(yxf 线. 例如,画出 的方向场.0)(,12ydx 输入 (1 则输出方向场的图形,从图中可以观察到, 当初始条件为 时, 这个微分方程的解介2/10y 于 和 1 之间, 且当 x 趋向于 或 时, 分别趋向于 与 1. )(xy 下面求解这个微分方程, 并在同一坐标系中画出方程的解与方向场的图解. 输入 sol=DSolveyx=1-yx2,y0=0,yx,x; g2=Plotsol1,1,2,x,-3,3,PlotStyle-Hue0.1,T
15、hickness0.005; Showg2,g1,Axes-None,Frame-True; 则输出微分方程的解 ,以及解曲线与方向场的图形. 从中可以看到, 微分方程xey21)( 的解与方向场的箭头方向相吻合. 实验内容 用 Dsolve 命令求解微分方程 例 2.1 求微分方程 的通解.2xey 输入 Clearx,y; DSolvey x+2x*yx=x*Exp-x2,yx,x 或 DSolveDyx,x+2x*yx=x*Exp-x2,yx,x 则输出微分方程的通解: 1Cex21y2 其中 C1是任意常数. 例 2.2 求微分方程 在初始条件 下的特解.0xyeyx21 输入 Cle
16、arx,y; DSolvex*y x+yx-Expx=0,y1=2 E,yx,x 126 则输出所求特解: xey 例 2.3 求微分方程 的通解.xey2cos52 输入 DSolvey x-2y x+5yx=Expx*Cos2 x,yx,x/Simplify 则输出所求通解: )x2Sin1c4x(2Cos)2c81(exyx 例 2.4 求解微分方程 , 并作出其积分曲线.xey 输入 g1=TablePlotEx+x3/3+c1+x*c2,x,-5,5, DisplayFunction-Identity,c1,-10,10,5,c2,-5,5,5; Showg1,DisplayFunc
17、tion-$DisplayFunction; 则输出积分曲线的图形. 例 2.5 求微分方程组 在初始条件 下的特解. 02yxdtet 0,10ttyx 输入 Clearx,y,t; DSolvex t+xt+2 yt=Expt, yt -xt- yt=0, x0=1,y0=0,xt,yt,t 则输出所求特解: )tSintCose(21ty,Costxt 例 2.6 求解微分方程 并作出积分曲线.,)2/5xdx 输入 RGBColor1,0,0,DisplayFunction-Identity; g2=PlotVectorField1,-2y/(x+1)+(x+1)(5/2),x,-0.
18、999,1,y,-4,4, Frame-True,ScaleFunction-(1 Showg1,g2,Axes-None,Frame-True,DisplayFunction-$DisplayFunction; 则输出积分曲线的图形. 127 用 NDSolve 命令求微积分方程的近似解 例 2.7 求初值问题: 在区间1.2,4上的近似解并作图.1,0)1()( 2.xyxy 输入 fl=NDSolve(1+x*yx)*yx+(1-x*yx)*yx=0,y1.2=1,y,x,1.2,4 则输出为数值近似解(插值函数)的形式: y-InterpolatingFunction1.2,4., 用
19、 Plot 命令可以把它的图形画出来.不过还需要先使用强制求值命令 Evalu-ate, 输入 PlotEvaluateyx/.fl,x,1.2,4 则输出近似解的图形. 如果要求区间1.2,4内某一点的函数的近似值, 例如 ,只要输入8.1xy y1.8/.fl 则输出所求结果 3.8341 例 2.8 求范德波尔(Van der Pel)方程 5.0,0)1(2 xxyyy 在区间0,20上的近似解. 输入 Clearx,y; NDSolveyx+(yx2-1)*yx+yx=0,y0=0,y0=-0.5,y,x,0,20; PlotEvaluateyx/.%,x,0,20 可以观察到近似解
20、的图形. 例 2.9 求出初值问题 0)(,1)0(cossin22yxx 的数值解, 并作出数值解的图形. 输入 NDSolveyx+Sinx2*yx+yx=Cosx2, y0=1,y0=0,yx,x,0,10 PlotEvaluateyx/.%,x,0,10; 则输出所求微分方程的数值解及数值解的图形 例 2.10 洛伦兹(Lorenz)方程组是由三个一阶微分方程组成的方程组.这三个方程看似简 单, 也没有包含复杂的函数, 但它的解却很有趣和耐人寻味. 试求解洛伦兹方程组 ,0)(,4)0(,12)(5)(6)(zyxtttztyxttt 并画出解曲线的图形. 输入 Cleareq,x,y
21、,z eq=Sequencext=16*yt-16*xt, yt=-xt*zt-yt+45xt,zt=xt*yt-4zt; sol1=NDSolveeq,x0=12,y0=4,z0=0, xt,yt,zt,t,0,16,MaxSteps-10000; g1=ParametricPlot3DEvaluatext,yt,zt/.sol1,t,0,16, PlotPoints-14400,Boxed-False,Axes-None; 128 则输出所求数值解的图形. 从图中可以看出洛伦兹微分方程组具有一个奇异吸引子, 这个吸 引子紧紧地把解的图形“吸”在一起. 有趣的是, 无论把解的曲线画得多长,
22、这些曲线也不 相交. 改变初值为 输入,10)(,)0(,6)(zyx sol2=NDSolveeq,x0=6,y0=-10,z0=10, xt,yt,zt,t,0,24,MaxSteps-10000; g2=ParametricPlot3DEvaluatext,yt,zt/.sol2,t,0,24, PlotPoints-14400,Boxed-False,Axes-None; ShowGraphicsArrayg1,g2; 则输出所求数值解的图形. 从图中可以看出奇异吸引子又出现了, 它把解“吸”在某个区域 内, 使得所有的解好象是有规则地依某种模式缠绕. 实验习题 1. 求下列微分方程的
23、通解: (1) ;0136yy (2) 24 (3) ;2sin5xe (4) .)(93y 2. 求下列微分方程的特解: (1) ;15,0,4 xxy (2) .2sin 3. 求微分方程 在初始条件 下的特解.分别求精确解cos)1( 10xy 和数值解 并作图.)0(x 4. 求微分方程组 的通解. tteyxdt235 5. 求微分方程组 的特解. 4,081ttyxt 6. 求欧拉方程组 的通解.32y 7. 求方程 在区间0,4上的近似解 .5,11 xx 实验 3 抛射体的运动 实验目的 通过微分方程建模和 Mathematica 软件,在项目一实验 5 的基础上,进一步研 究
24、在考虑空气阻力的情况下抛射体的运动. 问题 根据侦察,发现离我军大炮阵地水平距离 10km 的前方有一敌军的坦克群正以每小 时 50km 向我军阵地驶来,现欲发射炮弹摧毁敌军坦克群. 为在最短时间内有效摧毁敌军坦 克,要求每门大炮都能进行精射击,这样问题就可简化为单门大炮对移动坦克的精确射击 问题. 假设炮弹发射速度可控制在 0.2km/s 至 0.5km/s 之间,问应选择怎样的炮弹发射速度和 129 怎样的发射角度可以最有效摧毁敌军坦克. 说明 本节我们研究受到重力和空气阻力约束的抛射体的射程. 用 记抛射体)(,tyx 的位置, 其中 x 轴是运动的水平方向, y 轴是垂直方向. 通过在
25、 的约束下最大化 x, 可以0y 计算出使抛射体的射程最大的发射角. 假设 时抛射体(炮弹) 在原点 (0,0)以与水平线夹角0t 为 初始速度为 发射出去. 它受到的空气阻力为,0v (3.1).,dtyxkvFr 重力为 (3.2).,0(mgg 在推导 和 所满足的微分方程之前, 补充一点说明:虽然我们将位置变量)(txty 仅写作 t 的函数,但实际上位置变量还依赖于几个其它的变量或参数. 特别是,x 和 y),(txy 也依赖于发射角 、阻力系数 k、质量 m 及重力加速度 g 等. 为了推导 x 和 y 的方程 , 按照牛顿定律 并结合重力的公式(3.2)和空气阻力的公,aF 式(
26、3.1), 得到微分方程 (3.3)0)(txt (3.4)gyk 根据前面所述假设知, 满足下列初始条件),(tx , (3.5)0,)0(yx .sin)(,cos00vv 先求解 ,由方程(3.3),令 可将其化为一阶微分方程t .km 易求出其通解 )(tCetv 由 得 ,所以,cos)0(vxv cos0.)(tmkt 从 通过积分得到 x, 即, .cos)(0Devkttk 由 得 所以,0)(x,s0mD (3.6)1(co)(0tmkevktx 类似地,可从方程(3.4)解出 y. 令 方程化为一阶微分方程, 两端除以 m,得.g 再在上述方程两端乘以积分因子 得.tmke
27、,tmktt gev 130 即 ,)(tmktkgevd 两端积分得 .Ckettm 所以 .tmegv 利用初始条件 确定其中的常数 C 后, 积分 v 得到 y,再次利用初始sin)0(y 条 件 确定任意常数后,则得到0)(y (3.7).sin)1(0tmktmkevetkgy 下面我们利用公式(3.6)与(3.7) 来描绘炮弹运行的典型图形. 假定炮弹发射的初速度为 0.25km/s, 发射角为 , 输入5 Cleara,t,x,y,g,m,k xv_,a_,t_:=(m/k)*v*Cosa Pi/180*(1-Exp-(k/m)*t) yv_,a_,t_:=(g*m/k)(m/k
28、)-t-(m/k)*Exp-(k/m)*t)+ (m/k)*v*Sina Pi/180*(1-Exp-(k/m)*t) g=9.8;m=5.0;k=0.01; 炮弹飞行的时间由炮弹落地时的条件 所确定. 输入0y FindRooty350,55,t=0,t,50 则输出炮弹飞行的时间 t-57.4124 当发射角 时, 输入65 x350,55, 57.4124/N 则输出炮弹的最大射程为 10888.5 现在我们可以画出炮弹运行的典型轨迹了. 输入 ParametricPlotx350,55,t,y350,55,t, t,0,57.4124,PlotRange-0,11000,AxesLab
29、el-x,y 则输出图 3-1. 2000 4000 6000 8000 10000 x 2000 4000 6000 8000 10000 y 131 图 3-1 实验报告 在上述假设下,进一步研究下列问题: (1) 选择一个初始速度和发射角,利用 Mathematica 画出炮弹运行的典型轨迹. (2) 假定坦克在大炮前方 10km 处静止不动,炮弹发射的初速度为 0.32km/s,应选择什么 样 的发射角才能击中坦克?画出炮弹运行的几个轨迹图,通过实验数据和图形来说明你的结论 的合理性. (3) 假定坦克在大炮前方 10km 处静止不动,探索降低或调高炮弹发射的初速度的情况 下,应如何选
30、择炮弹的发射角?从上述讨论中总结出最合理有效的发射速度和发射角. (4) 在上题结论的基础上,继续探索,假定坦克在大炮前方 10km 处以每小时 50km 向大 炮 方向前进,此时应如何制定迅速摧毁敌军坦克的方案? 注:在研究过程中,还要包括适当改变阻力系数 k 与炮弹的质量 m 所带来的变化. 实验 4 蹦极跳运动 实验目的 利用 Mathematica 软件,通过微分方程建模,研究蹦极跳运动. 问题 在不考虑空气阻力和考虑空气阻力等多种情况下,研究蹦极跳运动中,蹦极者与蹦 极绳设计之间的各种关系. 说明 蹦极绳相当于一根粗橡皮筋或有弹性的绳子. 当受到张力使之超过其自然长度,绳 子会产生一
31、个线性回复力, 即绳子会产生一个力使它恢复到自然长度, 而这个力的大小与它 被拉伸的长度成正比. 在一次完美的蹦极跳过程中, 蹦极者爬上一座高桥或高的建筑物, 把 绳的一头系在自己身上, 另一头系在一个固定物体如桥栏杆上, 当他跳离桥时, 激动人心的 时刻就到来了. 这里要分析的是蹦极者从跳出那一瞬间起他的运动规律. 首先要建立坐标系. 假设蹦极者的运动轨迹是垂直的, 因此我们只要用一个坐标来确 定他在时刻 t 的位置. 设 y 是垂直坐标轴, 单位为英尺, 正向朝下, 选择 为桥平面, 时0y 间 t 的单位为秒, 蹦极者跳出的瞬间为 则 表示 t 时刻蹦极者的位置. 下面我们要求,0t)(
32、y 出 的表达式.)(y 由牛顿第二定律, 物体的质量乘以加速度等于物体所受的力. 我们假设蹦极者所受的力 只有重力、空气阻力和蹦极绳产生的回复力. 当然, 直到蹦极者降落的距离大于蹦极绳的自 然长度时, 蹦极绳才会产生回复力. 为简单起见, 假设空气阻力的大小与速度成正比, 比例 系数为 1, 蹦极绳回复力的比例系数为 0.4. 这些假设是合理的, 所得到的数学结果与研究所 做的蹦极实验非常吻合. 重力加速度 ./32sftg 现在我们来考虑一次具体的蹦极跳. 假设绳的自然长度为 蹦极者的体重为,20ftL 160lb ,则他的质量为 斯 . 在他到达绳的自然长度(即 前, 蹦5/160m
33、)20y 极者的坠落满足下列初值问题: ,1vmgdty.)0 利用 Mathematica 求解上述问题. 输入 132 g=32; m=5; L=200; v1t_,y1t_=vt,yt/. DSolvevt=-g-vt/m,yt=vt,v0=0,y0=0,v,y,t 则输出 )te5(e160),(e160 5/t/t5/t5/t 蹦极者坠落 L 英尺所用的时间为 t1=t/.FindRooty1t=-L,t,2 4.00609 现在我们需要找到当蹦极绳产生回复力后的运动初始条件. 当 时, 蹦极者的坠落1t 满足方程 )(4.01yLmvgdtv 初始条件为 解初值问题:.)(,)(t
34、tLy v2t_,y2t_=vt,yt/. DSolvevt=-g-vt/m-0.4*(L+yt)/m, yt=vt,vt1=v1t1,yt1=-L,v,y,t 则输出下列结果 )e)i45.187.20( _e)i42.3.0(e )i98.130.74()i65.1083.(e)i6.49.( ,42317 _ei7.16857.(e )i()9.29(e)i1.0.( ti6457.01()i98.6.( t t)i291.9.t5910(.0 ti. t)i5)i0. t52910.()i198.2406134t59(46 4t)i2. 这个解是用复指数函数来表示的. 现在蹦极者的位置
35、由命令 bungeeyt_=IftAll 则输出位置-时间图形(图 4-1) 10 20 30 40 -600 -400 -200 133 图 4-1 从上图可以看出, 蹦极者在大约 13s 内由桥面坠落 770ft, 然后弹回到桥面下 550ft, 上 下 振动几次, 最终降落到桥面下大约 600ft 处. 实验报告 1.在上述问题中( 求出需要多长时间蹦极者才能到达他运动轨迹上的),160,2wL 最低点, 他能下降到桥面下多少英尺? 2.用图描述一个体重为 195lb, 用 200ft 长绳子的蹦极者的坠落. 在绳子对他产生力之前, 他能做多长时间的“自由”降落? 3.假设你有一根 30
36、0ft 长的蹦极索, 在一组坐标轴上画出你所在实验组的全体成员的运 动轨迹草图. 4.一个 55 岁, 体重 185lb 的蹦极者, 用一根 250ft 长的蹦极索. 在降落过程中, 他达到的 最大速度是多少? 当他最终停止运动时, 他被挂在桥面下多少英尺? 5.用不同的空气阻力系数和蹦极索常数做实验, 确定一组合理的参数, 使得在这组参数 下, 一个 160lb 的蹦极者可以回弹到蹦极索的自然长度以上. 6.科罗拉多的皇家乔治桥(它跨越皇家乔治峡谷 )距谷底 1053ft, 一个 175lb 的蹦极者希 望 能正好碰到谷底, 则他应使用多长的绳子? 7.假如上题中的蹦极者体重增加 10lb,
37、 再用同样长的绳子从皇家乔治桥上跳下, 则当他 撞到乔治峡谷谷底时, 他的坠落速度是多少? 参 考 文 献 1 吴赣昌等. 高等数学多媒体学习系统, 海南出版社, 2005 2 吴赣昌等. 线性代数多媒体学习系统, 海南出版社, 2005 3 吴赣昌等. 概率论与数理统计多媒体学习系统, 海南出版社 , 2005 4 A.D.Andrew, G.L.Cain, S.Crum, T.D.Morley. 用 Mathematica 做微积分实验 . 俞正光, 章纪民译. 清华大学出版社, 2003 5 章栋恩,许晓革. 高等数学实验, 高等教育出版社, 2004 6 上海市教育委员会组编. 高等数学. 科学出版社, 1998 7 赵静 等. 工科数学实验. 高等教育出版社, 1999 8 乐经良, 向隆万, 李世栋. 数学实验. 高等教育出版社, 1999 9 李尚志, 陈发来等. 数学实验. 高等教育出版社, 1999 10 梁浩云. Mathematica 软件与数学实验. 华南理工大学出版社. 2001 11 张韵华. 符号计算系统 Mathematica 教程. 科学出版社. 2001 12 邓建松. Mathematica 使用指南. 彭冉冉译. 麦格劳希尔出版社. 2002
