《计算方法》 李晓红等 第五章 插值 例题 代码

#include "stdafx.h"

#include <math.h>

// 插值

 

double PI()

{

    return(atan(1.)*4);

}

 

// x | pi/6     pi/4        pi/3

//------------------------------------

// y | 1/2    sqrt(2)/2     sqrt(3)/2

void LagrangeInterpolation_Order1()

{

    printf("Lagrange 插值(1阶): \n");

    int n = 1;

    double x[2],y[2];

    double phi,x0,y0,l;

    int k,i;