用C语言编程解N元方程组
已经做过了,我都留着呢!
哈哈哈
#include<math.h>
#define M 10
float fdel(float *p,int n)
{ int i;
float sum=0;
for(i=0;i<n;i++)
sum+=*(p+i);
sum/=n;
return (fabs(sum));
}
void chu(float *p)
{int i;
for(i=0;i<M;i++)
*(p+i)=999;
return ;}
void print(float *p,int n)
{int i;
for(i=0;i<n;i++)
printf("x%d=%f ",i+1,*(p+i));
return ;}
/**void scan(float **p,int n)
{int i,j;
for(i=0;i<n;i++)
for(j=0;j<n;j++)
{printf(" a%d%d=",i+1,j+1);
scanf("%f",*(p+i)+j);}
return ;}**/
void den(float *a,float *b,int n)
{int i;
for(i=0;i<n;i++)
*(a+i)=*(b+i);
return ;}
void jie(float *p0,float a[M][M],float *p,float *b,int n)
{float sum=0;
int i,j;
for(i=0;i<n;i++)
{ for(j=0;j<n;j++)
{if(i!=j) sum=sum+a[i][j]*p0[j];
else {n=a[i][j];continue;}
}
*(p+i)=*(b+i)-sum/n;}
}
main()
{ float x0[M],xi[M][M],b[M],x[M];
int i,j,n;
chu(x0);
chu(b);
chu(x);
for(i=0;i<M;i++)
{printf("x0=%f ",x0[i]);
printf("b=%f ",b[i]);
printf("x=%f ",x[i]);}
for(i=0;i<M;i++)
for(j=0;j<M;j++)
xi[i][j]=0;
printf("shu ru wuizhishu geshu \nn=");
scanf("%d",&n);
printf( " shu ru fangchengzu gege xishu a[i][j]=\n");
for(i=0;i<n;i++)
for(j=0;j<n;j++)
{ printf(" a%d%d=",i+1,j+1);
scanf("%f",&xi[i][j]); }
printf("x[1][1]=%f",xi[1][1]);
for(i=0;i<n;i++)
print(*(xi+i),n);
for(i=0;i<n;i++)
printf("\nabs=%f",fdel(*(xi+i),n));
}