()
.
(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. 
