©2020-2024 ioDraw All rights reserved
实验三 牛顿插值多项式
实验内容：天安门广场升旗的时间是日出的时刻，而降旗的时间是日落时分，根据天安门广场管理委员会的公告，某年10月份升降旗的时间如下：
日期 1 15 22
升旗 6:09 6:23 6:31
降旗 17:58 17:36 17:26
根据上述数据构造Newton插值多项式，计算当年10月8日北京市的日照时长。
```c #include<stdio.h> #include<math.h> #define n 3 void main() { int
i,j,X=8,x[n]={1,15,22}; float
y[n]={17-6+(58-9)/60.0,17-6+(36-23)/60.0,17-6+(26-31)/60.0},
time=0.0,z[n][n];//z里面存储的是差商值 for(i=0;i<n;i++)//完成了0阶差商的赋值 z[i][0]=y[i];
for(j=1;j<n;j++)//j是控制我们计算几阶差商的，如果看成是二维数组的列号也可以
for(i=j;i<n;i++)//i只需要计算对角线下面的，也就是i>j的情况 z[i][j]=( ); //这个就是我们转换到数组后差商的计算公式
time=( ); //time里面就是我们最后的结果，它的形式就是newton插值多项式
printf("10月8号的日照时长为：%d小时%.0f分钟\n", (int)floor(time),60*fmod(time,1.0)); }