ВУЗ:
Составители:
88
if x<=0 then x:=Random*0.01;
if x>=1 then x:=1-Random*0.01;
CB:=b*CA0*Xa*x;
CC:=c*CA0*Xa*(1-x);
r2:=k2*(exp(nB*ln(CB)));
r4:=k4*(exp(nC*ln(CC)));
if (r1-r2/b)<>0 then
F:=1/(1+r3/(r1-r2/b)-r4/c/(r1-r2/b))
else F:=Fs;
end;
2:begin
if x<=0 then x:=Random*0.01;
if x>=1 then x:=1-Random*0.01;
XCp:=1/c*exp(1/c*ln(Kp2/Kp1))*exp(b/c*ln(b*x));
F:=1-XCp-1/Kp1*exp(b*ln(b*x));
end;
end; {end case}
end; {end function F}
Приложение 7
Решение системы нелинейных уравнений
модифицированным методом Ньютона-Рафсона
1. Подпрограмма решения системы нелинейных уравнений (Qbasic)
SUB Rafson (N, X() AS DOUBLE, F() AS DOUBLE, intX() AS DOUBLE, E)
REM Решение системы нелинейных уравнений
REM модифицированным методом Ньютона-Рафсона
REM Система уравнений вводится в подпрограмме SystemUrawnenie в виде
REM F(1) = f(x1,x2,...xn)=0
REM F(2) = f(x1,x2,...xn)=0
REM ....
REM F(N) = f(x1,x2,...xn)=0
REM Подпрограмма Rafson использует при своей работе
REM подпрограмму Rmatr
REM N - число уравнений
REM X(N) - массив значений X(i) (решение системы уравнений)
REM F(N) - массив значений F(i) (значения функций-уравнений для набора
X(i))
REM intX(N,2) - интервалы варьирования для каждой переменной X(i)
REM E - относительная погрешность
Nitro PDF Trial
www.nitropdf.com
Страницы
- « первая
- ‹ предыдущая
- …
- 86
- 87
- 88
- 89
- 90
- …
- следующая ›
- последняя »