/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:33:10 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv pf=,p_derfi s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <float.h>
#include <stdio.h>
#include <stdlib.h>
#include "p_derfi.h"
/*     program DRDERFI
 *>> 1996-06-27 DRDERFI  Krogh Minor change in formats for C conversion.
 *>> 1996-05-28 DRDERFI  Krogh Added external statement.
 *>> 1994-10-19 DRDERFI  Krogh  Changes to use M77CON
 *>> 1992-04-21 DRDERFI  CLL
 *>> 1991-11-19 CLL
 *>> 1987-12-09 DRDERFI  Snyder  Initial Code.
 *
 *     Demonstration driver for DERFI and DERFCI.
 *     ------------------------------------------------------------------ */
 
int main( )
{
	long int _l0, i;
	double floor, r, x;
 
	/*     ------------------------------------------------------------------
	 *--D replaces "?": DR?ERFI, ?ERFI, ?ERFCI
	 *     ------------------------------------------------------------------ */
	printf("\n           DRDERFI.. Demo driver for DERFI and DERFCI\n\n                X             DERFI(X) = DERFCI(1-X)\n");
	r = 1.0e-5;
	for (i = 1; i <= 4; i++)
	{
		x = 1.0e0 - r;
		r *= 10.0e0;
		printf("           %12.5g%29.14g\n", x, derfi( x ));
	}
	x = 0.9e0;
	for (i = 1; i <= 9; i++)
	{
		printf("           %12.5g%29.14g\n", x, derfi( x ));
		x -= 0.1e0;
	}
	printf("           %12.5g%29.14g\n", 1.0e-10, derfi( 1.0e-10 ));
	printf("           %12.5g%29.14g\n", 0.0e0, derfi( 0.0e0 ));
 
	printf("\n                X             DERFCI(X) = DERFI(1-X)\n");
	x = 1.0e-10;
	floor = DBL_MIN*1.0e11;
	for (i = 1; i <= 10; i++)
	{
		printf("           %12.5g%29.14g\n", x, derfci( x ));
		if (x < floor)
			goto L_40;
		x *= 1.0e-10;
	}
L_40:
	;
	exit(0);
} /* end of function */