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