数学实验报告_2.docVIP

数学实验报告 学生姓名： 陆莹 学 号： 20121315021 院 部： 数学与统计学院 专 业： 统计学 班 级： （1）班 任课教师： 费文龙 实验报告1 实验课程 大学数学实验 实验名称 Mathematica入门 实验日期 2013-9-11 指导老师 费文龙 专业 统计学 年级 2012 姓名 陆莹 学号 20121315021 得分 实验目的： 熟悉Mathematica软件包的使用。 实验内容： 用两种方式编写如下自定义函数，并求其导数f’(x)在x=-2.0,x=1.0,x=5.0处的值 分别用Plot3D, ParametricPlot3D函数画出（）的图像。 用Mathematica实现一个四人追逐问题，给出结果并划出追逐路线（如下图）。 实验要求： 1、撰写实验报告； 2、写出试验过程中所使用的Mathematica程序或语句和计算结果； 实验步骤： 1.Int[1]:=f[x_]:=E^x Sin[x]/;x?0; f[x_]:=Cos[x]/;0x?E; f[x_]:=Sin[x]*Cos[x]/;xE; g[x_]=D[f[x],x]; N[{g[-2.0],g[1.0],g[5.0]}] Int[2]:=Plot3D[1-x^2-y^2,{x,0,1},{y,0,1}] Int[3]:=x[s_,t_]:=Sin[s] Cos[t]; y[s_,t_]:=Sin[s] Sin[t]; z[s_,t_]:=Cos[s]; ParametricPlot3D[{x[s,t],y[s,t],z[s,t]},{s,0,Pi},{t,0,Pi/2}] 3. Int[4]:= v=1;t=18;dt=0.02;n=t/dt; T={{{0,10}},{{10,10}},{{10,0}},{{0,0}}}; d=Sqrt[(x2-x1)^2+(y2-y1)^2]; For[j=1,j?n,j++, For[i=1,i?4,i++,x1=T[[i,j,1]];y1=T[[i,j,2]] If[i?4,x2=T[[i+1,j,1]]; y2=T[[i+1,j,2]], x2=T[[1,j,1]];y2=T[[1,j,2 ] ] ]; x1=x1+v*dt*(x2-x1)/d;y1=y1+v*dt*(y2-y1)/d; T[[i]]=Append[T[[i]],{x1,y1}]]]; P=Graphics[{Line[T[[1]]],Line[T[[2]]],Line[T[[3]]],Line[T[[4]]],Line[{{0,10},{10,10},{10,0},{0,0},{0,10}}]}]; Show[P,AspectRatio?1] 实验结果： Out[1]={-0.179379,-0.841471,-0.839072} Out[2]= Out[3]= Out[4]= 实验报告2 实验课程 大学数学实验 实验名称 pi的近似求解 实验日期 2013-10-16 指导老师 费文龙 专业 统计学 年级 2012 姓名 陆莹 学号 20121315021 得分 实验目的： 练习的求解方法。 实验内容： 1、用反正切函数的幂级数展开式结合有关公式求，若要精确到以40位、50位数字，试比较简单公式和Machin公式所用的项数。 2、用数值积分计算，分别给出用梯形法和Simpson法精确到10位数字、用Simpson法精确到15位数字时所用的项数n及的近似值 3、用计算机模拟Buffon实验，给出n=1,000、10,000、1,000,000时的模拟结果。 实验要求： 撰写实验报告 写出试验过程中所使用的Mathematica程序或语句和计算结果 实验步骤： 1、 m=N[Pi,40] n=N[Pi,50] T[n_,x_]:=Sum[(-1)^(k-1)*x^(2 k-1)/(2 k-1),{k,1,n}]; N