/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:31:56 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include "dppnml.h"
#include <stdlib.h>
		/* PARAMETER translations */
#define	SQRT2	1.414213562373095048801688724209698078569672e0
		/* end of PARAMETER translations */
 
double /*FUNCTION*/ dppnml(
double u,
double mu,
double sigma)
{
	double dppnml_v;
 
	/* Copyright (c) 1996 California Institute of Technology, Pasadena, CA.
	 * ALL RIGHTS RESERVED.
	 * Based on Government Sponsored Research NAS7-03001.
	 *>> 1998-10-29 DPPNML Krogh  Moved external statement up for mangle.
	 *>> 1994-10-20 DPPNML Krogh  Changes to use M77CON
	 *>> 1994-08-05 DPPNML WV Snyder
	 *
	 *     Compute the Percentage Point of the Normal (Gaussian) Probability
	 *     Distribution, g**(-1) (u;mu,sigma).  This is given in terms of the
	 *     inverse of the co-error function, erfci, by x = mu -
	 *     sigma * sqrt(2) * erfci(2*u).
	 * */
	/*--D replaces "?": ?PPNML, ?ERFCI */
 
	dppnml_v = mu - sigma*SQRT2*derfci( u + u );
	return( dppnml_v );
} /* end of function */