大学数学实验报告.pptVIP

从理论上来讲筛法和试除法可以求所有的素数，但对于很大的素数是难以实现的， (七)、实验总结： (八)、指导教师评语栏： 数学实验报告 专业:数学与应用数学 班级:应数1班 姓名:陈苏梁 学号:091153014 数学实验报告一 (一)、实验名称：怎样计算π (二)、实验目的：确定几种计算π的方法 (三)、实验内容：寻找方法采用mathematica输入程序代码得出π的值 (四)、实验日期：2011.11.22用时2小时 (五)、实验原理：用数值积分计算π；用泰勒级数计算π；蒙特卡罗计算π (六)、实验过程： 用mathematica4.0软件输入以下程序： n=5000;y[x_]:=4/(1+x*x); s1=(Sum[y[k/n],{k,1,n-1}]+(y[0]+y[1])/2)/n; s2=(y[0]+y[1]+2*Sum[y[k/n],{k,1,n-1}]+4*Sum[y[(k-1/2) /n],{k,1,n}])/(6*n); Print[{N[s1,20],N[s2,30],N[pi,30]}] 实验结果： 用mathematica4.0软件输入以下程序： 的v T[x_,n_]:=Sum[(-1)^k*x^(2k+1)/(2k+1),{k,0,n}]; N[4*T[1,20000],20]//Timing T[x_,n_]:=Sum[(-1)^k*x^(2k+1)/(2k+1),{k,0,n}]; Print[N[4*(T[1/2,260]+T[1/3,170]),150]]; Print[N[16*T[1/5,110]-4*T[1/239,30],150]]; Print[N[Pi,150]] 实验结果： 用mathematica4.0软件输入以下程序： n = 10000; p = {}; Do[m = 0; Do[x = Random[]; y = Random[]; If[x^2 + y^2 = 1, m++], {k, 1, n}]; AppendTo[p, N[4m/n]], {t, 1, 10}]; Print[p]; Sum[p[[t]], {t, 1, 10}]/10 实验结果： (七)、实验总结： (八)、指导教师评语栏： 从实验结果可以看出，数值积分法和泰勒 级数法结果相近且准确。从实验结果可以 粗略地得出π=3.1415926…… 数学实验报告二 (一)、实验名称：迭代--分形 (二)、实验目的：以迭代的观点介绍分形的基本特征以及生成图形的基本方法 (三)、实验内容：采用mathematica作出Koch曲线，Minkowski“香肠”，树木花草，Weierstrass函数 Mandelbrot集以及它的局部放大，Julia集以及它的局部放大 (四)、实验日期：2011.11.23用时4小时 (五)、实验原理：部分与整体以某种方式相似的形体称为分形 (六)、实验过程： （1）Koch曲线 用mathematica4.0软件输入以下程序： redokoch[ptlist_List]:= Block[{tmp={},i,pnum=Length[ptlist]}, For[i=1,ipnum,i=i+1,tmp=Join[tmp,{ptlist[[i]], ptlist[[i]]*2/3+ptlist[[i+1]]/3, (ptlist[[i]]+ptlist[[i+1]])/ 2+{ptlist[[i]][[2]]-ptlist[[i+1]][[2]], ptlist[[i+1]][[1]]-ptlist[[i]][[1]]}*Sqrt[3]/6, ptlist[[i]]/3+ptlist[[i+1]]*2/3, ptlist[[i+1]]}]];tmp] lnko01={{0,0},{1,0}}; Show[Graphics[Line[Nest[redokoch,lnko01,5]], AspectRatio-Sqrt[3]/6]] 实验结果： （2）Minkowski 香肠 用mathematica4.0软件输入以下程序： redominkowski[ptlist_List]:= Block[{tmp={},tmp1, i, pnum=Length[ptlist]}, For[i =1 , i pnum, i = i +1, tmp1={ptlist[[i]][[2]]-ptlist[[i+1]][[2]], ptlist[[i+1]][[1]]-ptlist[[i]][[1]]}