align="justify"> program shag; a, b, h, x, x1, x2: real;, y2: real;: byte; y (x: real): real;: = -3371.9 * x +58413 ;; ('vvedi a, b, h'); (a, b, h);: = 0;: = a + h; x <= b do: = xh;: = x;: = y (x1) ;: = y (x2); y1 * y2 <0 then: = k +1; (x1: 5:2); (x2: 5:2);;: = x + h;; k = 0 then (' korney net ');.
Результат програми:
vvedi a, b, h
.1
.30
.40
У результаті обчислень кроковим методом отримуємо інтервал ізоляції кореня [17,3; 17,4]
Уточнення кореня методом Ньютона:
program nuton;
var a, b, e '); (a, b, e); u (a) <> 0 thenu (a) * z (a)> 0x: = ax: = b;: = 0; abs ( y (x))> e do: = i +1;: = xy (x)/u (x);; ('x =', x: 8:3);.
Результат програми:
vvedi a, b, e
.3
.4
.001
x = 17.323
стаціонарна швидкість Vст = 17,323 м/c
4.3 Знаходження часу розгону судна
Визначаємо час розгону судна - час за який судно виходить за стаціонарну швидкість:
t =? 1/F (V), де F (V) = T (V)-R (V)/m.
t =? 40400/(-3371,9 * x +58413)
Метод центральних прямокутників:
Program metodcentrpr; a, b, h, n, x, S: real;: integer;: array [1 .. 51] of real;; ('vvedite a'); (a) ; ('vvedite b'); (b);: = (ba)/50;: = (a + h)/2;: = 0; i: = 1 to 51 do [i]: = 40400/(- 3371.9 * x +58413);: = x + h;: = S + c [i];;: = h * S; ('S =', S: 6:3);;.
Результат програми:
vvedite a
vvedite b
.323
S = 46.290
Час розгону до Vст t = 46,29 сек.
4.4 Знаходження шляху розгону судна
Обчислюємо шлях розгону, шлях який пройде судно за час розгону.
Вирішуємо систему диференціальних рівнянь:
В
Метод Ейлера:
Program eiler;, xn, h, y0, N, D, z0, kp: real;, y, z, y1, V, z1: ​​real; f (x, z: real): real;: = (-3010.9 * x +56297)/40400;; ('Vvedi x0, xn, y0, z0, V, D'); (x0, xn, y0, z0, V, D);: = x0 ; y: = y0; z: = z0;: = V/D;: = xn/D; x <= xn do; ('t =', x: 8:4, 'S =', y: 8: 4, 'V =', z: 8:4);: = x + h;: = y + h * z;
z: = z + N;;
