.
. .
, ( ) , , .
: ; , .
- i .
, , , . , , .
1. : / .. . .: . -.1985. -448.
2. - / . , .. , .. , .. : , 2005. 46.
3. - / . , .. , .. .., : , 2005. 44.
4. - ( ) / . , .. , .. : , 2005. 86.
1. .., .., .., .. . . : . 1999. 250.
2. . Visual Basic 6: . .: . 1999. 576 .
3. .. : . .: . 1987. 240 .
4. : . .: - , 1978. 336 .
5. ., . Turbo Pascal 7.0 .: -. 1995. 288 .
6. ., . ( ). .: . 1974. 831 .
7. . Excel , . .: . 1999. 522 .
8. .. Visual Basic for Applications . .: . 1999. 512 .
9. . Word 2000: : . .: : BHV. 1999. 432 .
10. / .., .., .., .. . 1999. 44 .
|
|
- -
1. / . .. . .: , 1977. 415 .
2. / . .. . .: , 1975. 416 .
3. : / . .. , .: , 1987. 424 .
4. .: , 2000. 511 .
5. .., .., .. , . .: , 1974. 240 .
6. / .., .., .., .. . 1999. 44 .
7. .., .., .. . . .: , 1984. 400 .
8. ó : . . / .. , .. , .. , .. , .. ; ϳ . . .. . 2- . . . : , 2006 669.: .
1. 㳿 : . . . / .. , .. , .. .; . .. 2- , .-.: ., 1999 416 .: .
2. : . 2 . / .. , .. , .. ³ . .: , 1995.- 256 .: .
3. .., .., .. : : . . .-.: . ., 1984.-112 .
4. : / .. , .. , .. . // ϳ . . .. .- .: , 1990.- 400 .
5. .., .. : . . 5- .: , 1990. 448 .
6. .. : .: , 1984. 208 .
7. / ϳ . . ... .: , 1974.- 676 .
8. .., .. / ϳ . .. Գ, - 3- .: , 1988. 137 .
9. : , , / .. , .. , .. ./ ϳ . .. . 3- .: . 1980. 527.
10. / ϳ . .. . .2-. .: . 1975. 708.
11. ' / ϳ . .. . - .: . 1986. 448.
12. .., .. . - 2- .: . 1983. 223.
13. 3.1404-86. . .: - , 1986.- 56.
|
|
14. -. 2 . / ϳ .. .. . 4- - .: . 1985. .1. - 656.; .2. 495 .
15. .., .. . 4- - .: . 1975. 352.
16. .., .. . - .: . 1984. 288.
17. .., .. , . .: , 1988. 480 .
18. .., .. Գ . .: , 1987. 208 .
19. : / ϳ . .. . .: , 1986. 191 .
20. .. . 2- .: ³ ., 1986. 400 .
21. .. : . . .: , 1987. 320 .
22. .. : . .: , 1986. 272 .
23. : / .. , , .. . / ϳ . .. . 2- .: , 1990. 768 .
24. - / .. , .. , .. / ϳ . .. . .: , 1987. 560 .
25. .., .., .. . .: , 1988. 421 .
26. .., .. : . - .: , 1986. 320 .
27. .. : . . - .: , 1979. 288 .
28. : : 2 . / . : .. / ϳ . .. , .. . .: , 1984. . 1 - 582 .; . 2 - 495 .
()
II I
i i
ó
( )
:
:
: -09 ..
: ____________ ____________
: ____________ ____________
: ____ _________ 200__.
: _________________
− 2009
{
- }
program RK4;
uses
crt;
const
N = 2; { }
M = 5100; { }
{ X,Y1,Y2 }
type
Working_arrays = array [1..N] of Real;
Array_X = array [1..M] of Single;
Arrays_of_outcomes = array [1..N,1..M] of Single;
var
h,x,v,xk: Real;
i,j,l: Integer;
y,k,f,w: Working_arrays;
OX: Array_X;
OY: Arrays_of_outcomes;
{ }
procedure system_equations(var y,f:Working_arrays);
begin
f[1]:=y[2];
f[2]:=(y[1]/x-y[2])/x-y[1]
end;
{ . }
procedure Runge_Kutt_4(var i: Integer;
var OX: Array_X;
var OY: Arrays_of_outcomes);
begin
for j:=1 to n do y[j]:=w[j];
i:=0;
Repeat
i:=i+1;
system_equations(y,f);
for j:=1 to n do
begin
v:=h*f[j];
k[j]:=v;
y[j]:=w[j]+v/2
end;
x:=x+h/2;
system_equations(y,f);
for j:=1 to n do
begin
v:=h*f[j];
k[j]:=k[j]+2*v;
|
|
y[j]:=w[j]+v/2
end;
system_equations(y,f);
for j:=1 to n do
begin
v:=h*f[j];
k[j]:=k[j]+2*v;
y[j]:=w[j]+v
end;
x:=x+h/2;
OX[i]:=x;
system_equations(y,f);
for j:=1 to n do
begin
y[j]:=w[j]+(k[j]+h*f[j])/6;
OY[j,i]:=y[j];
w[j]:=y[j]
end
Until ((x > xk) or (i = M))
end;
procedure output(device:string);
var
d: Text;
begin
assign(d,device);
rewrite(d);
writeln(d,' :');
writeln(d);
for l:=1 to i do
begin
write(d,' x=',OX[l]:6:2);
for j:=1 to N do write(d,' y[',j:1,']=',
OY[j,l]:10:6);
writeln(d)
end;
writeln(d);
close(d)
end;
BEGIN { }
{ 1. }
ClrScr;
write(' h=> ');
readln(h);
write(' x=> ');
readln(x);
write(' x=> ');
readln(xk);
writeln(' Y:');
for j:=1 to n do
begin
write(' Y[',j:1,']=> ');
readln(w[j])
end;
{ 2. }
Runge_Kutt_4(i,OX,OY);
{ 3. }
ClrScr;
output('con');
output('result.txt');
write(' ');
repeat until UpCase(ReadKey) <>#0
END.