雅克比迭代法 c语言实现计算方法

#include “stdio.h“
#include “stdlib.h“
#include “math.h“
/*======================================================
// 函数名：jaco
// 功能描述：用雅克比迭代法解方程组
// 输入参数：a 解系数矩阵，b 常数矩阵，x返回的解向量
//           n 未知数个数，eps 精度要求。
//           iter 允许的最多迭代次数
// 返回值：整型。运行成功则返回迭代次数失败则返回0
=========================================================*/
int jaco（double *a double *bdouble *xint ndouble eps int iter）
{
    int ijpiter=0;
    double fr = 1.0 t *x2;
    if（（a==NULL）||（b==NULL）||（x==NULL））            /* 检测输入的指针是否为空*/
    {
      printf（“The pointer is NULL\n“）;
      return（0）;
    }
    x2 = （double *）malloc（n*sizeof（double））;       /* 分配空间并检测是否成功*/
    if（x2==NULL）
    {
        printf（“Memory alloc failed\n“）;
        return（0）;
    }
    for（i=0; i    if（fabs（a[i*n+i]） < 1e-12）               /* 因为要做除数，需要检查是否过小*/
    {
      printf（“A[%d%d] is too smallcannot converged“ii）;
      return（0）;
    }
    for（i=0; i