1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
*DECK DBDIFF
SUBROUTINE DBDIFF (L, V)
C***BEGIN PROLOGUE DBDIFF
C***SUBSIDIARY
C***PURPOSE Subsidiary to DBSKIN
C***LIBRARY SLATEC
C***TYPE DOUBLE PRECISION (BDIFF-S, DBDIFF-D)
C***AUTHOR Amos, D. E., (SNLA)
C***DESCRIPTION
C
C DBDIFF computes the sum of B(L,K)*V(K)*(-1)**K where B(L,K)
C are the binomial coefficients. Truncated sums are computed by
C setting last part of the V vector to zero. On return, the binomial
C sum is in V(L).
C
C***SEE ALSO DBSKIN
C***ROUTINES CALLED (NONE)
C***REVISION HISTORY (YYMMDD)
C 820601 DATE WRITTEN
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C***END PROLOGUE DBDIFF
C
INTEGER I, J, K, L
DOUBLE PRECISION V
DIMENSION V(*)
C***FIRST EXECUTABLE STATEMENT DBDIFF
IF (L.EQ.1) RETURN
DO 20 J=2,L
K = L
DO 10 I=J,L
V(K) = V(K-1) - V(K)
K = K - 1
10 CONTINUE
20 CONTINUE
RETURN
END
|
