重要的算法C语言实现源代码3
3.高斯列主元消去法,求解其次线性方程组C/C++ code
#include<stdio.h>
#include <math.h>
#define N 20
int main()
{ int n,i,j,k;
int mi,tmp,mx;
float a,b,x;
printf("\nInput n:");
scanf("%d",&n);
if(n>N)
{ printf("The input n should in(0,N)!\n");
getch();
return 1;
}
if(n<=0)
{ printf("The input n should in(0,N)!\n");
getch();
return 1;
}
printf("Now 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);}
printf("Now input b(i),i,j=0...%d:\n",n-1);
for(i=0;i<n;i++)
scanf("%f",&b);
for(i=0;i<n-2;i++)
{ for(j=i+1,mi=i,mx=fabs(a);j<n-1;j++)
if(fabs(a)>mx)
{ mi=j;
mx=fabs(a);
}
if(i<mi)
{ tmp=b;b=b;b=tmp;
for(j=i;j<n;j++)
{ tmp=a;
a=a;
a=tmp;
}
}
for(j=i+1;j<n;j++)
{ tmp=-a/a;
b+=b*tmp;
for(k=i;k<n;k++)
a+=a*tmp;
}
}
x=b/a;
for(i=n-2;i>=0;i--)
{ x=b;
for(j=i+1;j<n;j++)
x-=a*x;
x/=a;
}
for(i=0;i<n;i++)
printf("Answer:\n x[%d]=%f\n",i,x);
getch();
return 0;
}
#include<math.h>
#include<stdio.h>
#define NUMBER 20
#define Esc 0x1b
#define Enter 0x0d
float A ,ark;
int flag,n;
exchange(int r,int k);
float max(int k);
message();
main()
{
float x;
int r,k,i,j;
char celect;
clrscr();
printf("\n\nUse Gauss.");
printf("\n\n1.Jie please press Enter.");
printf("\n\n2.Exit press Esc.");
celect=getch();
if(celect==Esc)
exit(0);
printf("\n\n input n=");
scanf("%d",&n);
printf(" \n\nInput matrix A and B:");
for(i=1;i<=n;i++)
{
printf("\n\nInput a%d1--a%d%d and b%d:",i,i,n,i);
for(j=1;j<=n+1;j++) scanf("%f",&A);
}
for(k=1;k<=n-1;k++)
{
ark=max(k);
if(ark==0)
{
printf("\n\nIt's wrong!");message();
}
else if(flag!=k)
exchange(flag,k);
for(i=k+1;i<=n;i++)
for(j=k+1;j<=n+1;j++)
A=A-A*A/A;
}
x=A/A;
for( k=n-1;k>=1;k--)
{
float me=0;
for(j=k+1;j<=n;j++)
{
me=me+A*x;
}
x=(A-me)/A;
}
for(i=1;i<=n;i++)
{
printf(" \n\nx%d=%f",i,x);
}
message();
}
exchange(int r,int k)
{
int i;
for(i=1;i<=n+1;i++)
A=A;
for(i=1;i<=n+1;i++)
A=A;
for(i=1;i<=n+1;i++)
A=A;
}
float max(int k)
{
int i;
float temp=0;
for(i=k;i<=n;i++)
if(fabs(A)>temp)
{
temp=fabs(A);
flag=i;
}
return temp;
}
message()
{
printf("\n\n Go on Enter ,Exit press Esc!");
switch(getch())
{
case Enter: main();
case Esc: exit(0);
default:{printf("\n\nInput error!");message();}
}
}
关于本帖内容欢迎大家踊跃讨论，与在下交流！
页:
[1]