10个重要的算法C语言实现源代码(8-9-10-----秦九昭和幂法和高斯塞德尔)

#include <math.h>
float qin(float a[],int n,float x)
{    float r=0;
int i;
for(i=n;i>=0;i--)
r=r*x+a[i];
return r;
}
main()
{    float a[50],x,r=0;
int n,i;
do
{    printf("Input frequency:");
scanf("%d",&n);
}
while(n<1);
printf("Input value:");
for(i=0;i<=n;i++)
scanf("%f",&a[i]);
printf("Input frequency:");
scanf("%f",&x);
r=qin(a,n,x);
printf("Answer:%f",r);
getch();
}

#include<stdio.h>
#include<math.h>
#define N 100
#define e 0.00001
#define n 3
float x
={0,0,1};
float a

={{2,3,2},{10,3,4},{3,6,1}};
float y
;
main()
{ int i,j,k;
float xm,oxm;
oxm=0;
for(k=0;k<N;k++)
{ for(j=0;j<n;j++)
{ y[j]=0;
for(i=0;i<n;i++)
y[j]+=a[j][i]*x[i];
}
xm=0;
for(j=0;j<n;j++)
if(fabs(y[j])>xm) xm=fabs(y[j]);
for(j=0;j<n;j++)
y[j]/=xm;
for(j=0;j<n;j++)
x[j]=y[j];
if(fabs(xm-oxm)<e)
{ printf("max:%f/n/n",xm);
printf("v[i]:/n");
for(k=0;k<n;k++) printf("%f/n",y[k]);
break;
}
oxm=xm;
}
getch();
}

#include<math.h>
#include<stdio.h>
#define N 20
#define M 99
float a

;
float b
;
int main()
{    int i,j,k,n;
float sum,no,d,s,x
;
printf("/nInput dim of n:");
scanf("%d",&n);
if(n>N)
{ printf("Faild! Check if 0<n<N!/n "); getch();
return 1;
}
if(n<=0)
{ printf("Faild! Check if 0<n<N!/n ");getch();return 1;}
printf("Input a[i,j],i,j=0…%d:/n",n-1);
for(i=0;i<n;i++)
for(j=0;j<n;j++)
scanf("%f",&a[i][j]);
printf("Input b[i],i=0…%d:/n",n-1);
for(i=0;i<n;i++) scanf("%f",&b[i]);
for(i=0;i<n;i++) x[i]=0;
k=0;
printf("/nk=%dx=",k);
for(i=0;i<n;i++) printf("%12.8f",x[i]);
do
{ k++;
if(k>M){printf("/nError!/n”);getch();}
break;
}
no=0.0;
for(i=0;i<n;i++)
{ s=x[i];
sum=0.0;
for(j=0;j<n;j++)
if (j!=i) sum=sum+a[i][j]*x[j];
x[i]=(b[i]-sum)/a[i][i];
d=fabs(x[i]-s);
if (no<d) no=d;
}
printf("/nk=%2dx=",k);
for(i=0;i<n;i++)   printf("%f",x[i]);
}
while (no>=0.1e-6);
if(no<0.1e-6)
{ printf("/n/n answer=/n");
printf("/nk=%d",k);
for (i=0;i<n;i++)
printf("/n x[%d]=%12.8f",i,x[i]);
}
getch();
}