Список використаної літератури
1. Воробйов Г.Н., Бахвалов Н.С. " Чисельні методи". М .: Наука, 2006. 231 с.
. Єфімов А.В., Демидович Б.П. «Лінійна алгебра та основи математичного аналізу». М .: Наука, 2011. 386 с.
. Бараненко Г.С., Демидович Б.П. «Завдання і вправи з математичного аналізу для втузів». М .: Наука, 2010. 184 с.
. Абрамов С.А., Зима Е.В. «Початок програмування мовою Паскаль». М .: Наука, 2007. 8 с.
. Епанешников А.Є., Красильников Ю.І. «Програмування в середовищі турбо Паскаль». М .: Центр МІФІ СП Діалог, 2008. 3-6 с.
Додаток
Метод хорд: myhord;
interfacecrt; hord=object, b, c, cp, eps: real; init; shag; poisk: real ;; funct (x: real): real; funct (x: real): real ;:=2 * x - 3 * sin (2 * x) - 1; hord.init ;; ( Metod HORD ); ( Vvedite pogreshnost eps=); (eps); ( Vvedite nachalo otrezka: a=); (a); ( Vvedite konets otrezka: b=); (b);:=b;:=b - (funct (b) * (ab))/( funct (a) - funct (b)) ;; hord.shag; (funct (a) * funct (c) gt; 0) a:=cb:=c;:=c;:=b - (funct (b ) * (ab))/(funct (a) - funct (b)) ;; hord.poisk: real; (abs (c-cp) gt; eps) do (funct (c)=0) then ;;; ;:=c ;;.
Метод Зейделя: Dec_sys;
crt;=array [1 .. 10,1..10] of real;=array [1..10] of real; _sys=object: integer; Zeidel (a: matrix; b:vector; x: vector; e: real); Dec_Zeidel; Init; Done ;;
TDec_sys. Zeidel (a: matrix; b: vector; x: vector; e: real); i, j, flag: integer ;, s2, s, v, m: real;
:=1; i:=1 to n do:=0; j:=1 to n doj lt; gt; i then:=s + abs (a [i] [j]); s gt;=abs (a [i] [i]) then:=0 ;;
if (flag=1) then:=0; i:=1 to n do:=0;:=0;
j:=1 to i do:=s1 + (a [i] [j] * x [j]); j:=i + 1 to n do:=s2 + a [i] [j] * x [j];
:=x [i]; [i]:=x [i] - ((1/a [i] [i]) * (s1 + s2 - b [i]));
(abs (vx [i]) gt; m) then:=abs (v - x [i]) ;;; until (m gt;=e);
( ђҐіҐ ЁҐ бЁб⥬л: ); i:=1 to n do ( x raquo ;, i, = raquo ;, x [i]: 5: 5 );
writeln ( Ёб⥬ Ґ б室Ёвбп! ) ;; TDec_sys. Dec_Zeidel; eps: real;: matrix ;, x: vector ;, j: integer ;; ( ђҐіҐ ЁҐ бЁб⥬л га ў Ґ Ё©:raquo;);(laquo;a11ъx1+a12ъx2+a13ъx3=b1raquo;);(laquo;a21ъx1+a22ъx2+a23ъx3=b1raquo;);(laquo;a31ъx1+a32ъx2+a33ъx3=b1raquo;);
; (, ўҐ¤ЁвҐ в®з ®бвм: ); readln (eps);
; (, ўҐ¤ЁвҐ н «Ґ¬Ґ вл а ᳨ॠ® © ¬ ваЁжл бЁб⥬л: ); i:=1 to n doj:=1 to n do ( a raquo ;, i, j, =); (a [i] [j]) ;;;
; (, ўҐ¤ЁвҐ н «Ґ¬Ґ вл бв®« Ўж бў®Ў®¤ ле з «Ґ ®ў: ); i:=1 to n do ( b raquo ;, i, =); (b [i]) ;;
i:=1 to n do [i]:=0 ;;;
(a, b, x, eps);
;
;
TDec_sys. Init;:=3 ;;
TDec_sys. Done ;;;
.
Метод Рунге-Кутта: rk; crt; TEiler=object, y, h, b, dx1, dx2, dx3, dx4, ddx: real; init; shag; poisk ;; proizv (f, g: real): real; proizv (f, g: real): real;:=g * gf * f ;; TEiler.init ;;:=0;:=0.5;:=0.1;:=1 ;; TEiler.shag;:=proizv (x, y);:=proizv (x + h/2, y + h * dx1/2);:=proizv (x + h/2, y + h * dx2/2); :=proizv (x + h, y + h * dx3);:=(dx1 + 2 * dx2 + 2 * dx3 + dx4)/6;:=y + ddx * h;:=x + h ;; TEiler. poisk; f: text; (f, 1.txt ); (f); (x lt;=b) do ( x= raquo ;, x: 3: 3, raquo ;, y = raquo ;, y: 3: 3); (f, x, raquo ;, y) ;;; (f) ;;
end.
