小生编了个数值程序，求教

[jinjiachen]

#include<stdio.h>
int main(void)
{
int i,max=1000;
float epsi=0.000001,x[3],y[3],a[4],b[4],c[4];
printf("形如a1x1+a2x2+a3x3=a4\n b1x1+b2x2+b3x3=b4\n c1x1+c2x2+c3x3=c4\n");
printf("请输入系数");
for(i=0;i<4;i++)
{
scanf("%f",&a);
}

for(i=0;i<4;i++)
{
scanf("%f",&b);
}

for(i=0;i<4;i++)
{
scanf("%f",&c);
}
for(i=0;i<3;i++)
{
scanf("%f",&x);
y=x;
}
for(i=1;i<=max;i++)
{
x[0]=-a[1]/a[0]*x[1]-a[2]/a[0]*x[2]+a[3]/a[0];
x[1]=-b[0]/b[1]*x[0]-b[2]/b[1]*x[2]+b[3]/b[1];
x[2]=-c[0]/c[2]*x[0]-c[1]/c[2]*x[1]+c[3]/c[2];
if((x[0]-y[0])<epsi&&(x[1]-y[1])<epsi&&(x[2]-y[2])<epsi)
printf("x1=%f x2=%f x3=%f ",x[0],x[1],x[2]);
else
{
for(i=0;i<3;i++)
y=x;
}
if(i==max)
printf("发散");
}
return(0);
}
以上是我用c编的一个三元一次方程组的数值解法，用的是高斯迭代，有些地方不够完善，有待改进，不过我现在的疑问是：
开头的几个循环是输入系数，不过我运行的时候一直让我输入，明明系数几经满了，还在让我输入。。。。

[jinjiachen]

问题解决了，嘿嘿，我刚开始认为是程序的问题，原来是我写的程序不够好，导致运行出错，哈哈

[jinjiachen]

#include<stdio.h>
int main(void)
{
int i,max=10;
float epsi=0.000001,x[3],y[3],a[4],b[4],c[4];
printf("形如a1x1+a2x2+a3x3=a4\n b1x1+b2x2+b3x3=b4\n c1x1+c2x2+c3x3=c4\n");
printf("请输入系数\n");
for(i=0;i<4;i++)
{
scanf("%f",&a);
}

for(i=0;i<4;i++)
{
scanf("%f",&b);
}

for(i=0;i<4;i++)
{
scanf("%f",&c);
}
printf("请赋初值");
for(i=0;i<3;i++)
{
scanf("%f",&x);
y=x;
}
printf("输入完毕，请等候\n");
for(i=1;i<max;i=i+1)
{
x[0]=-a[1]/a[0]*x[1]-a[2]/a[0]*x[2]+a[3]/a[0];
x[1]=-b[0]/b[1]*x[0]-b[2]/b[1]*x[2]+b[3]/b[1];
x[2]=-c[0]/c[2]*x[0]-c[1]/c[2]*x[1]+c[3]/c[2];
if((x[0]-y[0])<epsi&&(x[1]-y[1])<epsi&&(x[2]-y[2])<epsi)
{
printf("x1=%f x2=%f x3=%f \n",x[0],x[1],x[2]);
break;
}
else
{
for(i=0;i<3;i++)
y=x;
}
if(i==max)
{
printf("发散");
}
}
return(0);
}
改进了一下，收敛得方程可以显示解了，不过发散怎么也不行，改过很多次，但还是不行，各位大虾帮帮我

[b33e]

这东西，matlab搞定

[Large_Windmill]

Matlab几句话搞定，比这好弄多。

[月下叹逍遥]

建议使用牛顿下山迭代，用C大才小用了。。。

[jinjiachen]

这东西，matlab搞定

嗯嗯，我知道那东西强，不过我的目的并不是要解方程，一是对计算方法的理解，二是对c的训练

[jinjiachen]

建议使用牛顿下山迭代，用C大才小用了。。。

牛顿下山要求导，通用性不好，我只是想自己编程实现一些常用算法，你能帮帮我看下哪里不足吗？我的现在的问题是：收敛的方程可以输出解，发散的就一直在算，不能输出“发散”，我觉得好像逻辑没有错啊

[lainme]

for(i=0;i<3;i++)
y=x;

指标重复了吧

目测 if (i==max) 这段要放到那个 for 循环外边。缩进一下就看出来了。

[jinjiachen]

for(i=0;i<3;i++)
y=x;

指标重复了吧

目测 if (i==max) 这段要放到那个 for 循环外边。缩进一下就看出来了。

对哦，i重复了，我去修改下试试，谢谢