来自「收集的C语言算法程序」· 代码 · 共 1 行

#include<iostream.h>#include<math.h>//确定要计算的角度值#define N 0.57891//利用宏定义线性插值函数#define L(a) (a-x_1)/(x_0-x_1)*y_0+(a-x_0)/(x_1-x_0)*y_1//利用宏定义函数，该函数用以将x0与x1逐步接近#define X(M,n) int(M*pow(10,n)/10)/pow(10,n-1)void main(){       double x_0=0.5,x_1=0.6,y_0=0.47943,y_1=0.56464,temp=N,temp1;       for(int i=3;i<=6;i++)       {              temp1=X(temp,i) ;              y_0=L(temp1);              if(i==6)                     break;              x_0=X(temp,i);       }       //输出运算结果       cout<<"经用线性插值法计算知sin(0.57891)的近似值为："<<y_0<<endl;       //输出截断误差限       cout<<"截断误差限为："<<sin(x_1)*(fabs((N-x_0)*(N-x_1))/2)<<endl;}