数值分析试验报告.docVIP

 《数值分析》课程实验报告 （院）系 数理系 ___________ 专业 ______ _信息与计算科学_________________ 班级 1002班学号_ 201010010232__学生姓名 阳 丹_ 指导教师_ 聂存云 ______ ___________ 实验 一 实验时间：2012.11.15 实验名称：逐次分半加速方法求积分 实验目的 ：逐次分半加速方法求积分 题目：计算，要求计算误差小于10^(-5). 程序源程序： package math; public class Test1 { public static void main(String[] args) { for (int i = 3; i = 10; i++) { double s = countSqrt(1, 9, 8.0 / Math.pow(2, i + 1)) - countSqrt(1, 9, 8.0 / Math.pow(2, i)); if (judgeError(s)) { System.out.println(误差小于10^（-5）); double I = countSqrt(1, 9, 8.0 / Math.pow(2, i + 1)) + 1.0 / 3* s; System.out.println(在允许误差内I= + I); break; } } } //复化梯形公式计算I public static double countSqrt(double lowLimit, double upLimit, double h) { double s = 0; for (double i = lowLimit + h; i = upLimit - h; i = i + h) { s = s + Math.sqrt(i); } double I = h / 2 * (1 + 2 * s + 3); return I; } // 判断误差是否大于10^（-5），如果小于10^（-5）返回true，其他情况返回false public static boolean judgeError(double a) { if (Math.abs(a) Math.pow(10, -5)) { return true; } else return false; } } 运行结果： 误差小于10^（-5） 在允许误差内I=17.333333333325616 实验 二 实验时间：2012.11.17 实验名称：高斯消去法 实验目的 ：用高斯消去法解方程组 题目：设有线性方程组 分别用Gauss消去，LU分解法，追赶法解方程组 源程序： package math; /*Gauss消去求解方程组*/ import java.util.Scanner; public class Test_Gauss { static double a[][]; // 系数矩阵 static double b[]; static double x[]; // x的值 static int n; static int n2; // 记录换行的次数 public static void main(String[] args) { Scanner as = new Scanner(System.in); System.out.print(输入方程组的元数：); n = as.nextInt(); a = new double[n + 1][n + 1]; b = new double[n + 1]; x = new double[n + 1]; System.out.println(输入方程组的系数矩阵a：); for (int i = 1; i = n; i++) for (int j = 1; j = n; j++) a[i][j] = as.nextDouble(); System.out.println(输入方程组矩阵b：); for (int i = 1; i = n; i++) b[i] = as.nextDouble(); solveEquation(); back(); print1(); } // 打印二维矩阵 public static void print2() { System.out.println(增广矩阵为：