()
.
(1.1)
, , .
, , , ( ) . ( ), , . , , .
.
() . .
(. 1.1).
. 1.1. .
: (); ();
(1.2)
- ( );
(1.3)
;
, , ; (1.4)
; ; . (1.5)
- :
(1.6)
;
(1.7)
;
(1.8)
.
. .
- ( ),
. (1.9)
:
. (1.10)
, (1.1) - :
(1.11)
(1.1). . , :
(1.12)
, (1.13)
; ; ,
(1.14)
, ; , ,
, ; , ; (1.15)
; , ; . (1.16)
, .
(1.13) (1.1) - . , .
(1.1) , , , .
.
, . :
|
|
,
. (1.17)
..
(1.18)
.
,
.. (1.18) .
.
.. (1.19)
.
, , , .
. :
.. . (1.20)
, . .
1.4. .
(1.11) ,
(1.21)
,
(1.22)
.
, ,
; ; ,
, , , , , .
-
function krz_odu2
a=input(' a=');
b=input(' b=');
al1=input(' alfa1='); bt1=input(' beta1=');
gm1=input(' gamma1=');
al2=input(' alfa2='); bt2=input(' beta2=');
gm2=input(' gamma2=');
n=input(' n=');
l=b-a; h=l/(n-1); xi=(a:h:b)'; x=xi(2:n-1);
P=3*x;
q=2+abs(x);
f=2-x;
tm=1/h^2-P/(2*h);
t0=q-2/h^2;
tp=1/h^2+P/(2*h);
A=diag([tm;-bt2/h],-1)+diag([al1-bt1/h;t0;al2+bt2/h])+diag([bt1/h;tp],1);
disp(' '),disp(A)
b=[gm1;f;gm2];
disp(' '),disp(b)
y=A\b;
disp(':'),disp(y)
plot(xi,y),grid on,title('reshenie krajevoj zadachi y=y(x)')
n=9
a=-3
b=3
alfa1=1
beta1=2
gamma1=0
alfa2=1
beta2=0
gamma2=2
n=9
-1.67 2.67 0 0 0 0 0 0 0
6.28 0.69 -2.72 0 0 0 0 0 0
0 4.78 -0.06 -1.22 0 0 0 0 0
0 0 3.28 -0.81 0.28 0 0 0 0
0 0 0 1.78 -1.56 1.78 0 0 0
0 0 0 0 0.28 -0.81 3.28 0 0
0 0 0 0 0 -1.22 -0.06 4.78 0
0 0 0 0 0 0 -2.72 0.69 6.28
0 0 0 0 0 0 0 0 1.00
4.25
3.5
2.75
1.25
0.5
-0.25
:
-1.0629
-0.66432
-4.1819
-5.2704
43.962
44.862
7.6812
11.67
>>
n=17
a=-3
b=3
alfa1=1
beta1=2
gamma1=0
alfa2=1
beta2=0
gamma2=2
n=17
..
:
0.58
0.47
0.32
0.09
-0.14
0.81
5.26
13.34
20.52
22.21
17.38
10.81
6.50
4.44
3.32
2.56
2.00
>>
n=31
|
|
a=-3
b=3
alfa1=1
beta1=2
gamma1=0
alfa2=1
beta2=0
gamma2=2
n=31
..
:
0.57
0.51
0.44
0.36
0.27
0.16
0.07
0.06
0.30
1.06
2.65
5.33
9.06
13.41
17.60
20.76
22.34
21.95
19.83
16.65
13.24
10.21
7.85
6.15
4.96
4.12
3.50
3.01
2.62
2.28
2.00
>>
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.