1      program DRDERF
2c>> 1996-06-17 DRDERF Krogh Minor change in formats for C conversion.
3c>> 1994-10-19 DRDERF Krogh  Changes to use M77CON
4c>> 1992-05-13 DRDERF CLL
5c>> 1991-10-16 DRDERF CLL add demo of DERFCE
6c>> 1987-12-09 DRDERF   Lawson  Initial Code.
7c--D replaces "?": DR?ERF, ?ERF, ?ERFC, ?ERFCE
8c
9c     Demonstration driver for DERF, DERFC, and DERFCE
10c
11      external         D1MACH, DERF, DERFC, DERFCE
12      double precision D1MACH, DERF, DERFC, DERFCE, X, XMAX, YE, YC, YCE
13c
14      if(log10(D1MACH(1)) .lt. -65.0d0) then
15         XMAX = 12.0d0
16      else
17         XMAX = 9.0d0
18      endif
19      X= -6.0D0
20      print '(4x,''X'',9x,''DERF'',12x,''DERFC'',10x,''DERFCE''/1x)'
21   20 if(X .le. XMAX)then
22         YE= DERF(X)
23         YC= DERFC(X)
24      if(X .lt. 0.0D0)then
25            print '(1x,f5.1,1x,2e17.8)', X,YE,YC
26      else
27            YCE = DERFCE(X)
28            print '(1x,f5.1,1x,3e17.8)', X,YE,YC,YCE
29      end if
30      if(X .lt. 6.0D0)then
31            X = X + 0.5D0
32      else
33            X = X + 1.0D0
34      end if
35      go to 20
36      end if
37      stop
38      end
39