标题:迭代法求解AX=B方程组
只看楼主
小白笑小白
Rank: 1
等 级:新手上路
帖 子:4
专家分:0
注 册:2022-6-29
结帖率:100%
 问题点数:0 回复次数:1 
迭代法求解AX=B方程组

请问这个程序为什么能运行,啥结果也没有
double Jacobi(double **a, double **g, double *x1,double *x0)//雅克比迭代法
{
    double sum,err=0.0;
    x0=(double *)malloc(sizeof(double )*n7);
    x1=(double *)malloc(sizeof(double )*n7);

    for(i=0;i<=(n+1)*(n+1);++i)
    {
        x0[i]=0.0;
        x1[i]=0.0;
    }
    for (k=0;k<maxn;k++)//k为迭代次数
    {
         max_error=0.0;
        for (i=0;i<=(n+1)*(n+1)-1;i++)
        {
            sum=0.0;
            for (j=0;j<=(n+1)*(n+1)-1;j++)
            {
                if(j!=i)
                    sum+=a[i][j]*x0[j];
            }
            x1[i]=(b[i]-sum)/a[i][i];
        }
        for(i=0;i<=(n+1)*(n+1)-1;i++)
        {
            err=fabs(x1[i]-x0[i]);//绝对误差
            if(err>max_error)
                max_error=err;
                x0[i]=x1[i];
        }
        if (err<pow(10,-7))
            break;
   }
     printf("经过%d次雅克比迭代得到以下解向量\n", k+1);

}
搜索更多相关主题的帖子: for 方程组 double sum 迭代法 
2022-07-11 00:47
op123
Rank: 6Rank: 6
等 级:贵宾
威 望:21
帖 子:170
专家分:461
注 册:2022-6-4
得分:0 
就不能把所有代码和输入样例都打出来吗!
2022-07-31 14:24



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-509536-1-1.html




关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.209039 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved