数值计算程序大放送-线性代数方程组

Word-BREAK: break-all; TEXT-INDENT: 24px; LINE-HEIGHT: normal; HEIGHT: 200px; WORD-WRAP: break-word" onload="this.style.overflowX='auto';">

选自<<徐世良数值计算程序集(C)>>

每个程序都加上了适当地注释，陆陆续续干了几个月才整理出来的啊。

今天都给贴出来了

#include "stdlib.h"
#include "math.h"
#include "stdio.h"
// 全选主元高斯消去法
//a-n*n 存放方程组的系数矩阵，返回时将被破坏
//b-常数向量
//x-返回方程组的解向量
//n-存放方程组的阶数
//返回0表示原方程组的系数矩阵奇异
int cagaus(double a[],double b[],int n,double x[])
{
int *js,l,k,i,j,is,p,q;
    double d,t;
    js=malloc(n*sizeof(int));
    l=1;
    for (k=0;k<=n-2;k++)
{
  d=0.0;
        for (i=k;i<=n-1;i++)
  {
   for (j=k;j<=n-1;j++)
            {
    t=fabs(a[i*n+j]);
    if (t>d)
    {
     d=t;
     js[k]=j;
     is=i;
    }
            }
  }
  if (d+1.0==1.0)
  {
   l=0;
  }
  else
  {
   if (js[k]!=k)
   {
    for (i=0;i<=n-1;i++)
    {
     p=i*n+k;
     q=i*n+js[k];
     t=a[p];
     a[p]=a[q];
     a[q]=t;
    }
   }
   if (is!=k)
   {
    for (j=k;j<=n-1;j++)
    {
     p=k*n+j;
     q=is*n+j;
     t=a[p];
     a[p]=a[q];
     a[q]=t;
    }
    t=b[k];
    b[k]=b[is];
    b[is]=t;
   }
  }
  if (l==0)
  {
   free(js);
   printf("fail\n");
   return(0);
  }
  d=a[k*n+k];
  for (j=k+1;j<=n-1;j++)
  {
   p=k*n+j;
   a[p]=a[p]/d;
  }
  b[k]=b[k]/d;
  for (i=k+1;i<=n-1;i++)
  {
   for (j=k+1;j<=n-1;j++)
   {
    p=i*n+j;
    a[p]=a[p]-a[i*n+k]*a[k*n+j];

[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] 下一页

复制本页网址和标题，发送给你QQ/Msn的好友一起分享

上一篇:浙江台州互联网行业浅析－－问题篇

下一篇:数值计算程序大放送-特征值和特征向量