Оптимизация химико-технологических процессов органического синтеза. Бочкарев В.В - 86 стр.

UptoLike

Составители: 

86
a = 1 + Sign * RND
IF kcicle > 20 THEN
PRINT равнение x = F(x) не имеет решения"
GOTO m2
END IF
END IF
x1 = (a * x0 - xm1) / (a - 1)
F1 = F(x1, NF)
IF (ABS(x1 - F1) > .00001) THEN
IF (ABS(x0 - F0) < ABS(xm1 - Fm1)) THEN
xm1 = x0: Fm1 = F0: x0 = x1: F0 = F1: GOTO m1
ELSE
x0 = x1: F0 = F1: GOTO m1
END IF
END IF
m2:
END SUB
FUNCTION F! (x, NF)
SELECT CASE NF
CASE 1
F = x ^ 2 + x - 9
CASE 2
F = SIN(x) + x - .5
END SELECT
END FUNCTION
2. Процедура нахождения решение уравнениях = F(x) (Turbo Pascal)
PROCEDURE Urawnenie(var x1:real; x:real;NF:integer);
{ Решение уравнения x=f(x) }
{Процедура составлена Бочкаревым В.В., ТПУ, кафедра ТООС}
{ Вид функции f(x) задается функцией F(x,NF) }
{ х - начальное приближение; х1 - решение
NF - номер функции в FUNCTION F(x,NF)
Абсолютная погрешность вычислений - 0.00001 }
var
xm1,x0,Fm1,F0,F1,a:real;
kcicle, Sign :integer;
BEGIN
xm1:=x;
Fm1:=F(xm1,NF);
Nitro PDF Trial
www.nitropdf.com