guass-xiao.c

数值分析的几个C实现,比较小.

#include<stdio.h>

main()
{int k,i,j,l,n=3;double z;
 double a[3][3]={{2,-1,3},{4,2,5},{1,2,0}},b[3]={1,4,7},d,t;



 for(k=0;k<n;k++)
 { for(j=k+1;j<n;j++)
     a[k][j]/=a[k][k];
   b[k]/=a[k][k];
   for(i=k+1;i<n;i++)
   { for(j=k+1;j<n;j++)
        a[i][j]-=a[i][k]*a[k][j];
     b[i]-=a[i][k]*b[k];
   }
  }
 for(i=n-2;i>=0;i--)
 { z=0;
   for(j=i+1;j<n;j++)
     z+=a[i][j]*b[j];
   b[i]-=z;
  }
 for(i=0;i<n;i++)
  printf("%.6f\n",b[i]);

 getch();
 }