【ソースコード】
#include <stdio.h>
#include <stdlib.h>
#include <math.h>

/* Interpolation with Lagrange method */

main()
{

  double  x[10],y[10];
  double  xx,yy;
  double  pn,li;
  double  deg=3.141592/180.0;

  int i,j;
  int n;

  n=10;
  for (i=0;i<n;i++){
    x[i]=(double)i *10.0*deg;  y[i]=sin(x[i]);
  }

  xx=45.0*deg;
  yy=sin(xx);

  pn=0.0;
  for (i=0;i<n;i++){
    li=1.0;
    for (j=0;j<n;j++){
      if(i != j){
        li *= (xx-x[j])/(x[i]-x[j]);
      }
    }
    pn+=y[i]*li;
  }

  printf(" 45 deg:%10.5f rad   Pn:%7.5f   True:%7.5f \n",xx,pn,yy);
  
  return(0);
}


【実行結果】
　45 deg:   0.78540 rad   Pn:0.70711   True:0.70711