/*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_derfce s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_derfce.h"
/*     DRDERFCE
 *>> 1996-05-28 DRDERFCE Krogh Moved formats up.
 *>> 1994-10-19 DRDERFCE Krogh  Changes to use M77CON
 *>> 1992-04-29 DRDERFCE CAO Replaced '1' in format.
 *>> 1992-01-29 DRDERFCE cao/cll eliminate underflow error message.
 *>> 1992-01-29 DRDERFCE cao/cll better output labels.
 *>> 1991-11-20 DRDERFCE CLL Edited for Fortran 90
 *>> 1987-12-09 DRDERFCE Lawson  Initial Code.
 *     DEMONSTRATION DRIVER FOR  DERFCE
 *     ------------------------------------------------------------------
 *--D replaces "?": DR?ERFCE, ?ERFC, ?ERFCE
 *     ------------------------------------------------------------------ */
 
int main( )
{
	long int j;
	double w, x, y, z;
 
 
	/*     ------------------------------------------------------------------
	 * */
	x = 0.0e0;
	printf("    X        DERFCE       DERFC*(X*X)      DIFF\n \n");
	for (j = 1; j <= 18; j++)
	{
		w = derfc( x )*exp( x*x );
		y = derfce( x );
		z = y - w;
		printf(" %5.1f   %15.8g%15.8g %9.2e\n", x, y, w, z);
		x += 0.5e0;
	}
	for (j = 1; j <= 18; j++)
	{
		y = derfce( x );
		printf(" %5.1f   %15.8g\n", x, y);
		x += 1.0e0;
	}
	exit(0);
} /* end of function */