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