拉格朗日函数c++

浏览数：23 / 时间：2015年06月08日 
  1. #include <iostream>
  2. using namespace std;
  3. //拉格朗日插值求解方程
  5. double * xs;   //all input x
  6. double * ys;   //all input y
  7. int n;      //size
  9. //1.1275 1.1503 1.1735 1.1972
  10. //0.1191 0.13954 0.15932 0.17903
  11. //1.130
  12. void init()
  13. {
  14.     cout << "please input n " << endl;
  15.     cin >> n;
  16.     xs = new double[n];
  17.     ys = new double[n];
  19. //input x
  20.     cout << "please input the x‘s value !" << endl;
  21.     for(int i=0; i<n; i++)
  22.     {
  23.         cin >> xs[i];
  24.     }
  26. //input y
  27.     cout << "please input the y‘s value !" << endl;
  28.     for(int i=0; i<n; i++)
  29.     {
  30.         cin >> ys[i];
  31.     }
  32. }
  34. double everyItem(int num, double x)
  35. {
  36.     double y = ys[num];
  37.     double up = 1;
  38.     double down = 1;
  39.     for(int i=0; i<n; i++)
  40.     {
  41.         if(i != num)
  42.         {
  43.             up *= (x - xs[i]);
  44.             down *= (xs[num] - xs[i]);
  45.         }
  46.     }
  48.     return y*up/down;
  49. }
  51. double lagrange(double x)
  52. {
  53.     double total = 0;
  54.     for(int i=0; i<n; i++)
  55.     {
  56.         total = total + everyItem(i , x);
  57.     }
  59.     return total;
  60. }
  64. int main()
  65. {
  67.     init();
  68.     double input;
  69.     cout << "please input you x !" <<endl;
  70.     cin >> input;
  72.     double result = lagrange(input);
  73.     cout << "the result is " << result << endl;
  74.     return 0;
  75. }





附件列表

     

    郑重声明：本站内容如果来自互联网及其他传播媒体，其版权均属原媒体及文章作者所有。转载目的在于传递更多信息及用于网络分享，并不代表本站赞同其观点和对其真实性负责，也不构成任何其他建议。

    标签： 编程语言 Java C C# Javascript js Python Ruby Spring 程序 VB

    相关文章

    随机文章

    您可能还喜欢

    您可能还喜欢