2 | MRS - belka, nb
r[i_] » 11-2 - 4 flml ♦ 6 f, - 4 ft.x ♦ f- p h4
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
t0 - 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 f0 ♦ 6 fj - 4 f2 ♦ fj - 0 0.01 ♦ fo - 4 fr . 6 t II 0.01 ♦ ti - 4 f* ♦ 6f,- 4 f4 ♦ U ** 0
0.01 ♦ U - 4 fj • 6 f4 - 4 f5 ♦ f6 * 0
0.01 + fj - 4 f4 ♦ 6fj- 4 f4 ♦ f, « 0
0.01 ♦ f« - 4 fs ♦ 6 f*- 4 f, + f, » 0
0.01 ♦f5-4fe.6f7-4fj.f,= 0 0.01 ♦ f j - 4 f7 ♦ 6 fj - 4 fj » fio = 0 0.01 ♦ f7 - 4 f* . 6 f»- 4 fł0 • fu 0
U2 ■ Ul /. f.i -♦ fi /. ffi.i -* - fB.i; U2 // Tabl©Form
0.01 - 4 f0 ♦ 7 ft - 4 f2 «• fj = 0
0.01 ♦ f0 - 4 fi ♦ 6 f2 - 4 f3 ♦ f« » 0
0.01 ♦ fj - 4 f2 * 6 f3 - 4 f4 * fs = 0
0.01 ♦ fj- 4 fj ♦ 6 f« - 4 f5 ♦ f«» 0
0.01 ♦ f3- 4 U * 6 fs- 4 ft ♦ f7 » 0
0.01 ♦ f« - 4 f5 ♦ 6 f« - 4 f7 ♦ f* - 0
0.01 ♦ f5 - 4 ft ■» 6 fT-4 f* ♦ f, • 0
0.01*f*-4f?*6f$-4f»*fio- 0
0.01 ♦ f7 - 4 f* - 5 f, - 4 f,0 « 0
wb* {f0 = 0, fn = 0);
UR = Join[wb, U2); UR // TableForm
f0 — 0 fio - 0
0.01 - 4 f0 ♦ 7 f! - 4 fj ♦ f3 « 0
0.01 ♦ f0 - 4 ft * 6 f2 - 4 fj ♦ f« 0
0.01 <• fi - 4 f2 * 6 fj - 4 f« ♦ fj « 0
0.01 ♦ f2 - 4 fj ♦ 6 f« - 4 fs ♦ f» = 0
0.01 ♦ fj - 4 f4 * 6 f5 - 4 f* * f7 » 0
0.01 ♦ fĄ - 4 f5 ♦ 6 f«- 4 f, ♦ f8 = 0
0.01 ♦ f5 - 4 U * 6 f7- 4 fs ♦ f, * 0
0.01»fj-4f7»6f|-4f}»f;
0.01 ♦ f7 - 4 f, ♦ 5 f,- 4 fio « 0