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