/*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_derf s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <float.h>
#include <stdio.h>
#include <stdlib.h>
#include "p_derf.h"
 
int main( )
{
	long int _l0;
	double x, xmax, yc, yce, ye;
 
	/*>> 1996-06-17 DRDERF Krogh Minor change in formats for C conversion.
	 *>> 1994-10-19 DRDERF Krogh  Changes to use M77CON
	 *>> 1992-05-13 DRDERF CLL
	 *>> 1991-10-16 DRDERF CLL add demo of DERFCE
	 *>> 1987-12-09 DRDERF   Lawson  Initial Code.
	 *--D replaces "?": DR?ERF, ?ERF, ?ERFC, ?ERFCE
	 *
	 *     Demonstration driver for DERF, DERFC, and DERFCE
	 * */
 
	if (log10( DBL_MIN ) < -65.0e0)
	{
		xmax = 12.0e0;
	}
	else
	{
		xmax = 9.0e0;
	}
	x = -6.0e0;
	printf("    X         DERF            DERFC          DERFCE\n \n");
L_20:
	if (x <= xmax)
	{
		ye = derf( x );
		yc = derfc( x );
		if (x < 0.0e0)
		{
			printf(" %5.1f %17.8e%17.8e\n", x, ye, yc);
		}
		else
		{
			yce = derfce( x );
			printf(" %5.1f %17.8e%17.8e%17.8e\n", x, ye, yc, yce);
		}
		if (x < 6.0e0)
		{
			x += 0.5e0;
		}
		else
		{
			x += 1.0e0;
		}
		goto L_20;
	}
	exit(0);
} /* end of function */