/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:33:18 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv pf=,p_spsi s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_spsi.h"
 
int main( )
{
	long int i, ie;
	float e, p, x;
 
	/*>> 1996-06-19 DRSPSI Krogh  Minor change in format for C conversion.
	 *>> 1994-10-19 DRSPSI Krogh  Changes to use M77CON
	 *>> 1994-09-01 DRSPSI CLL
	 *>> 1994-07-21 DRSPSI WVS Original Code
	 *--S replaces "?": DR?PSI, ?PSI, ?PSIE */
	printf("          X        SPSI(X)     FLAG  EST.ABS.ERR\n");
	x = -4.5e0;
	for (i = 1; i <= 14; i++)
	{
		p = spsi( x );
		spsie( &e, &ie );
		printf("  %14.8g%14.8g%4ld  %14.8g\n", x, p, ie, e*fabsf( p ));
		x += 1.0e0;
	}
	exit(0);
} /* end of function */