一、写在前面
实验目的
(1) 熟悉求解线性方程组的有关理论和方法;
(2) 能编程实现雅可比及高斯-塞德尔迭代法、列主元高斯消去法、约当消去,追赶法
(3) 通过测试,进一步了解各种方法的优缺点
(4) 根据不同类型的方程组,选择合适的数值方法
实验内容
**本次实验参考公式**
二、实验过程
【参考代码】
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
//求解线性方程组
#define N_EQUATIONS 3 //方程组个数,即未知元个数
#define N_DETERMINANT 4 //行列式阶数
double* Gausss_Seidel(double A[][N_EQUATIONS], double B[], double e, int N)
// 输入参数: A 系数矩阵, B 初始向量, e 误差限, N 最大迭代次数
{
int i, j;
int k = 0; //记录当前迭代次数
double ACC; //累加器
double MAX = 0.0;
double Y[N_EQUATIONS];
double X[N_EQUATIONS];
for<