3462403626

Jacek Złydach (JW2)

Metody Numeryczne - Zadanie Dodatkowe III

dist(n,w) dist(n,w) 2dist(n,w)

V

- 1

^A2n,(2n)+1 ^A(2n)+l,2n ^A(2n)+l,(2n)+l ^A(2n)+1 ,(2n)+3

return A

Poniższa funkcja tworzy fragment macierzy współczynników odpowiedzialny za ostatni wielomian składów' SPLINE-2 i ostatni węzeł interpolacyjny. Przyjęty jest wariant z parametryzowaną pochodną wostatnim węźle.

Rost(n,w) :=

2n,2n

- dist(n,w)

<- dist(n,w)

(2n+l),(2n)+l return A

Poniższa funkcja tworzy fragment macierz wyrazów wolnych B odpowiedzialny za n-tą parę równań ciągłości.

B(n,rmax, w,y) :=

2n,0 *n+l ^n

^Brmax-I ₀    ® Zapewnienie właściwego rozmiaru macierzy potrzebnego do

sumowania.

^B(2n)+1,0    ⁰

return B

Poniższa funkcja oblicza współczynniki krzywej SPLINE-2 dla podanych węzłów w i wartości y w tych węzłach oraz wartości ostatniej pochodnej d przez bezpośrednie rozwiązanie równania macierzowego.

last(w)-2

S_mat(w,y,d) :=

^    R(i, 2 last(w), w) + Rost(last(w) - 1, w)

last(w)-l

'y'    B(i,21ast(w),w,y)

Parametryzacja pochodnej

r <— A B return r

str. 5/19