C++でのラグランジュの補間

このチュートリアルでは、lagrangesの補間式の結果を見つけるプログラムを作成します。

プログラムのロジックを書く必要はありません。数式をコードに変換するだけです。コードを見てみましょう。

例

#include<bits/stdc++.h>
using namespace std;
struct Data {
   int x, y;
};
double interpolate(Data function[], int xi, int n) {
   double result = 0;
   for (int i = 0; i < n; i++) {
      double term = function[i].y;
      for (int j = 0; j < n; j++) {
         if (j != i) {
            term = term * (xi - function[j].x) / double(function[i].x - function[j].x);
         }
      }
      result += term;
   }
   return result;
}
int main() {
   Data function[] = {{0,3}, {1,2}, {6,9}, {10,17}};
   cout << interpolate(function, 3, 5) << endl;
   return 0;
}

出力

上記のコードを実行すると、次の結果が得られます。

3

結論

チュートリアルに質問がある場合は、コメントセクションにそのことを記載してください。