计算方法差商与牛顿插值

最新推荐文章于 2022-11-26 13:47:49 发布

Authur_gyc

最新推荐文章于 2022-11-26 13:47:49 发布

阅读量1.3k

点赞数

分类专栏： # 计算方法

本文链接：https://blog.csdn.net/WHY995987477/article/details/106174791

版权

计算方法专栏收录该内容

5 篇文章

订阅专栏

题目

在这里插入图片描述

牛顿插值

在这里插入图片描述

求差商代码

#include <iostream>
#include <algorithm>
#include <vector>
#include <map>
#include <string>
#include <string.h>
#include <cstdio>
#include <cstdlib>
#include <queue>
using namespace std;
const int maxn = 1e2;

double ans[maxn][maxn];
struct node{
    double x, fx;
}num[maxn];

int main(){
    printf("个数:");
    int n;///5
    cin >> n;
    for(int i = 0; i < n; i++){
        cin >> num[i].x;
    }
    for(int i = 0; i < n; i++){
        cin >> num[i].fx;
        ans[i][0] = num[i].fx;
    }
    for(int i = 1; i <= n; i++){
        for(int j = i; j <= n; j++)
        {
            ans[j][i] = (ans[j - 1][i - 1] - ans[j][i - 1]) / (num[j - i].x - num[j].x);
        }
    }
    for(int i = 0; i < n; i++){
        for(int j = 0; j <= i; j++){
            printf("%f ", ans[i][j]);
        }
        printf("\n");
    }
    return 0;
}
/*
5
0 1 2 3 4
0 16 46 88 0
*/