trch0.c

数值分析的常用算法

  #include "stdio.h"
  #include "trch.c"
  main()
  { int n,i,j,k;
    static double t[6]={10.0,5.0,4.0,3.0,2.0,1.0};
    static double tt[6]={0.0,-1.0,-2.0,-3.0,-4.0,-5.0};
    double b[6][6];
    n=6;
    i=trch(n,t,tt,b);
    if (i>0)
      { printf("B=inv(T):\n");
        for (i=0; i<=5; i++)
          { for (j=0; j<=5; j++)
              printf("%12.6e ",b[i][j]);
            printf("\n");
          }
        printf("\n");
      }
  }