/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:31:59 */
/*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 "sppnml.h"
#include <stdlib.h>
		/* PARAMETER translations */
#define	SQRT2	1.414213562373095048801688724209698078569672e0
		/* end of PARAMETER translations */
 
float /*FUNCTION*/ sppnml(
float u,
float mu,
float sigma)
{
	float sppnml_v;
 
	/* Copyright (c) 1996 California Institute of Technology, Pasadena, CA.
	 * ALL RIGHTS RESERVED.
	 * Based on Government Sponsored Research NAS7-03001.
	 *>> 1998-10-29 SPPNML Krogh  Moved external statement up for mangle.
	 *>> 1994-10-20 SPPNML Krogh  Changes to use M77CON
	 *>> 1994-08-05 SPPNML 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).
	 * */
	/*--S replaces "?": ?PPNML, ?ERFCI */
 
	sppnml_v = mu - sigma*SQRT2*serfci( u + u );
	return( sppnml_v );
} /* end of function */