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