/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:33:08 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv pf=,p_dcdnml s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_dcdnml.h"
 
int main( )
{
	double c, mu, p, sigma, x;
 
	/*>> 2001-05-25 DRDCDNML Krogh Added comma to format.
	 *>> 1996-05-28 DRDCDNML Krogh  Changes to use M77CON
	 *>> 1994-07-06 DRDCDNML WV Snyder JPL set up for CHGTYP
	 *>> 1994-04-12 DRDCDNML WV Snyder JPL repair format to display sign
	 *
	 *     Evaluate the Cumulative Normal Distribution using DCDNML.
	 *
	 *--D replaces "?": DR?NML, DR?CDNML, ?CDNML, ?PPNML */
 
	x = -4.0e0;
	mu = 0.0e0;
	sigma = 1.0e0;
	printf("        X       C=DCDNML(X)     DPPNML(C)\n");
L_30:
	if (x <= 4.0e0)
	{
		c = dcdnml( x, mu, sigma );
		p = dppnml( c, mu, sigma );
		printf(" %14.7g%14.7g  %14.7g\n", x, c, p);
		x += 0.5e0;
		goto L_30;
	}
	exit(0);
} /* end of function */