newton chazhi.c

C语言程 反对萨芬记录卡上的交付离开的撒娇发了啊

#include<stdio.h>
#include<conio.h>
#include<alloc.h>
#define N 4
void Difference(float *x,float *y,int n)
 {   float *f;   int k,i;
   f=(float *)malloc(n*sizeof(float));
   for(k=1;k<=n;k++)
     {  f[0]=y[k];
         for(i=0;i<k;i++)
           f[i+1]=(f[i]-y[i])/(x[k]-x[i]) ;
             y[k]=f[k];}
     }
main()
{ int i;
    float varx1=0.895,b,varx2=0.596;
    float x[N+1]={0.4,0.55,0.65,0.8,0.9};
    float y[N+1]={0.41075,0.57815,0.69675,0.88811,1.02652}   ;
    Difference(x,(float *)y,N);
    clrscr();      b=y[N];
    for(i=N-1;i>=0;i--)  b=b*(varx1-x[i])+y[i];
     printf("NInterp(%f)=%f\n",varx1,b);
    for(i=N-1;i>=0;i--)  b=b*(varx2-x[i])+y[i];
     printf("NInterp(%f)=%f",varx2,b);
   getch();
 }

NInterp(0.895000)=1.019368
NInterp(0.596000)=0.632016