Мэйпл 2023 здесь работает хорошо.
Код:
CurveFitting:-PolynomialInterpolation([seq([-Pi + Pi/7*j, sin(-Pi + Pi/7*j)], j = 0 .. 14)], z, form = Lagrange):
evalf(%);
5.180448965*10^(-6)*(z + 3.141592654)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.00006067650209*(z + 3.141592654)*(z + 2.692793703)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.0003026494422*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.0008322859662*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.001334883046*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.001111206303*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.001111206303*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.001334883046*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.0008322859662*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.795195802)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) - 0.0003026494422*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 2.243994753)*(z - 2.692793703)*(z - 3.141592654) + 0.00006067650209*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.692793703)*(z - 3.141592654) - 5.180448965*10^(-6)*(z + 3.141592654)*(z + 2.692793703)*(z + 2.243994753)*(z + 1.795195802)*(z + 1.346396852)*(z + 0.8975979011)*(z + 0.4487989507)*z*(z - 0.4487989507)*(z - 0.8975979011)*(z - 1.346396852)*(z - 1.795195802)*(z - 2.243994753)*(z - 3.141592654)