newcsfy.c

来自「拉格朗日插值 用拉格朗日插值来计算函数值 使用说明： 一般」· C语言 代码 · 共 30 行

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

void Newcsfy(double x[],double y[],int n,double R[])
{  double y1[10][10];
   double x1[10],Ax[10];
   double x2,y2,Nx,wn;
   int i,j,k,t;
   for(i=0;i<n;i++)                                    /*将二元数组y[n][n]赋零;*/
      { for(j=0;j<n;j++)
            {y1[i][j]=0;}
      }
   for(i=0;i<n;i++)                                    /*将y值赋给二元数组y[0][n];*/
      { y1[0][i]=y[i];
        x1[i]=x[i];}
   k=1;
   for(j=0;j<n;j++)
      {
       {for(i=j;i<n;i++)
        y1[j+1][i+1]=(y1[j][i+1]-y1[j][i])/(x[i-j+k]-x[i-j]);
        }
        k++;
    }
   for(i=0;i<n;i++)                                   /*Ax是N(x)函数方程系数*/
      {Ax[i]=y1[i][i];
       }
   R[0]=Ax[n-1];
 
}