数 学 实 验
报 告
实验五 迭代(一)——方程求解
学院:数学与信息科学学院
班级:09级数学(4)班
姓名:***
学号:***
实验五 素数
实验名称迭代(一)——方程求解实验目的用迭代法求解方程及线性方程组实验环境Mathematica4.0系统实验的基本理论与方法方程求根
给定实数域上光滑的实值函数f(x)以及初值定义数列
(1)
称为f(x)的一个迭代序列。
给定迭代函数f(x)以及一个初值利用(1)迭代得到数列如果数列收敛于一个,则有
. (2)
即是方程x=f(x)的解。由此启发我们用如下的方法球方程g(x)=0的近似解。
将方程g(x)=0改写为等价的方程x=f(x), (3)
然后选取一初值利用(1)做迭代。迭代数列收敛的极限就是方程g(x)=0的解。
线性方程组的迭代求解
给定一个n元线性方程组
(6)
或写成距阵的形式Ax=b, (7)
其中是n阶方程, 及均为n维列向量。熟知,当距阵A的行列式非零时方程(7)有唯一的解。迭代法是求解这些问题的有效方法之一。用迭代的方法求解线性方程组的思想与上节介绍的方程求解的方法是类似的。假设我们可以将方程组(7)改写成x=Mx+f
其中是n 阶距阵,是n维列向量。任意给定初始向量,由迭代 (9)
确定向量序列如果收敛到向量,则有
即是方程组(7)的解。实验的内容与步骤一、方程求根
1、考察用迭代函数f(x)=2sin(x)求解方程g(x)=2sin(x)-x=0的解的情况。
(1)在同一直角坐标系中,画出y=f(x)及y=x的图像。从图上观察,方程x=2sin(x)有几个解?
在计算机中打开Mathematica4.0系统;
点击鼠标进入工作区后,输入以下语句
按Shift和Enter键运行。
运行结果:
3、素数的判别
(1)对n=2,3,… ,100,观察2n-1被n整除所得的余数;再取其他的整数m(如3,4,5),观察mn-1被n整除的情况:观察当n为素数时的结果。
m=2时
输入程序:
运行结果:
从运行结果可以发现:当n为素数时,2^n-1被n整除所得的余数都是1。
m=3时
输入程序:
运行结果:
m=4时
输入程序:
运行结果:
(2)对n=2,3,..,300,判断哪些Mersenne数Mn=2n-1是素数?
输入程序:
运行结果:
二、素数的分布
1、素数沿数轴的分布
将素数从小到大排序.用表示相邻的素数间的间隔.计算并画图
n=1000时画出(pn,dn)的图:
t:=
Table[{Prime[i],Prime[i+1]-Prime[i]},{i,1,1000}]
ListPlot[t,PlotStyle?RGBColor[1,0,0]]
n=10000时画出(pn,dn)的图:
t:=
Table[{Prime[i],Prime[i+1]-Prime[i]},{i,1,10000}]
ListPlot[t,PlotStyle?RGBColor[1,0,0]]
2、在二维坐标平面上标出点列观察趋于无穷的趋势,并和函数,以及比较.类似观察,以及, ,和
对任意整数计算, , 看看哪个更接近