数学实验上机实验报告
西南交通大学数值分析上机实验报告
数值分析上机实习报告学号：姓名：专业：联系电话：任课教师：序 (3)一、必做题 (4)1、问题一 (4)1.1 问题重述 (4)1.2 实验方法介绍 (4)1.3 实验结果 (5)2、问题二 (7)2.1 问题重述 (7)2.2 实验原理 (7)雅各比算法：将系数矩阵A分解为：A=L+U+D，则推到的最后迭代公式为： (8)2.3 实验结果 (8)二、选做题 (10)3、问题三 (10)3.1 问题重述 (10)3.2 实验原理 (10)3.3 实验结果 (11)总结 (11)序伴随着计算机技术的飞速发展，所有的学科都走向定量化和准确化，从而产生了一系列的计算性的学科分支，而数值计算方法就是解决计算问题的桥梁和工具。
数值计算方法，是一种研究并解决数学问题的数值近似解方法，是在计算机上使用的解数学问题的方法。
为了提高计算能力，需要结合计算能力与计算效率，因此，用来解决数值计算的软件因为高效率的计算凸显的十分重要。
数值方法是用来解决数值问题的计算公式，而数值方法的有效性需要根据其方法本身的好坏以及数值本身的好坏来综合判断。
数值计算方法计算的结果大多数都是近似值，但是理论的严密性又要求我们不仅要掌握将基本的算法，还要了解必要的误差分析，以验证计算结果的可靠性。
数值计算一般涉及的计算对象是微积分，线性代数，常微分方程中的数学问题，从而对应解决实际中的工程技术问题。
在借助MA TLAB、JA V A、C++ 和VB软件解决数学模型求解过程中，可以极大的提高计算效率。
本实验采用的是MATLAB软件来解决数值计算问题。
MA TLAB是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，其对解决矩阵运算、绘制函数/数据图像等有非常高的效率。
本文采用MATLAB对多项式拟合、雅雅格比法与高斯－赛德尔迭代法求解方程组迭代求解，对Runge-Kutta 4阶算法进行编程，并通过实例求解验证了其可行性，使用不同方法对计算进行比较，得出不同方法的收敛性与迭代次数的多少，比较各种方法的精确度和解的收敛速度。
离散数学上机实验报告
离散数学上机实验报告————————————————————————————————作者：————————————————————————————————日期：ﻩ《离散数学》实验报告姓名:学号:班级：ﻬ实验一连结词逻辑运算一.实验目的实现二元合取、析取、蕴涵和等价表达式的计算。
熟悉连接词逻辑运算规则,利用程序语言实现逻辑这几种逻辑运算。
二.实验内容从键盘输入两个命题变元P和Q的真值，求它们的合取、析取、蕴涵和等价四种运算的真值。
要求对输入内容进行分析，如果不符合0、1条件需要重新输入，程序有良好的输入输出界面。
三．实验环境使用Miｃroｓoft Visual C+＋6.０为编程软件，采用称C/C＋+语言为编程语言实现。
四.实验过程１.算法分析:合取：p，q都为1的时候为1,其他为0析取：p,q都为0的时候为０，其他为1蕴含：ｐ为1，q为0时为0,其他为１等价:p,q同真同假2.程序代码：#iｎclｕdｅ<stｄio.h>inｔmaiｎ（){ﻩinｔ P,Ｑ,a，b，c，d，p,ｑ;printf（" P的值"）;for(P＝0;P<2；P++)ﻩ{ﻩfｏr(Q=0;Q＜２;Q++)ﻩﻩprintf("\ｔ%d＂,P);ﻩ}printf("\n Q的值");ｆor(Ｐ=0;P<2;P++）ﻩ{ﻩﻩfｏr(Q=０;Q<２;Q＋+)ﻩprintf("＼t％ｄ",Q);ﻩ｝printf("\n 非P的值");for（P=0;P<２;P++){ﻩｆor(Ｑ＝0;Q<2;Ｑ＋+）ﻩ{ﻩﻩｉf(Ｐ＝=０）／*判断非P的值*/ﻩﻩﻩp=1;ﻩﻩeｌsｅﻩﻩp＝0；ﻩﻩﻩｐｒintｆ("＼t%d",p)；ﻩ｝ﻩ}ﻩprintf（"\n 非Q的值＂);ﻩfor(Ｐ＝０;P＜2；P++）ﻩ{ﻩfor（Ｑ＝0；Q<2;Q++)ﻩﻩ{ﻩﻩif(Q==1)／*判断非Q的值*/ﻩﻩﻩq＝０;ﻩelseﻩﻩq=1;ﻩpｒｉｎｔf("\ｔ％ｄ",q);}ﻩ}ﻩprinｔf("\n Ｐ与Q的值"）;foｒ(Ｐ=0;P<2;P++)ﻩ{ﻩfｏr(Ｑ＝0;Q<２;Ｑ++)ﻩ｛ﻩﻩﻩif(Q=＝0||P==0）／*判断P与Q的值*/ﻩa=0;ﻩﻩeｌsｅﻩﻩa=1;ﻩｐrintf("＼t%d",ａ);ﻩ｝ﻩ｝ﻩprintｆ("\n Ｐ或Q的值");ﻩfoｒ（P=０;P<２;P+＋)ﻩ{for(Q＝０;Q<2;Q++)ﻩﻩ{ﻩif(Q=＝1||P＝=1)／＊判断P或Ｑ的值*/ ﻩﻩﻩb=1;ﻩelｓeﻩﻩb＝0；ﻩprintf（"\t%ｄ",b);}ﻩ}ﻩprｉntf（"＼nＰ蕴含Ｑ的值＂);ﻩfｏr(P＝０;P<２;P++)ﻩ{ﻩｆor（Q=０;Q＜２;Q++)ﻩﻩ｛ﻩﻩif(P==1&&Q==0)/＊判断P蕴含Q的值*/ﻩﻩc=0；ﻩﻩelseﻩｃ＝1;ｐｒintｆ（"\t%d＂,ｃ);｝}ｐrintｆ(＂\nP等价Q的值");for(P＝0;P<2;P＋+)ﻩ{ﻩﻩｆｏr(Q=0;Q<2;Ｑ++)ﻩ{ﻩﻩｉｆ（P＝＝Q)/＊判断P等价Q的值*/ ﻩd=1;ﻩﻩｅlｓｅﻩﻩｄ=0;ﻩﻩｐｒiｎtf（＂\ｔ%d",d);}ﻩ｝prｉntｆ("\ｎ")；ｒetｕrｎ 0；}3.实验数据及结果分析:实验二关系的复合运算及逆运算一．实验目的熟悉关系的复合运算和逆运算，编程实现关系复合运算和逆运算算法。
数值分析第一次上机练习实验报告
数值分析第一次上机练习实验报告一、实验目的本次实验旨在通过上机练习，加深对数值分析方法的理解，并掌握实际应用中的数值计算方法。
二、实验内容1. 数值计算的基本概念和方法在本次实验中，我们首先回顾了数值计算的基本概念和方法。
数值计算是一种通过计算机进行数值近似的方法，其包括近似解的计算、误差分析和稳定性分析等内容。
2. 方程求解的数值方法接下来，我们学习了方程求解的数值方法。
方程求解是数值分析中非常重要的一部分，其目的是找到方程的实数或复数解。
我们学习了二分法、牛顿法和割线法等常用的数值求解方法，并对它们的原理和步骤进行了理论学习。
3. 插值和拟合插值和拟合是数值分析中常用的数值逼近方法。
在本次实验中，我们学习了插值和拟合的基本原理，并介绍了常见的插值方法，如拉格朗日插值和牛顿插值。
我们还学习了最小二乘拟合方法，如线性拟合和多项式拟合方法。
4. 数值积分和数值微分数值积分和数值微分是数值分析中的两个重要内容。
在本次实验中，我们学习了数值积分和数值微分的基本原理，并介绍了常用的数值积分方法，如梯形法和辛卜生公式。
我们还学习了数值微分的数值方法，如差商法和牛顿插值法。
5. 常微分方程的数值解法常微分方程是物理和工程问题中常见的数学模型，在本次实验中，我们学习了常微分方程的数值解法，包括欧拉法和四阶龙格-库塔法。
我们学习了这些方法的步骤和原理，并通过具体的实例进行了演示。
三、实验结果及分析通过本次实验，我们深入理解了数值分析的基本原理和方法。
我们通过实际操作，掌握了方程求解、插值和拟合、数值积分和数值微分以及常微分方程的数值解法等数值计算方法。
实验结果表明，在使用数值计算方法时，我们要注意误差的控制和结果的稳定性。
根据实验结果，我们可以对计算结果进行误差分析，并选择适当的数值方法和参数来提高计算的精度和稳定性。
此外，在实际应用中，我们还需要根据具体问题的特点和条件选择合适的数值方法和算法。
四、实验总结通过本次实验，我们对数值分析的基本原理和方法有了更加深入的了解。
数学的上机实验报告
实验题目：线性代数求解方程组一、实验目的1. 理解线性代数中方程组的求解方法。
2. 掌握利用计算机求解线性方程组的算法。
3. 熟悉数学软件（如MATLAB、Python等）在数学问题中的应用。
二、实验内容本次实验主要利用数学软件求解线性方程组。
线性方程组是线性代数中的一个基本问题，其求解方法有很多种，如高斯消元法、矩阵求逆法等。
本实验以高斯消元法为例，利用MATLAB软件求解线性方程组。
三、实验步骤1. 编写高斯消元法算法程序。
2. 输入方程组的系数矩阵和常数项。
3. 调用程序求解方程组。
4. 输出解向量。
四、实验代码及分析1. 高斯消元法算法程序```matlabfunction x = gaussElimination(A, b)[n, m] = size(A);assert(n == m, 'The matrix A must be square.');assert(n == length(b), 'The length of b must be equal to the number of rows in A.');% 初始化解向量x = zeros(n, 1);% 高斯消元for i = 1:n-1% 寻找最大元素[~, maxIdx] = max(abs(A(i:n, i)));maxIdx = maxIdx + i - 1;% 交换行A([i, maxIdx], :) = A([maxIdx, i], :);b([i, maxIdx]) = b([maxIdx, i]);% 消元for j = i+1:nfactor = A(j, i) / A(i, i);A(j, i:n) = A(j, i:n) - factor A(i, i:n); b(j) = b(j) - factor b(i);endend% 回代求解for i = n:-1:1x(i) = (b(i) - A(i, i+1:n) x(i+1:n)) / A(i, i); endend```2. 输入方程组的系数矩阵和常数项```matlabA = [2, 1, -1; 1, 2, 1; -1, 1, 2];b = [8; 5; 2];```3. 调用程序求解方程组```matlabx = gaussElimination(A, b);```4. 输出解向量```matlabdisp('解向量为：');disp(x);```五、实验结果与分析实验结果：```解向量为：2-13```实验分析：通过高斯消元法，我们成功求解了给定的线性方程组。
实验报告三
攀枝花学院实验报告实验课程：Visual C#.NET程序设计教程实验项目：上机实验5 实验日期：系：数计学院成绩：一、实验目的1、区分静态类与非静态类，掌握静态字段、静态方法和静态构造函数的定义方法。
2、理解类的继承性与多态性，掌握其应用方法。
3、理解抽象类、接口的概念，掌握抽象类与接口的定义及使用方法。
4、理解分部类和命名空间的概念，掌握分部类和命名空间的使用方法。
二、实验要求1.熟悉Visual 的基本操作方法。
2.认真阅读本章相关内容，尤其是案例。
3.实验前进行程序设计，完成源程序的编写任务。
4.反复操作，直到不需要参考教材、能熟练操作为止。
三、实验步骤1、设计一个Windows应用程序，在该程序中首先构造一个学生基本类，再分别构造小学生、中学生、大学生等派生类，当输入相关数据，单击不同的按钮（小学生、中学生、大学生）将分别创建不同的学生对象，并输入当前的学生总人数、该学生的姓名、学生类型和平均成绩。
具体要求如下：（1）每个学生都有的字段为：姓名、年龄。
（2）小学生的字段还有语文、数学，用来表示这两科的成绩。
（3）中学生在此基础上多了英语成绩。
（4）大学生只有必修课和选修课两项成绩。
（5）学生类具有方法来统计自己的总成绩，并输出。
（6）通过静态成员自动记录学生总人数。
（7）成员初始化能通过构造函数完成。
源程序如下：using System;using ;using ;using ;using ;using ;using Test3_1{public partial class Form1 : Form{public Form1(){InitializeComponent();}public abstract class Student{protected string name;protected int age;protected static int number;public Student(string name, int age){= name;= age;number++;}public string Name { get { return name; } }public virtual string type{get { return"学生"; }}public abstract double total();public abstract double Average();public string getInto(){string result = ("总人数：{0}，姓名：{1}，{2}，{3}岁", number, Name, type, age);if (type == "小学生")result += ("，平均成绩为{0:N2}:\n", total() / 2);else if(type=="中学生")result += ("，平均成绩为{0:N2}:\n", total() / 3);elseresult += ("，总学分为{0:N2}:\n", total());return result;}}public class Pupil : Student{protected double chinese;protected double math;public Pupil(string name, int age, double chinese, double math): base(name, age){= chinese;= math;}public override string type{get{return"小学生";}}public override double total(){return chinese + math;}}public class Middle : Student{protected double chinese;protected double math;protected double english;public Middle(string name, int age, double chinese, double math,double english) : base(name, age){= chinese;= math;= english;}public override string type{get{return"中学生";}}public override double total(){return chinese + math+english;}}public class University : Student{protected double majors;protected double elective;public University(string name, int age, double majors, double elective): base(name, age){= majors ;= elective;}public override string type{get{return"大学生";}}public override double total(){return majors+elective;}}private void btpupil_Click(object sender, EventArgs e){int age = ;double chinese = ;double math = ;Pupil p = new Pupil, age, chinese, math);+= ();}private void btMiddle_Click(object sender, EventArgs e){int age = ;double chinese = ;double math = ;double english = ;Middle p = new Middle, age, chinese, math,english);+= ();}private void btUniversity_Click(object sender, EventArgs e) {int age = ;double chinese = ;double math = ;University u = new University, age, chinese, math);+= ();}}}运行结果如图所示：2、完善上机实验4-3设计的银行帐户管理系统，增加一个VIP账户的管理。
哈工大_数学实验报告
数学实验报告实验一Matlab的使用1.上机实验各种数据输入方法：程序语句：a=[1 2 3;4 5 6 ;7,8,9] 程序语句：linspace(1,10,5) 等等…………计算结果：a = 计算结果：ans =1 2 34 5 6 1.0000 3.2500 5.5000 7.7500 10.00007 8 92.(1) (a)方法：(b) 方法：程序语句：程序语句：a=[-3 5 0 8;1 -8 2 -1;0 -5 9 3;-7 0 -4 5]; a=[-3 5 0 8;1 -8 2 -1;0 -5 9 3;-7 0 -4 5];b=[0;2;-1;6]; b=[0;2;-1;6];inv(a)*b a\b计算结果：计算结果：ans = ans =-0.6386 -0.6386-0.4210 -0.4210-0.3529 -0.35290.0237 0.0237(2) 4个矩阵的生成语句：矩阵a 的生成语句：e=eye(3,3); a=[e r;o s]r=rand(3,2); 验证语句：o=zeros(2,3); a^2s=diag([1,2]);%此为一个任取的2X2 矩阵b=[e r+r*s; o s^2]计算结果相同：ans =1.0000 0 0 1.9003 1.45790 1.0000 0 0.4623 2.67390 0 1.0000 1.2137 2.28630 0 0 1.0000 00 0 0 0 4.00003.生成多项式的语句：poly ([2,-3,1+2i,1-2i,0,-6])计算结果：ans = 1 5 -9 -1 72 -180 0 计算x＝0.8,-x=-1.2 之值的指令与结果：指令：polyval([1,5,-9,-1,72,-180,0],0.8) 结果：ans= -100.2179指令：polyval([1,5,-9,-1,72,-180,0],-1.2) 结果：ans= 293.29004.求a的指令与结果：指令：a=compan([1,0,-6,3,-8])结果：a =0 6 -3 81 0 0 00 1 0 00 0 1 0求a的特征值的指令与结果：roots(p)的指令与结果为：指令：eig(a) 指令：roots([1,0,-6,3,-8])结果：结果：ans = ans =-2.8374 -2.83742.4692 2.46920.1841 + 1.0526i 0.1841 + 1.0526i0.1841 - 1.0526i 0.1841 - 1.0526i结论：利用友元阵函数a=company(p) 和eig(a) 可以与roots(p)有相同的作用，结果相同。
数据结构实验一 约瑟夫环问题实验报告电子版
for(i = 1;i<length;i++){
tmp = (Node *)malloc(sizeof(Node));
tmp->number = num[i];
tmp->pass = pas[i];
pri->next = tmp;
pri = tmp;
pri->next = head;
for(i=0;i<time;i++){ //找到要删除的结点
tmp = tmp->next;
}
printf("%d ",tmp->number);
timeห้องสมุดไป่ตู้= tmp->pass - 1;
deleteFromList(&head,tmp);//删除结点
tmp = tmp->next;//从下一个结点又开始计算
initList(head);
createFromTail(head,num,pas,sizeof(num)/sizeof(num[0]));
p = head;
printf("\n约瑟夫计数前，每个数和他的密码:\n");
for(i = 0;i<sizeof(num)/sizeof(num[0]);i++){
}
}
// 从链表中删除
void deleteFromList(List *head,Node *tmp)
{
Node *tmp1;
Node *tmp2;
tmp1 = *head;
tmp2 = tmp1;
//如果链表剩了一个元素
《数学建模与数学实验》上机实验报告
成都信息工程大学《数学建模与数学实验》上机实验报告专业信息与计算科学班级姓名学号实验日期成绩等级教师评阅日期[问题描述]下表给出了某一海域以码为单位的直角坐标Oxy 上一点（x，y）（水面一点）以英尺为单位的水深z，水深数据是在低潮时测得的，船的吃水深为5英尺，问在矩形区域（75，200）x （-50，150）里那些地方船要避免进入。
[模型]设水面一点的坐标为（x,y,z），用基点和插值函数在矩形区域（75，200）*（-50，150）内做二维插值、三次插值，然后在作出等高线图。
[求解方法]使用matlab求解：M文件：water.mx=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5];y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.584 -33.5];z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9];cx = 75:0.5:200;cy = -50:0.5:150;[cx,cy]=meshgrid(cx,cy);作出曲面图：代码如下：>> water>> cz=griddata(x,y,z,cx,cy,'cubic');>> meshz(cx,cy,cz)>> xlabel('X'),ylabel('Y'),zlabel('Z')>>作出等高线图：代码如下：>> water>> cz=griddata(x,y,z,cx,cy,'cubic');>> figure(2)>> contour(cx,cy,cz,[-5,-5],'r')>> hold on>> plot(x,y,'*')>> xlabel('X'),ylabel('Y')[结果]插值结果等值图：[结果分析及结论]根据等值图可看出：红色区域为危险区域，所以船只要避免进入。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
y= [1.97687,2.17002,2.34158,2.46389,2.71512,3.06045,3.27829,3.51992,3.8215,4.2 435,4.55188,4.88753]; p = polyfit(x,y,2); xx = min(x):.1:max(x); yy = polyval(p,xx);
plot(x,y,'*'); hold on; plot(xx,yy,'r'); legend('Ô-ʼÊý¾Ý','2½×ÄâºÏÇúÏß',2); hold off; r = vpa(poly2sym(p),2); title(char(r))
4、讨论与分析 答：通过本次试验，我了解到数值计算与 MATLAB 之间的关系，也可以使用 MATLAB 来计算函数，求得函数值，并且可以通过数值来反推函数。
1847
1933
2021
2111
2203
2297
49 至 54 列
2393
2491
2591
2693
2797
2903
55 至 60 列
3011
3121
3233
3347
3463
3581
61 至 66 列
3701
3823
3947
4073
4201
4331
67 至 72 列
4463
4597
4733
4871
5011
end end j
2、判断 n^2-79*n+160 所求是否是素数和判断在 1000 之内，由 n^2-79*n+160 算得的素 数有多少
V= 1至6列
1601 7 至 12 列
1163 13 至 18 列
797 19 至 24 列
503 25 至 30 列
281 31 至 36 列
131 37 至 42 列
实验报告
课程名称
数学实验
实验项目
数学实验上机实践
课程学时 上机地点 实验类型
48
逸夫楼 315，数学科学学院实验室
综合类
实验学时
班级
20161123
学号
2016112301
姓名
作业成绩
A
B
C
D
(共五档)
实验 完成情况 （按要求完 成实验在括 号内打√，没 有按要求完 成实验在括 号内打×）
实验一：函数绘图实验（ ） 实验二：微积分实验（ ） 实验三：数值计算实验（ ） 实验四：怎样计算 Pi？（ ） 实验五：素数实验（ ） 实验六：函数迭代实验（ ） 实验七：最佳分数近似值实验（ ） 实验八：分形几何实验（ ） 实验九：混沌实验（ ） 实验十：概率统计实验（ ）
实验三：用 MATLAB 画 y=sinx/x
程序：x=0:0.001:20*pi
y=sin(x)./x;
plot(x,y);
4、讨论与分析 答：通过此次试验，我学会了使用 MATLAB 软件来绘制函数图像，虽然还不是非常熟 练，但是已经明白了具体的使用方法。但还是有不足之处，还没有深刻的理解很多程序。
实验二：微积分实验
1、实验目的 熟悉并了解使用数学软件，进行微积分问题计算的相关数学软件命令，让学生通过
实验理解微积分，解决微积分计算上的问题。 2、实验内容
⑴求函数 y 2x3 6x2 18x 7 的极值。
⑵计算二重积分 xydxdy ,其中 D 是由 y2 x, y x 2 所围成的区域。
实验二： 1、 做出函数覆盖区域， 2、 求出交点为（4,2），（1，-1）
3、 将二重积分改写为∫（-1→2）ydy∫(y^2→y+2)xdx。
程序：ezplot('x=y^2');
hold on ezplot('x=y+2') smys x y z fistans=int(int(x*y,x,y^2,y+2),y,-1,2);
7351 8413 9547
7523 8597 9743
7697 8783 9941
7873 8971
8051 9161
j = 85
程序：clear all
n=0:100; V=power(n,2)+n+41 P=isprime(V') j=0; for i=1:100
if(V(1,i)<=10000) if (P(i,1)==1) j=j+1; else j=j; end
实验五：素数实验
1、实验目的 通过实验理解素数理论，找出素数理论的某些规律。
2、实验内容 本实验是教材《数学实验》中“实验五 素数”的练习 12，参见教材第 43 页，原
题如下：对 n 0,1, ,100 ，计算 n2 n 41，它们能否都给出素数？在 10000 以内的素
数中，由公式 n2 n 41 给出的素数占多少？类似地，对公式 n2 79n 1601以及公式
53 43 至 48 列
47 49 至 50 列
113
1523 1097
743 461 251 113 47 53 131
P=11111 11111 11111 11111 11111 11111 11001 10111 10
j = 31
程序：clear all n=0:50; V=power(n,2)+n+41 P=isprime(V') j=0; for i=1:50
if(V(1,i)<=1000) if (P(i,1)==1) j=j+1; else j=j; end
end end j
1447 1033
691 421 223 97 43 61
1373 971 641 383 197 83
41 71
1301 911 593 347 173 71
41 83
1231 853 547 313 151
4、讨论与分析 答：通过今天的试验，我学会使用 MATLAB 计算微积分，这样可以更加快速、便捷的 计算出微积分，也可以让我更加深刻的理解微积分。
实验三：数值计算实验
1、实验目的 应用数学软件进行数值计算方面的相关实验，熟悉程序设计方法。通过求数值积分、
数据拟合等问题，并结合函数绘图，去理解数学、应用数学。 2、实验内容
61 43 97
3、判断由 n^3+4*n+31 算得的素数有多少？
13 至 18 列
197
223
251
281
313
347
19 至 24 列
383
421
461
503
547
593
25 至 30 列
641
691
743
797
853
911
31 至 36 列
971
1033
1097
1163
1231
1301
37 至 42 列
1373
1447
1523
1601
1681
1763
43 至 48 列
p=p+((1/16)^n)*((4/(8*n+1))-(2/(8*n+4))-(1/(8*n+5))-(1/(8*n+6))); end vpa(p,20) absolute=abs(p-pi)
4、讨论与分析 答：通过本次试验，更加深刻的认识到了π的计算难度，而且也让我认识到了中国古代 数学家计算出 pi 的伟大
32 曹琰惠
E
实验一：函数绘图实验
1、实验目的 利用数学软件绘制数学函数曲线及曲面，通过实验了解函数图形的绘制方法。
2、实验内容 ⑴在同一个图形中，绘制双曲线 x2 y2 1，以及的双曲线 2 条渐近线。 ⑵在同一个图形中，绘制球面 x2 y2 z2 1与锥面 z x2 y2 相交的曲面。 ⑶自选题目：绘制一个或者多个平面图形、空间曲面图形。
6n2 6n 31做同样的判别，你自己能否给出一个或者多个类似的公式？
3、程序设计及运行结果
1、判断 n^2+n+41 所求是否是素数和判断在 10000 之内，由 n^2+n+41 算得的素数有多
少
结果：V =
1至6列
41
43
47
53
61
71
7 至 12 列
83
97
113
131
151
173
5153
73 至 78 列
5297
5443
5591
5741
5893
6047
79 至 84 列
6203
6361
6521
6683
6847
7013
85 至 90 列 7181
91 至 96 列 8231
97 至 100 列 9353
P=111111 111111 111111 111111 111111 111111 111100 110111 101111 110111 111110 111111 111101 111001 011010 101111 0111
D
⑶自选题目：计算一个或者多个微积分习题。 3、程序设计及运行结果 实验一：
程序：close all; clear all;
clcsyms x f=@(x)2*x^3-6*x^2-18*x+7; [xmin]=fminbnd(f,-5,5), [fxmin]=double(subs(f,x,[xmin])) [xmax]=fminbnd(inline(-f(x)),-5,5) [fxmax]=double(subs(f,x,[xmax])) ezplot(f,[-5,5]); hold on plot([xmin,xmax],subs(f,x,[xmin,xmax]),'ro','LineWidth',5) text(xmin,double(subs(f,x,xmin)+0.2),'¼«Ð¡Öµ') text(xmax,double(subs(f,x,xmax)+0.2),'¼«´óÖµ')