IER = 0
IR=0
CRIT1 = TEN**(-NSIG)
SQREPS = SQRT(SEPS)
DO 30 I=1,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 CONTINUE
15 FXI = F(XI)
AFXI = ABS(FXI)
C TEST FOR CONVERGENCE
IF (AFXI .LE. EPS) GO TO 25
DI = SQREPS
IF (AXI .GE. P1) DI = SQREPS*AXI
HI = DMIN1(AFXI,DI)
! HI = 1.0E-06
FXIPHI = F(XI + HI)
DER = (FXIPHI - FXI)/HI
ADER = ABS(DER)
IF (ADER .GE. ONE) GO TO 16
IF (AFXI .GE. SINF*ADER) GO TO 20
16 XIPI=FXI/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 CONTINUE
ITMAX = IC
IF(IER.EQ.0) GO TO 9005
IF(IR.LE.1) GO TO 9000
IER=35
9000 CONTINUE
CALL UERTST(IER,6HZREAL2)
9005 RETURN
END
uj5u.com熱心網友回復：
跪求大神幫忙uj5u.com熱心網友回復：
Fortran轉matlab，發到Delphi版~~~uj5u.com熱心網友回復：
把這個下載下來嗎？這個東西好用嗎？uj5u.com熱心網友回復：
@早打大打打核戰爭轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/31023.html
標籤：語言基礎/算法/系統設計
上一篇：delphi真是讓人抓狂,矛盾的存在，完全不給1兩年經驗的delphi程式員一點作業機會
下一篇：求助Matlab作圖