帮忙看一段代码！谢谢!

flameboy

请教一下异步迭代法求解线性方程组的程序。
以下是我写的程序，算不出结果，但我不知错误在哪里。
初学者刚学写代码，请大家别笑话啊。
方程组形式:
a[i][j]*x[i]=b[i]

#include<iostream>
using namespace std;

const int size=2;
const int times =10;

float a[size][size]={1,1,1,2};
float b[size]={1,2,};

int main(){
        float* x=new float[size];
        float* sum=new float[size];
        for(int i=0;i<size;i++){
                x[i]=i;
                sum[i]=0;
        }
        for(int i=0;i<times;i++){
                for(int row=0;row<size;row++){
                        for(int line=0;line<size;line++){
                                if(line=row)
                                        continue;
                                sum[row]+=a[row][line]*x[row];
                        }
                }
                for(row=0;row<size;row++)
                        x[row]=(b[row]-sum[row])/a[row][row];
        }
        for(i=0;i<size;i++)
                cout<<"x"<<i<<"= "<<x[i]<<" ";
        delete []x;
}

boywhp

好恐怖的代码啊？最好说明你的酸法的思路啊

jiacat

if(line=row)－＿－！这个是什么意思？

boywhp

if(line=row)
continue;
这个我用VB的都懂呢！

jiacat

楼上的比较强-_-!

Hush

确实比较强...不亏是用VB的 Orz

boywhp

哈哈！过奖啊 但是好象还不服气呢！嘿嘿

flameboy

算法思路：
row代表行，line 代表列；


#include<iostream>
using namespace std;

const int size=2;
const int times =10;

float a[size][size]={1,1,1,2};//系数矩阵2×2；
float b[size]={1,2,};//方程形式是：ax=b;

int main(){
float* x=new float[size];
float* sum=new float[size];//一个存放求和结果的数组。
for(int i=0;i<size;i++){
x[i]=i;//赋初值
sum[i]=0;
}
for(int i=0;i<times;i++){            //迭代次数
for(int row=0;row<size;row++){   
//从第一行开始，处于每一列的系数与相应的x相乘，并作和。但不包括  
//a[row][line](row=line)与x[row]的乘积。因为x[row]是每次迭代要解出
//的结果，作为下一次迭代时的待入值。这就是sum[row]。
//而x[row]=(b[row]-sum[row])/a[row][row];
for(int line=0;line<size;line++){
if(line=row)
continue;
sum[row]+=a[row][line]*x[row];
}
}
for(row=0;row<size;row++)
x[row]=(b[row]-sum[row])/a[row][row];
}
for(i=0;i<size;i++)
cout<<"x"<<i<<"= "<<x[i]<<" ";
delete []x;
}

boywhp

看见了没有！大家现在可以笔试了，呵呵！我数学不懂，只能看看语法了

DreamSnow

数学学的不是很好，研究算法实在是费尽，找点语法错误吧！
if(line=row)
    continue;
我想这个应该不对吧？line=row是一个赋值，这样写每次运行到这里都会continue的，下面的代码段是执行不到的，应该是这样吧？
if(line==row)
    continue;