2 | MRS - belka, nb

r[i_] » 11-2 - 4 f_lml ♦ 6 f, - 4 f_t._x ♦ f- p h⁴

0.01 ♦ f.2.1 - 4 f-i.i ♦ 6 fi - 4 fi»i ♦ £2.1

Na lewo i na prawo musimy dodać dodatkowe (fikcyjne) węzły, ale ich przemieszczenia wyznaczamy z warunków brzegowych

t₀ - 0 f.i » fi f„ = 0

f n« 1 ^>ł * £n-l

f o “ 0

f.i = fi fio = o fu “-f*

Budujemy cały układ równań

Ul = Table [r[i] = 0, {i, 1. n - 1} ] ; Ul / / TableForm

0.01 ♦ f.i - 4 f₀ ♦ 6 fj - 4 f₂ ♦ fj - 0 0.01 ♦ fo - 4 fr . 6 t II 0.01 ♦    ti - 4    f* ♦ 6f,- 4    f₄ ♦    U **    0

0.01 ♦    U - 4    fj • 6 f₄ -    4    f₅ ♦    f₆ *    0

0.01 +    fj - 4    f₄ ♦ 6fj- 4    f₄ ♦    f, « 0

0.01 ♦    f« - 4    f_s ♦ 6 f*- 4    f, +    f, » 0

0.01 ♦f₅-4f_e.6f₇-4f_j.f,= 0 0.01 ♦ f j - 4 f₇ ♦ 6 fj - 4 fj » fio = 0 0.01 ♦ f₇ - 4 f* . 6 f»- 4 f_ł0 • f_u 0

U2 ■ Ul /. f.i -♦ fi /. ffi.i -* - f_B.i; U2 // Tabl©Form

0.01 -    4 f₀ ♦    7 f_t - 4 f₂    «•    fj =    0

0.01 ♦    f₀ - 4    fi ♦ 6 f₂ -    4    f₃ ♦    f« »    0

0.01 ♦    fj - 4    f₂ * 6 f₃ -    4    f₄ *    f_s =    0

0.01 ♦    fj- 4    fj ♦ 6 f« -    4    f₅ ♦    f«» 0

0.01 ♦    f₃- 4    U * 6 f_s-    4    f_t ♦    f₇ » 0

0.01 ♦    f« - 4    f₅ ♦ 6 f« -    4    f₇ ♦    f* - 0

0.01 ♦    f₅ - 4    f_t ■» 6 f_T-4    f* ♦    f, • 0

0.01*f*-4f?*6f$-4f»*fio- 0

0.01 ♦ f₇ - 4 f* - 5 f, - 4 f,₀ « 0

wb* {f₀ = 0, f_n = 0);

UR = Join[wb, U2); UR // TableForm

f₀ — 0 fio - 0

0.01 -    4 f₀ ♦    7 f! - 4 fj    ♦    f₃ «    0

0.01 ♦    f₀ - 4    f_t * 6 f₂    -    4    fj ♦    f«    0

0.01 <•    fi - 4    f₂ * 6 fj    -    4    f« ♦    fj «    0

0.01 ♦    f₂ - 4    fj ♦ 6 f«    -    4    f_s ♦    f» =    0

0.01 ♦    fj - 4    f₄ * 6 f₅    -    4    f* *    f₇ »    0

0.01 ♦    f_Ą - 4    f₅ ♦ 6 f«- 4    f, ♦    f₈ = 0

0.01 ♦    f₅ - 4    U * 6 f₇- 4    f_s ♦    f, * 0

0.01»fj-4f₇»6f|-4f}»f;

0.01 ♦ f₇ - 4 f, ♦ 5 f,- 4 fio « 0