lit
    JINRLIB

DCAUCH -

D104

: K.S.Koelbig :

popa c pa

d103_1

aao ooc.

:

: - FUNCTION
: - DGAUSS
: - F - /-

:

B o apco pa DGAUSS(F,A,B,EPS) a po a paa I oo ooc.
apap:

F - opopa-, cocao ooa c opao .
F oa ocaa opaopo EXTERNAL popa, a DGAUSS.
A,B - (REAL*8) p p ppoa o X.
EPS - (REAL*8) aaa ooc c (C.pa "ooc").

:

oo paa [A,B] opc G8(A,B) G16(A,B) a 8 16-oo p

d103_1

oo apapo op acca cc a:

d103_2

Toa:

d103_3

, aa c X0 = A oa Xk = B, o
XI (I=1,2,....) pca a: XI = XI-1 + LAMBDA(B-XI-1)
LAMBDA paa po ocoaoc 1, 1/2, 1/4... , oopo:
R(XI-1 ,XI ) < EPS.

Ec a ao-o c pocca o:

d103_4

ac acoo a, o 1.D0+0.005D0*Q ca o o 1.D00 -a ao ooc, o popaa ppa pao, DGAUSS pcaac a 0.D0

:

c cpooo copa oo pa a F(X) paa [A,B], a apa EPS oo ca p p ooco opoc p ABS(I)>1, a p ABS(I)<1 - op p acoo opoc. o oo, c -co opao, p apoca (c.o), , c

d103_5

oa oo:

d103_6

o ca c, DGAUSS apa pao coo o o. F, cpoc a op [A,B], ooc oo , ocaa c.

:

. acoc, B o A.

:

D103.1: pa ooc o oca (c. "o").
oaa pao oo o o.

:

a F o oa paa (A,B) pc. Coao, popa o paoa, oa a op.


:
       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       EXTERNAL F
       . . .
       Y=DGAUSS(F,1.D0,2.D0,1.D-12)
       WRITE(*,*) 'Y=',Y
       . . .
       DOUBLE PRECISION FUNCTION F(X)
       IMPLICIT REAL*8 (A-H,O-Z)
       F=X**3
       RETURN
       END
:
       Y=3.750000000000000


home up e-mail