内蒙古农业大学数值计算方法作业

数值计算方法学院: 水利与土木建筑工程学院学号: 姓名: 班级: 任课教师: 贾德彬目录1.非线性方程求根 .31.1 迭代法 .31.2 牛顿法 .31.3 弦截法 .41.4 二分法 .42.插值 .52.1 线性插值 .52.2 二次插值 .52.3 拉格朗日插值 .62.4 分段线性插值 .72.5 分段二次插值 .73.数值积分 .83.1 复化矩形积分法 .83.2 复化梯形积分法 .93.3 复化辛甫生积分法 .93.4 变步长梯形法 .104.线性方程组数值解法 .104.1 约当消去法 .104.2 高斯消去法 .114.3 三角分解法 .124.4 雅可比迭代法 .135.常微分方程数值解法 .145.1 显式欧拉公式 .145.2 欧拉公式的预测校正系统 .145.3 两步欧拉公式 .155.4 改进欧拉公式 .165.5 四阶龙格库塔法 .171.非线性方程求根1.1 迭代法1.1.1 程序代码Private Sub Command1_Click()Dim x1!, x0!, ep!x0 = InputBox(“请输入初始值 x0“)ep = InputBox(“请输入误差限 ep“)f = 0While f = 0x1 = (Exp(x0) - x0) / 2If (x1 - x0) epx = (x1 + x2) / 2f = (Exp(-(x / 4) * (2 - x) - 1f1 = (Exp(-(x1 / 4) * (2 - x1) - 1If f * f1 0 Thenx1 = xElseIf f = 0 ThenPrint xElseIf f * f1 j Thent = t * (xx - x(j) / (x(i) - x(j)End IfNext jSum = Sum + t * y(i)Next iPrint “x=“; xx, “y=“; Sumf = InputBox(“是否继续(0/1)?“)WendEnd Sub2.3.2 例题已知（4.1，2.0248） ， （4.5，2.1213） ，求点 4.3 的值2.3.3 计算结果2.4 分段线性插值2.4.1 程序代码Private Sub Command1_Click()Dim x(), y()n = InputBox(“请输入插值结点数 n:“)ReDim x(n), y(n)For i = 0 To nx(i) = InputBox(“请输入 x(“ + Str(i) + “)=“)y(i) = InputBox(“请输入 y(“ + Str(i) + “)=“)Next if = 0While f = 0c = InputBox(“请输入未知点的自变量值 X:“)i = 1While i n - 1If c = n - 1If x0 x(i) ThenIf x0 - x(i - 1) x(i) - x0 Thenk = i - 1f = 1Elsek = if = 1End IfElsei = i + 1End IfLoopIf f = 0 Thenk = n - 1End Iff1 = (x0 - x(i + 1) * (x0 - x(i) / (x(i - 1) - x(i + 1) * (x(i - 1) - x(i)f2 = (x0 - x(i + 1) * (x0 - x(i - 1) / (x(i) - x(i + 1) * (x(i) - x(i - 1)f3 = (x0 - x(i) * (x0 - x(i - 1) / (x(i + 1) - x(i) * (x(i + 1) - x(i - 1)y0 = f1 * y(i - 1) + f2 * y(i) + f3 * y(i + 1)a = MsgBox(“