Lagrange interpolation (dup): Difference between revisions
From JSXGraph Wiki
| A WASSERMANN (talk | contribs) No edit summary | A WASSERMANN (talk | contribs) No edit summary | ||
| Line 17: | Line 17: | ||
|          p[1] = board.createElement('point', [3,-1], {style:6}); |          p[1] = board.createElement('point', [3,-1], {style:6}); | ||
|          var polynomial = function(x) { |          var polynomial = function(x) { | ||
|                 var i; | |||
|                  var y = 0.0; |                  var y = 0.0; | ||
|                 var xc = []; | |||
|                 for (var i=0;i<p.length;i++) { | |||
|                   xc[i] = p[i].X(); | |||
|                 } | |||
|                  for (var i=0;i<p.length;i++) { |                  for (var i=0;i<p.length;i++) { | ||
|                    var t = p[i].Y(); |                    var t = p[i].Y(); | ||
|                    for (var k=0;k<p.length;k++) { |                    for (var k=0;k<p.length;k++) { | ||
|                      if (k!=i) { |                      if (k!=i) { | ||
|                        t *= (x- |                        t *= (x-xc[k])/(xc[i]-xc[k]); | ||
|                      } |                      } | ||
|                    } |                    } | ||
Revision as of 22:07, 10 December 2008
Constructs a polynomial of degree n through n+1 given points. Points can be added by clicking on "Add point".
        board = JXG.JSXGraph.initBoard('box', {originX: 250, originY: 250, unitX: 50, unitY: 25});
        // Axes
        b1axisx = board.createElement('axis', [[0,0], [1,0]], {});
        b1axisy = board.createElement('axis', [[0,0], [0,1]], {});
        var p = [];
        p[0] = board.createElement('point', [-1,2], {style:6});
        p[1] = board.createElement('point', [3,-1], {style:6});
        
        var polynomial = function(x) {
                var y = 0.0;
                for (var i=0;i<p.length;i++) {
                  var t = p[i].Y();
                  for (var k=0;k<p.length;k++) {
                    if (k!=i) {
                      t *= (x-p[k].X())/(p[i].X()-p[k].X());
                    }
                  }
                  y += t;
                }
                return y;
            };
        graph = board.createElement('curve', ['x', polynomial, 'x', -10, 10], {curveType:'graph'});
        function addPoint() {
          p.push(board.createElement('point',[(Math.random()-0.5)*10,(Math.random()-0.5)*3],{style:6}));
          board.update();
        }
