/*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_serf s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <float.h>
#include <stdio.h>
#include <stdlib.h>
#include "p_serf.h"
 
int main( )
{
	long int _l0;
	float x, xmax, yc, yce, ye;
 
	/*>> 1996-06-17 DRSERF Krogh Minor change in formats for C conversion.
	 *>> 1994-10-19 DRSERF Krogh  Changes to use M77CON
	 *>> 1992-05-13 DRSERF CLL
	 *>> 1991-10-16 DRSERF CLL add demo of SERFCE
	 *>> 1987-12-09 DRSERF   Lawson  Initial Code.
	 *--S replaces "?": DR?ERF, ?ERF, ?ERFC, ?ERFCE
	 *
	 *     Demonstration driver for SERF, SERFC, and SERFCE
	 * */
 
	if (log10f( FLT_MIN ) < -65.0e0)
	{
		xmax = 12.0e0;
	}
	else
	{
		xmax = 9.0e0;
	}
	x = -6.0e0;
	printf("    X         SERF            SERFC          SERFCE\n \n");
L_20:
	if (x <= xmax)
	{
		ye = serf( x );
		yc = serfc( x );
		if (x < 0.0e0)
		{
			printf(" %5.1f %17.8e%17.8e\n", x, ye, yc);
		}
		else
		{
			yce = serfce( 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 */