最近在上计算方法这门课，要求是用MATLAB做练习题，但是我觉得C语言也很棒棒啊~

题目：

和直接法不同，迭代法是一种逐次逼近的方法，将复杂问题简单化，求比较大的方程组时一般都不会用直接法。迭代法有好几种，这里使用了Jacobi迭代与Gausse_Seidel迭代法。

使用VS2017，代码如下：

//使用Jacobi迭代法与Gausse_Seidel迭代法计算线性方程组

#include "stdafx.h"

#include#include "math.h"

//根据用户输入的行列数，生成二维矩阵A L U D 向量 b x y ，并全部初始化为0

double **A, *b, *x, *y, **L, **U,**D;

double calculate_e = 0.0001;//默认精度为10^-4

unsigned int RANK = 4;

unsigned int makematrix()

{

unsigned int r, c;

printf("请输入矩阵行列数，用空格隔开：");

scanf_s("%d %d", &r, &c);

A = (double**)malloc(sizeof(double*)*r);//创建一个指针数组，把指针数组的地址赋值给a ,*r是乘以r的意思

for (int i = 0; i < r; i++)

A[i] = (double*)malloc(sizeof(double)*c);//给第二维分配空间

for (int i = 0; i < r; i++) {

for (int j = 0; j < c; j++)

A[i][j] = 0.0;

}

b = (double*)malloc(sizeof(double)*r);

for (int i = 0; i < r; i++)

{

b[i] = 0.0;

}

x = (double*)malloc(sizeof(double)*c);

for (int i = 0; i < c; i++)

{

x[i] = 0.0;

}

L = (double**)malloc(sizeof(double*)*r);//创建一个指针数组，把指针数组的地址赋值给a ,*r是乘以r的意思

for (int i = 0; i < r; i++)

L[i] = (double*)malloc(sizeof(double)*c);//给第二维分配空间

for (int i = 0; i < r; i++) {

for (int j = 0; j < c; j++)

L[i][j] = 0.0;

}

U = (double**)malloc(sizeof(double*)*r);//创建一个指针数组，把指针数组的地址赋值给a ,*r是乘以r的意思

for (int i = 0; i < r; i++)

U[i] = (double*)malloc(sizeof(double)*c);//给第二维分配空间

for (int i = 0; i < r; i++) {

for (int j = 0; j < c; j++)

U[i][j] = 0.0;

}

D = (double**)malloc(sizeof(double*)*r);//创建一个指针数组，把指针数组的地址赋值给a ,*r是乘以r的意思

for (int i = 0; i < r; i++)

D[i] = (double*)malloc(sizeof(double)*c);//给第二维分配空间

for (int i = 0; i < r; i++) {

for (int j = 0; j < c; j++)

D[i][j] = 0.0;

}

y = (double*)malloc(sizeof(double)*c);

for (int i = 0; i < c; i++)

{

y[i] = 0.0;

}

return r;

}

//提示用户输入一个方阵的内容 还有常数向量 计算精度

void getmatrix(void)//输入矩阵并呈现

{

printf("请按行从左到右依次输入系数矩阵A，不同元素用空格隔开\n");

for (int i = 0; i < RANK; i++)

{

for (int j = 0; jget_add)? fabs(x[i] - y[i]): get_add;

}

if (fabs(get_add) <= calculate_e)

{

printf ("迭代次数为：%d",k + 1);

break;

}

if (k == 99)//失效

{

return false;

}

}

for (int i = 0; i < RANK; i++) //交换xy

{

double temp;

temp = x[i];

x[i] = y[i];

y[i] = temp;

}

printf("求解x，解得：\n");

for (int i = 0; iget_add) ? (x[i] - y[i]) : get_add;

}

if (fabs(get_add) <= calculate_e)

{

printf("迭代次数为：%d", k + 1);

break;

}

if (k == 99)//失效

{

return false;

}

}

printf("求解x，解得：\n");

for (int i = 0; i

按设计的提示为所欲为 老老实实输入题目的系数矩阵和常数向量后，得到运行结果：

可以看到，通过不断迭代可以达到非常高的精度。