Как заставить эту программу посчитать число Лангранжа и вывести его?

Question

[KetrinP]

import numpy as np
import matplotlib.pyplot as plt

x = np.array([-6, 7, 4, 3], dtype=float)
y = np.array([-297, 249, 33, 9], dtype=float)


def lagranz(x, y, t):
    z = 0
    for j in range(len(y)):
        p1 = 1;
        p2 = 1
        for i in range(len(x)):
            if i == j:
                p1 = p1 * 1;
                p2 = p2 * 1
            else:
                p1 = p1 * (t - x[i])
                p2 = p2 * (x[j] - x[i])
        z = z + y[j] * p1 / p2
    return z


xnew = np.linspace(np.min(x), np.max(x), 100)
ynew = [lagranz(x, y, i) for i in xnew]
plt.plot(x, y, 'o', xnew, ynew)
plt.grid(True)
plt.show()

Comments

[KetrinP]

если бы она поддавалась этим пинкам. А как это прописать?
Dr. Bacon,

[KetrinP]

Dr. Bacon, График то она спокойно выводит, а нужно еще число

[Developer]

KetrinP, график это и есть числа. Вот их и выводи

[Алан Гибизов]

Начни с правильного оформления кода, как принято на этом ресурсе.
Ты просишь помощи, но ты просишь без уважения.

[KetrinP]

Алан Гибизов, Я бы с радостью начала с правильного оформления. Но я не знаю как его делать. И понимаю, что это ужасно

[Алан Гибизов]

KetrinP, наверху окна редактора есть кнопочки, помогающие это сделать. Смелее тыкайте в них и пробуйте.

Answer 1

[Алан Гибизов]

Может, что-то такое:

spoiler

import numpy as np
import matplotlib.pyplot as plt

x = np.array([-6, 7, 4, 3], dtype=float)
y = np.array([-297, 249, 33, 9], dtype=float)


def lagranz(x, y, t):
    z = 0
    for j in range(len(y)):
        p1 = 1;
        p2 = 1
        for i in range(len(x)):
            if i == j:
                p1 = p1 * 1;
                p2 = p2 * 1
            else:
                p1 = p1 * (t - x[i])
                p2 = p2 * (x[j] - x[i])
        z = z + y[j] * p1 / p2
    return z


xnew = np.linspace(np.min(x), np.max(x), 100)
ynew = [lagranz(x, y, i) for i in xnew]

counter = 0
step = 9
for xx, yy, tex in zip(xnew, ynew, xnew):
    if counter % step == 0:
        t = plt.text(xx, yy, round(tex, 2), horizontalalignment='right' if xx < 0 else 'left', 
                 verticalalignment='center', fontdict={'color':'red' if xx < 0 else 'green', 'size':14}) 
    counter += 1

plt.plot(x, y, 'o', xnew, ynew)
plt.grid(True)
plt.show()

(Содрано в инфернете, кривовато допилено в части шага отображения).