1 (313)

pytania w

™’T^! }m'*^Uy ,kryPt "'⁰*° “^ftWlemć °^B10WO w>^irowad!wno B do których odnoszą |

interpolacja.py

1    import numpy

2

3    x • numpy.linspaceC-l, 1, 6)

4    y *• numpy.exp(x)

5

6    n ■ x.shape[0]

7

8    def rozwin_wielomian(dwumiany):

9    n ■ dwumiany. ahape [0] + 1

10    wielomian ■ numpy. zer os ((n))

11    wielomian [0] ■ dwumiany [0]

12    wielomian[1] * 1.0

13    for i in rangę(1 ,n-l):

14    nowy ■ numpy.zeros((n))

15    for | In rangę(i+1) :

16    nowy[j+lQ| +* wielomian [j]

17    nowytj] +- wielomian [j] 1 dwumiany [i]

18    wielomian * nowyC:]

19    return wielomian

20

21    wielomian * numpy .zerosC(n))

22

23 for i in rangę (n) :

mianownik ■ 1 for k in rangę (n) : ifk !« i:

mianownik li Cx[i]

licznik = [J for k in range (n) : ifk I- i:

licznik.append(-x[k]) tablica = numpy.array(licznik)

licznik = rozwin.wielomian(tablica) L - licznik / mianownik * y[i] wielomian +* L

24

25

x[k])

26

27

28

29

30

31

32

33

34

35

36

37    napis * ^,,P^,/,d(x) = '/,.4f" l (n, wielomian[iJ)

38    for i in rangę (1, n) :

39    napis +* " + 5(.41 * sTM¹'' l CwielomianCi], i)

40

41    print napis