IMPLICIT REAL*8 (A-H,O-Z)
EXTERNAL FQQQQ
COMMON/BPRINT/IPT,NFILE,NDIG,NPUNCH,JPT,MFILE
COMMON/BNWRP1/IVFLIM,INVLIM
COMMON/BNWRP2/EPCOR,EPSVRT
COMMON/BFIDIF/FDFRAC,FDMIN
COMMON/BSTACK/AINT(1000)
COMMON/BSTAK/NQ,NTOP
DIMENSION X(25),FPD(25),SPD(25,25)
CHARACTER*8 ALABEL(25)
open (unit=8,file='newtst.out',status='unknown')
call dflt
NQ=1000
C You should verify that the solution to the equations specified in FQQQQ is
C X(1)=X(2)=X(3)=X(4)=X(5)=1.0
X(1)=1.5
X(2)=1.5
X(3)=1.5
X(4)=1.5
X(5)=1.5
ACC=1.D-6
FDFRAC=1.D-4
FDMIN=1.D-6
ITERL=50
CALL LABEL(ALABEL,NP)
CALL NEWRAP(X,NP,ACC,ITERL,ALABEL,FQQQQ,IER)
WRITE (6,1234) IER
1234 FORMAT(' IER = ',I3)
STOP
END
SUBROUTINE FQQQQ(X,NP,FPD,*)
IMPLICIT REAL*8 (A-H,O-Z)
DIMENSION X(NP),FPD(NP)
C The following step is needed because X(1), X(2), X(3), X(4) will be
C raised to a power and X(3) will also have the log taken---the RETURN 1
C ensures that such a step will not be accepted
IF(X(1).LE.0..OR.X(2).LE.0..OR.X(3).LE.0..OR.X(4).LE.0.) RETURN 1
TEMP=DEXP(X(2))*DLOG(X(3))+X(4)
IF(TEMP.LE.0.) RETURN 1
FPD(1)=X(1)*TEMP**X(5)-1.
FPD(2)=X(1)+X(2)*X(3)+X(4)*X(5)-3.
FPD(3)=X(1)-X(2)+2.*X(3)**2-X(4)-X(5)
FPD(4)=X(1)**X(2)+X(2)**X(1)-X(4)**X(5)-X(3)**X(4)
FPD(5)=X(1)+X(2)+X(3)+X(4)+X(5)-5.
RETURN
END
**********************************************************************
OUTPUT FILE FOLLOWS
**********************************************************************
NEWTON-RAPHSON BEGINS
IVFLIM = 7
INVLIM = 7
EPCOR = 0.100000E+00
EPSVRT = 0.100000E-14
ACC = 0.100000E-05
ISTOP = 0 0
ITER TYPE N
NO
X-VECTOR
0 NORM 0 0.150000E+01 0.150000E+01 0.150000E+01 0.150000E+01 0.150000E+01
G-VECTOR
0 NORM 0 0.806239E+01 0.300000E+01 0.150000E+01 0.000000E+00 0.250000E+01
X-VECTOR
1 NORM 0 0.750000E+00 0.146451E+01 0.142857E+01 0.626829E+00 0.198009E+01
G-VECTOR
1 NORM 0 0.247634E+01 0.108333E+01 0.760204E+00 0.340340E+00 0.125000E+01
X-VECTOR
2 NORM 0 0.133418E+00 0.218893E+01 0.146254E+01 0.235875E+00 0.160424E+01
G-VECTOR
2 NORM 0 0.550688E-01 0.713215E+00 0.382409E+00-0.700161E-01 0.625000E+00
X-VECTOR
3 NORM 0 0.327121E+00 0.182040E+01 0.125892E+01 0.335220E+00 0.125834E+01
G-VECTOR
3 NORM 0-0.335220E+00 0.406774E-01 0.829207E-01 0.142733E-01-0.442424E-12
X-VECTOR
4 NORM 0 0.838395E+00 0.111671E+01 0.107576E+01 0.921765E+00 0.104736E+01
G-VECTOR
4 NORM 0-0.339845E-01 0.513697E-02 0.670920E-01-0.695599E-01-0.563327E-12
X-VECTOR
5 NORM 0 0.964080E+00 0.103734E+01 0.101571E+01 0.993107E+00 0.989763E+00
G-VECTOR
5 NORM 0-0.528702E-03 0.657190E-03 0.721262E-02-0.100433E-01 0.207723E-12
X-VECTOR
6 NORM 0 0.983931E+00 0.102115E+01 0.100645E+01 0.998534E+00 0.989944E+00
G-VECTOR
6 NORM 0-0.123789E-03 0.151024E-03 0.171734E-03-0.588473E-03 0.399680E-14
X-VECTOR
7 NORM 0 0.991930E+00 0.101059E+01 0.100323E+01 0.999295E+00 0.994953E+00
G-VECTOR
7 NORM 0-0.281314E-04 0.377555E-04 0.207027E-04-0.170505E-03 0.766054E-14
CONVERGENCE REACHED
NORM IS LESS THAN ACC
X VECTOR
X1 9.9193032716E-01
X2 1.0105942720E+00
X3 1.0032278421E+00
X4 9.9929490151E-01
X5 9.9495265725E-01
Return to
|Sect. E|Beginning|