lit
    JINRLIB

TRAPER -

D108

: K.S.Koelbig :

F(x) aaa co a Yi xi (i=1,2,...,N). Yi a c c caap oo epsiloni ( apa ap).

popaa TRAPER po co a paa

d108_1

a oa Yi , co op pa.

popaa c a caapoo oo

d108_2

a F(A) F(B) c c oo o po.

:

: - SUBROUTINE
: - TRAPER

:

CALL TRAPER(X,Y,E,N,A,B,RE,SD), :

X,Y,E - (REAL*8) acc o N, cooc Xi, Yi, epsiloni;
N - (INTEGER) co a ;
A,B - (REAL*8) p ppoa;
RE - (REAL*8) po a paa;
SD - (REAL*8) a caapoo oo.

Ec a epsiloni aac, acc E o o .

Opa:

Xo opa a p A B ( acoc, B o A), o cc , o o oa X(1) o X(N). o ca a F(A) F(B) oa o cpao Y(1) Y(2) Y(N-1) Y(N) cooco, o o pc pa paa. A=B N<2, RE SD 0. poaa, c Xi pa, o pop oo poo.

aa:

TRAPER cooaa oo ocao aa. o o aa ppoa o cooa popa DGAUSS(D103).


:

B d108_3, F(x) 5 .
       . . .
       IMPLICIT REAL*8 (A-H,O-Z)
       DIMENSION X(5),Y(5),E(5)
       DATA X/0.0D0, 1.5708D0, 3.1416D0, 4.7124D0, 6.2832D0/
       DATA Y/0.0D0, 1.0D0,    0.0D0,    1.0D0,    0.0D0/
       DATA E/0.1D0, 0.15D0,   0.2D0,    0.12D0,   0.26D0/
       CALL TRAPER(X,Y,E,5,0.0D0,6.2832D0,RE,SD)
       . . .
:
       RE=   3.141600    SD=   .4874544


home up e-mail