/*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_sei s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_sei.h"
 
int main( )
{
	long int j;
	float y, z;
	static float x[14]={-80.e0,-20.e0,-5.e0,-1.e0,-.4e0,-.3e0,-.001e0,
	 .001e0,.3e0,.4e0,1.e0,5.e0,20.e0,80.e0};
		/* OFFSET Vectors w/subscript range: 1 to dimension */
	float *const X = &x[0] - 1;
		/* end of OFFSET VECTORS */
 
	/*>> 1996-05-28 DRSEI    Krogh Added external statement.
	 *>> 1994-10-19 DRSEI    Krogh  Changes to use M77CON
	 *>> 1992-03-16 DRSEI    CLL
	 *>> 1990-11-29 CLL
	 *>> 1987-12-09 DRSEI    Lawson  Initial Code.
	 *     ------------------------------------------------------------------
	 *--S replaces "?": DR?EI, ?EI, ?E1
	 *     ------------------------------------------------------------------ */
	/*     ------------------------------------------------------------------ */
 
	printf("    X             SEI(X)              SE1(X)\n\n");
 
	for (j = 1; j <= 14; j++)
	{
		y = sei( X[j] );
		z = se1( X[j] );
		printf(" %7.3f     %15.8g     %15.8g     \n", X[j], y, z);
	}
	exit(0);
} /* end of function */