fy"> 40 = X - V = (X - W) * (FX - FV) = P * (FX - FW) = P * Q - (X - W) * R = 2.0 * ( Q - R)
IF (Q. GT.0. D0) P =-P = Abs (Q) = E = D
Continue
IF (Abs (P). GE. Abs (0.50 * Q * R)) GoTo 40
IF (P. LE. Q * (A - X)) GoTo 40
IF (P. GE. Q * (B - X)) GoTo 40 = P/Q = X + D
IF ((U - A). LT. TOLM) D = Sign (TOLM, XH-X)
IF ((B - U). LT. TOLM) D = Sign (TOLM, XH-X)
GoTo 50
Continue
IF (X. GE. XH) E = A - X
IF (X. LT. XH) E = B - X = C * E
Continue = Abs (D)
IF (FU. LT. TOLM) U = X + Sign (TOLM, D)
IF (FU. GE. TOLM) U = X + D = F (U)
IF (FU. GT. FX) GoTo 60
IF (U. GE. X) A = X
IF (U. LT. X) B = X = W = FW = X = FX = U = FU
GoTo 20
Continue
IF (U. LT. X) A = U
IF (U. GE. X) B = U
IF (FU. LE. FW) GoTo 70
IF (W. EQ. ​​X) GoTo 70
IF (FU. LE. FV) GoTo 80
IF (V. EQ. ​​X) GoTo 80
IF (V. EQ. ​​W) GoTo 80
GoTo 20
Continue = W = FW = U = FU
GoTo 20
Continue = U = FU
GoTo 20
Continue = X = FX
RETURN
END Subroutine FMIN
program Console3
Real * 8 AX, BX, F, TOL, XM, FM
External F = -1 = 0 = 0.005
CALL FMIN (AX, BX, F, TOL, XM, FM)
write (*, *) 'MINIMUM NAIDEN NA'
write (*, *) 'X', XM
write (*, *) 'Y', FM
end program Console3
Результат виконання програми:
В
Метод "Гауса" для рішень системи лінійних алгебраїчних рівнянь на FORTRAN
Налагодження і тестування мовою FORTRAN
Вихідний код
SUBROUTINE DSIMQLL (A, B, N, TOL, IER)
DIMENSION A (1), B (1)
INTEGER N, IER, JJ, IQS, NY, J, JY, IT, I, IMAX, IXJ, IXJX
DOUBLE PRECIS...
