摘要

描述Lagrange polynomial.svg	English:   本W3C状态不明的矢量图使用Mathematica创作。 This image shows, for four points ((−9, 5), (−4, 2), (−1, −2), (7, 9)), the (cubic) interpolation polynomial L(x) (dashed, black), which is the sum of the scaled basis polynomials y1ℓ1(x), y2ℓ2(x), y3ℓ3(x) and y4ℓ4(x) (colors). The interpolation polynomial passes through all four control points, and each scaled basis polynomial passes through its respective control point and is 0 where x corresponds to the other three control points.
日期	2015年12月26日
来源	Self-made, based on Image:Lagrangepolys.png
作者	User:Glosser.ca
其他版本	Image:Lagrangepolys.png

Source

colors = ColorData[97] /@ Range[4];
pts = {{\[Minus]9, 5}, {\[Minus]4, 2}, {\[Minus]1, \[Minus]2}, {7, 9}};
{xs, ys} = Transpose[pts];

Basis[pts_, x_] := Times @@ (x - DeleteCases[pts, #])/Times @@ (# - DeleteCases[pts, #]) & /@ pts

Plot[
 {
  Evaluate[ys*Basis[xs, x]],
  Total[ys*Basis[xs, x]]
 }, {x, -9, 7},
 PlotStyle -> {Automatic, Automatic, Automatic, Automatic, Directive[Black, Dashed]},
 Frame -> True,
 Axes -> False,
 ImageSize -> Large,
 GridLines -> Automatic,
 GridLinesStyle -> Directive[LightGray, Dashed],
 PlotStyle -> Thick,
 Epilog -> Table[{PointSize[0.02], ColorData[97][i], Point[pts[[i]]]}, {i,4}],
 PlotLegends -> Placed[
  LineLegend[{
   Subscript[\[ScriptL], 1][x], 
   Subscript[\[ScriptL], 2][x],
   Subscript[\[ScriptL], 3][x], 
   Subscript[\[ScriptL], 4][x],
   L[x]
  }, LegendLayout -> {"Row", 1}],
  Top
 ],
 FrameStyle -> Black
]

许可协议

我，本作品著作权人，特此采用以下许可协议发表本作品：

已授权您依据自由软件基金会发行的无固定段落及封面封底文字（Invariant Sections, Front-Cover Texts, and Back-Cover Texts）的GNU自由文件许可协议1.2版或任意后续版本的条款，复制、传播和／或修改本文件。该协议的副本请见“GNU Free Documentation License”。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	本文件采用知识共享署名-相同方式共享 3.0 未本地化版本许可协议授权。
	您可以自由地： 共享 – 复制、发行并传播本作品 修改 – 改编作品 惟须遵守下列条件： 署名 – 您必须对作品进行署名，提供授权条款的链接，并说明是否对原始内容进行了更改。您可以用任何合理的方式来署名，但不得以任何方式表明许可人认可您或您的使用。 相同方式共享 – 如果您再混合、转换或者基于本作品进行创作，您必须以与原先许可协议相同或相兼容的许可协议分发您贡献的作品。
	本许可协议标签作为GFDL许可协议更新的组成部分被添加至本文件。http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

您可以选择您需要的许可协议。