34
35
36
37#include "implicit_f.inc"
38
39
40
41 INTEGER IPT
42
44 . h(*), p1(*), p2(*), p3(*)
45
46
47
48
50 . g(3,8)
51 DATA g/-.5773502691896,-.5773502691896,-.5773502691896,
52 1 .5773502691896,-.5773502691896,-.5773502691896,
53 2 .5773502691896, .5773502691896,-.5773502691896,
54 3 -.5773502691896, .5773502691896,-.5773502691896,
55 4 -.5773502691896,-.5773502691896, .5773502691896,
56 5 .5773502691896,-.5773502691896, .5773502691896,
57 6 .5773502691896, .5773502691896, .5773502691896,
58 7 -.5773502691896, .5773502691896, .5773502691896/
59
60 CALL basis8(g(1,ipt),g(2,ipt),g(3,ipt),h,p1,p2,p3)
61 RETURN
subroutine basis8(r, s, t, h, pr, ps, pt)