Вторая пачка часть 111
§ 1 и 2 в гл. 2). Однако их при¬ме-ня¬ют в прак¬ти¬ческих рас¬че¬тах до¬стବточ¬но час¬то, так как они удобны, прос¬ты, не требуют сп嬬¬ци¬аль¬-ных про¬грамм. И, что осо¬бен¬но важ¬но, позволяют бы¬стро сделать пред¬ва¬ри¬тель¬ные (пр謬¬¬ки¬доч¬¬ные) расчеты.
Процедуры NEW1, NEW2, NEWTON из § 1 и 2 (гл.2) впо묬не пригодны для вычисления про-из¬вод¬ных. Но с уче¬том ре¬шаемой задачи они должны быть не¬сколь¬ко мди¬ф謬ц謬ро¬ва¬ны, хотя формальные па¬ра¬мет¬ры ос¬та¬ют¬ся у прце¬дур без изменений.
{***ПЕРВАЯ ИНТЕРПОЛЯЦИОННАЯ ФОРМУЛА НЬЮТОНА**}
FUNCTION NEW1 (KEY:INTEGER;X,Y,R1,R2,R3,R4: MAS;
X1:REAL):REAL;
LABEL 30;
VAR I,J,K : INTEGER; Q : REAL;
BEGIN
I := 1;
J := 14;
IF X1 < X[I] THEN
GOTO 30;
{ ****МЕТОД ДВОИЧНОГО ПОИСКА ИНТЕРВАЛА****}
REPEAT
K := (I+J) DIV 2;
IF X1< X[K] THEN
J := K;
IF X1>=X[K] THEN
I := K;
UNTIL J <= I+1;
30:
Q := (X1 - X[I]) / (X[2]-X[1]);
{ ****ВЫБОР ТИПА РАБОТЫ ПРОЦЕДУРЫ****}
CASE KEY OF
0:NEW1:=Y[I]+Q*(R1[I]+(Q-1)*(R2[I]/2.0+
(Q-2)*R3[I]/6.0));
1: NEW1 := R1[I] + 0.5*(2*Q-1)*R2[I] +
((3*Q-6)*Q+2)*R3[I]/6+
Окончание таблицы 2.1
№ Вариант А Вариант Б № Вариант А Вариант Б
Параметр i Х Y X Y Параметр i Х Y X Y
1 0.11000 9.00000 0.10000 2.00000 Расчетные 14 0.26708 3.71492 0.37083 5.27809
2 0.12208 8.13717 0.12083 2.60852 точки, 15 0.27917 3.54559 0.39167 5.23644
3 0.13417 7.40129 0.14167 3.08333 через 16 0.29125 3.38354 0.41250 5.12278
Расчетные 4 0.14625 6.77476 0.16250 3.46364 которые 17 0.30333 3.22789 0.43333 4.93141
точки, 5 0.15833 6.24154 0.18333 3.77930 проходит 18 0.31542 3.07873 0.45417 4.66007
через 6 0.17042 5.78712 0.20417 4.05183 интерполя- 19 0.32750 2.93704 0.47500 4.31085
которые 7 0.18250 5.39849 0.22500 4.29535 ционная 20 0.33958 2.804663 0.49583 3.89127
проходит 8 0.19458 5.06405 0.24583 4.51758 кривая 21 0.35167 2.68421 0.51667 3.41516
интерполя- 9 0.20667 4.77360 0.26667 4.72085 22 0.36375 2.57910 0
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа