各位大虾帮小弟看看这个程序哪有问题（用高斯法解放程）

[coldfishc]

这是一道用高斯法解放程的题
我写的程序如下：
#include<iostream.h>
double fun(double a,double b);
void main()
{
int row,col;
double times;
double
matrix[4][5]={1,2,-12,8,275,4,7,-2,4,-3,7,9,5,11,6,-12,-8,3,49};
double x[4]={0,0,0,0};


row=0,col=0;
for(int i=0;i<4;++i)
{
times=fun(matrix[row][col],matrix[row+i][col]);
for(int j=0;i<4;++j)
{
matrix[row+j][col]-=matrix[row+j][col+1]*times;
}
}


double othsum=0;
for(i=0;i<4;++i)
{
x[3-i]=(matrix[3-i][5]-othsum);
for(int j=4-i;j<4;++j)
{
othsum+=(matrix[3-i][j]*(x[4-i]));
}
}



for(i=0;i<4;++i)
{
cout<<x[i];
}
}


double fun(double a,double b)
{
return (b/a);
}

可以编译通过。
qq:154584468

[izdlywn]

建议你单步运行一下

我能看出的错误:

你定义是:double
matrix[4][5]={1,2,-12,8,275,4,7,-2,4,-3,7,9,5,11,6,-12,-8,3,49};
double x[4]={0,0,0,0};

那么你的matrix[3][4]就是最后一个元素.

但是你的: for(i=0;i<4;++i)
{
x[3-i]=(matrix[3-i][5]-othsum);
for(int j=4-i;j<4;++j)
{
othsum+=(matrix[3-i][j]*(x[4-i]));
}
}
中的: x[3-i]=(matrix[3-i][5]-othsum); 这句

matrix[3-i][5]必然是越间的.

[freefalcon]

row=0,col=0;
for(int i=0;i<4;++i)
{
times=fun(matrix[row][col],matrix[row+i][col]);
for(int j=0;j<5;++j) // <<<<<<<<<<<<这里你原来是i<4，估计是笔误吧，:)

[coldfishc]

谢谢

[Kudeet]

你具体出的是什么问题