#include<iostream>
using namespace std;
const int n=3;
void jacobi();
double A[n][n]={{10,-1.0,-2.0},{-1.0,10.0,-2.0},{-1.0,-1.0,5.0}};///在 main() 函数前写,不然在调用jacobi()时会出错
double B[n]={72,83,42};
int main()
{
jacobi();
}
void jacobi()
{
double X[n]={0,0,0,};
for (int k=0;k<10;k++)
{
for(int i=0;i<n;i++)
{
double sum=0;
for(int j=0;j<n;j++)
{
if(j==i) continue;
sum+=A[i][j]*X[j];
}
X[i]=(B[i]-sum)/A[i][i];
}
}
for (int i=0;i<n;i++)
cout<<X[i]<<'\t';
cout<<endl;
}C++理解全局变量在调用函数中的应用
最新推荐文章于 2023-02-28 17:04:28 发布
本文介绍了一个使用C++实现的雅可比迭代法求解线性方程组的例子。该程序定义了系数矩阵A和常数向量B,并通过迭代计算得到未知数的近似值。
294

被折叠的 条评论
为什么被折叠?



