高斯赛德尔迭代c语言

‘壹’ 19.7192x^3 +295.788x^2 = 673.7341 一元三次方程求解，最好过程

假如给我们一个一般的三次方程： ax3+3bx2+3cx+d=0 （1） 如果令 x=y-b/a 我们就把方程（1）推导成 y3+3py+2q=0 （2） 其中 3p=c/a-b2/a2，2q=2b3/a3-3bc/a2+d/a 。 借助于等式 y=u-p/u 引入新变量u 。把这个表达式带入（2），得到： （u3）2+2qu3-p3=0 （3） 由此得 u3=-q±√（q2+p3）， 于是 y=3√（-q±√（q2+p3））-p/3√（-q±√（q2+p3）） 。 =3√（-q+√（q2+p3））+3√（-q-√（q2+p3）） 。 （最后这个等式里的两个立方根的积等于-p 。）
用此法代入求解，过程复杂，不宜展示

‘贰’ 谁能帮我设计一个C语言程序,求解方程组的,急啊,哪位高手帮下咯,万分感谢

假定你要的是线性方程组，下面的
float *GauseSeidel(float *a,int n)
是高斯赛德尔法求解线性方程组的通用子程序。
N 是 迭代次数极限。

main()里写了调用的例子。

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define N 100
float *GauseSeidel(float *a,int n)
{
int i,j,nu=0;
float *x,dx;
x=(float *)malloc(n*sizeof(float));
for(i=0;i<=n-1;i++)
x[i]=0.0;
do {
for(i=0;i<=n-1;i++) {
float d=0.0;
for(j=0;j<=n-1;j++)
d+=*(a+i*(n+1)+j)*x[j];
dx=(*(a+i*(n+1)+n)-d)/(*(a+i*(n+1)+i));
x[i]+=dx;
}
if(nu>=N)
{
printf("iter divergence\n");
exit(1);
}
nu++;
}
while(fabs(dx)>1e-6);
return x;
}
void main()
{
int i;
float *x;
float c[12]={5,2,1,8,2,8,-3,21,1,-3,-6,1};
float *GauseSeidel(float *,int);
x=GauseSeidel(c,3);
for(i=0;i<=2;i++)
printf("x[%d]=%f\n",i,x[i]);
getch();
}