--- ray/src/rt/rayinit.cal	1991/08/22 08:57:19	1.1
+++ ray/src/rt/rayinit.cal	1991/10/09 11:36:45	1.3
@@ -3,10 +3,6 @@
 {
 	Initialization file for Radiance.
 
-	4/14/86
-}
-
-{
 	The following are predefined:
 
 	Dx, Dy, Dz			- ray direction
@@ -134,3 +130,24 @@ turbulenceb(x,y,z,s) = if( s-1.01, 0,
 turbulencec(x,y,z,s) = if( s-1.01, 0,
 			sgn(noise3(x/s,y/s,z/s))*noise3c(x/s,y/s,z/s) +
 			turbulencec(x,y,z,2*s) );
+
+			{ Normal distribution from uniform range (0,1) }
+
+un2`private(t) : t - (2.515517+t*(.802853+t*.010328))/
+		(1+t*(1.432788+t*(.189269+t*.001308))) ;
+un1`private(p) : un2`private(sqrt(log(1/p/p))) ;
+
+unif2norm(p) : if( .5-p, un1`private(p), -un1`private(1-p) ) ;
+
+nrand(x) = unif2norm(rand(x));
+
+			{ Local (u,v) coordinates for planar surfaces }
+crosslen`private = Nx*Nx + Ny*Ny;
+			{ U is distance from origin in XY-plane }
+U = if( crosslen`private - FTINY,
+		(Py*Nx - Px*Ny)/crosslen`private,
+		Px);
+			{ V is defined so that N = U x V }
+V = if( crosslen`private - FTINY,
+		Pz - Nz*(Px*Nx + Py*Ny)/crosslen`private,
+		Py);