Подпрограмма CALL DBCPOL. 
Вычисление целевой функции fc2

Входные параметры (аргументы) подпрограммы: nperUC — количество оптимизируемых переменных; WRK2 — массив размерности nperUC, содержащий в WRK2(1) зарядное напряжение на конденсаторе и в WRK2(2) его ёмкость.

Выходной параметр fc2- вычисленное по входным параметрам значение целевой функции fc2. Целевую функцию fc 1 вычисляет пп DBCPOL из математической библиотеки IMSL (при вызове этой подпрограммы в качестве выходного параметра целевая функция fcl обозначена ValueBCPOL). Алгоритм, используемый этой подпрограммой, описан в п. 2.8.

Ниже приводится текст пп CALL_DBCPOL.

SUBROUTINE CALL_DBCPOL (nperUC, WRK2, ValueBCPOL).

use MSIMSL; use DESCRIBE; use DESCRIBE2; use OPTUC; IMPLICIT NONE.

INTEGER:maxfcn, nperUC.

REAL (8): WRK2(nperUC), ftol, fc2, fcl, tsr, ValueBCPOL EXTERNAL CALLJ DM UcO=WRK2(l); Emk=WRK2(2).

nVizov=0; maxfcn=1000; ftol=ld-7; nVizov2=0; fcMin=ld20.

kVizBCPOL=kVizBCPOL+l.

IF (optimizacia) THEN.

100 FORMATf fc2 UcO Emk, mkF Tsr, ms WcO, Dj kViz VdOpt, m/s Wkin, Dj').

IF (kVizBCPOL > kDopViz2) THEN; kVizBCPOL=0; write (*,*)'kDopViz2(DBCPOL)='₍kDopViz2; Enter 0->STOP or 1->C0NTINUE' write (9,*)'kDopViz2(DBCPOL)=', kDopViz2,' Enter 0->STOP or 1->C0NTINUE'.

READ (*,*) i; if (i==0) STOP ENDIF.

CALL DBCPOL (CALL_IDM, npar, AR, 0, ZNIZWRK, ZVRXWRK, ftol, maxfcn, DTX, ValueBCPOL).

ENDIF l (optimizacia).

10 tsr=ValueBCPOL; fc2=ABS (ldO-tsr/tzad); j=0.

DO 1=1,7; IF (OPTPAR (i)) THEN; j=j+l; OPTXAR (i)=DTX (j);ENDIF; ENDDO.

IF (fc2.

WRITE (9," (' n fcl nwK twKr twKz rKatvn &.

& rDnar rDvn hDisk x v rKatNar')").

CALL CALL_IDM (npar, DTX, fcl).

WRITE (8,110)fc2,UcO, Emk, fcl, WcO, maxfcn+l, Vd_fclOpt, massPodv*Vd_fclOpt**2/2dO.

IVdjTlOpt-скор.диска в момент сраб. ИДМ.

WRITE (9,100);

WRITE (9,HO)fc2,UcO, Emk, fcl, WcO, maxfcn+l, Vd_fclOpt, massPodv*Vd_fclOpt**2/2dO lmaxfcn-кол.обращ.к CALLJDM из DBCPOL. +1 -т.к.после CALLJDM есть ещё 1 вызов WRITE!*, 100).

WRITE!*, H0) fc2, UcO, Emk, fcl, WcO, maxfcn+l, Vd_fclOpt, massPodv*Vd_fclOpt**2/2dO 110 FORMAT (3pF8.3,0pF8.1,6pF8.1,3pF8.B, 0pF8.3,i6,F9.2,F10.2).

END SUBROUTINE CALL_DBCPOL.