global shapeParameterMu = 1;
global shapeParameterSigma  = .5;
shapeParameterSigma:>0.0001;

sqrt2pi = Sqrt(8*Arctan(1));
category categoryPlaceholder = (rateClassCount, EQUAL, MEAN, 
							 Exp((_x_-shapeParameterMu)(shapeParameterMu-_x_)/(2*shapeParameterSigma*shapeParameterSigma))/(sqrt2pi__*shapeParameterSigma)/ZCDF(shapeParameterMu/shapeParameterSigma), /* density */
							 (1-ZCDF((shapeParameterMu-_x_)/shapeParameterSigma))/ZCDF(shapeParameterMu/shapeParameterSigma), /*CDF*/
							  0, 				   /*left bound*/
							 1e25, 			   /*right bound*/
							 (shapeParameterMu*(1-ZCDF(-shapeParameterMu/shapeParameterSigma)-ZCDF((shapeParameterMu-_x_)/shapeParameterSigma))+
			    			 shapeParameterSigma*(Exp(-shapeParameterMu*shapeParameterMu/(2*shapeParameterSigma*shapeParameterSigma))-
			    			 Exp((_x_-shapeParameterMu)(shapeParameterMu-_x_)/(2*shapeParameterSigma*shapeParameterSigma)))/sqrt2pi__)/ZCDF(shapeParameterMu/shapeParameterSigma)
		  	 				);
