IR=0
* * (CRIT1=TEN - NSIG)
SQREPS=SQRT (SEPS)
DO I=1, 30 N
IC=1
XI=X (I)
5 AXI=ABS (XI)
IF (i. EQ. 1) GO TO 15
NM1=I - 1
DO 10 J=1, NM1
IF (ABS (XI - X (J)). LT EPS2) XI=XI + ETA
10 the CONTINUE
15 looking=F (XI)
AFXI=ABS (looking)
C TEST FOR CONVERGENCE
IF (EPS) AFXI. LE. GO TO 25
DI=SQREPS
IF (AXI. GE. (P1) DI=SQREPS * AXI
HI=DMIN1 (AFXI, DI)
! HI=1.0 e-06
FXIPHI=F (XI) + HI
DER=(FXIPHI - looking)/HI
ADER=ABS (DER)
IF (ADER. GE. ONE) GO TO 16
IF (AFXI. GE. SINF * ADER) GO TO 20
16 XIPI=looking/DER
XI=XI - XIPI
C TEST FOR CONVERGENCE
ERR1=ABS (XIPI)/DMAX1 (P1, AXI)
IF (ERR1. LE. CRIT1) GO TO 25
IC=IC + 1
IF (IC) LE) ITMAX) GO TO 5
C ROOT NOT FOUND, NO CONVERGENCE
X (I)=111111. D0
IR=IR + 1
IER=33
GO TO 30
C ROOT NOT FOUND, DERIVATIVE=0.
20 X (I)=222222. D0
IR=IR + 1
IER=34
GO TO 30
25 X (I)=XI
30 the CONTINUE
ITMAX=IC
IF (IER. EQ. 0) GO TO 9005
IF (IR. LE. 1) GO TO 9000
IER=35
9000 the CONTINUE
CALL UERTST (IER, 6 hzreal2)
The RETURN of 9005
END
CodePudding user response:
A great god genuflect is begged for helpCodePudding user response:
Fortran, matlab, to Delphi version ~ ~ ~CodePudding user response:
Download this? This stuff works?CodePudding user response:
@ early play big kind of nuclear war