6gmqr.c

c语言常用算法

  #include "4maqr.c"
  #include "stdlib.h"
  int gmqr(a,m,n,b,q)
  int m,n;
  double a[],b[],q[];
  { int i,j;
    double d,*c;
    int maqr(double [],int,int,double []);
    c=malloc(n*sizeof(double));
    i=maqr(a,m,n,q);
    if (i==0) { free(c); return(0);}
    for (i=0; i<=n-1; i++)
      { d=0.0;
        for (j=0; j<=m-1; j++)
          d=d+q[j*m+i]*b[j];
        c[i]=d;
      }
    b[n-1]=c[n-1]/a[n*n-1];
    for (i=n-2; i>=0; i--)
      { d=0.0;
        for (j=i+1; j<=n-1; j++)
          d=d+a[i*n+j]*b[j];
        b[i]=(c[i]-d)/a[i*n+i];
      }
    free(c); return(1);
  }