zhuigan.c

来自「两种曲线拟合的算法」· C语言 代码 · 共 27 行

#include <stdio.h>
#define N 4
main()
{float sz[N][N+1],p[N],q[N];
 int i,j;
 printf("请输入增广矩阵：\n");                  //输入
 for(i=0;i<N;i++)
     for(j=0;j<=N;j++)
         scanf("%f",&sz[i][j]);
 p[0]=sz[0][0];
 q[0]=sz[0][1]/p[0];
 for(i=1;i<N;i++)                             //分别求p[i] q[i]
     {p[i]=sz[i][i]-sz[i][i-1]*q[i-1];
      if(i<(N-1))
      q[i]=sz[i][i+1]/p[i];
     }
 p[0]=sz[0][N]/p[0];
 for(i=1;i<N;i++)                              //求y[i]
     p[i]=(sz[i][N]-sz[i][i-1]*p[i-1])/p[i];
 q[N-1]=p[N-1];
 for(i=N-2;i>=0;i--)
     q[i]=p[i]-q[i]*q[i+1];
 printf("结果为：");
 for(j=0,i=N;j<N;i--,j++)
     printf("X%d=%5.6f ",i,q[j]);
 printf("\n");
}