diff -r -c1 recipes_c-kr/doc/ADDRESS ../recipes_c-kr/doc/ADDRESS
*** recipes_c-kr/doc/ADDRESS	Wed Oct 21 01:00:18 1992
--- ../recipes_c-kr/doc/ADDRESS	Fri Nov 30 17:54:58 2001
***************
*** 1,3 ****
       Our address is: Numerical Recipes Software
!                      P.O. Box 243
                       Cambridge, MA 02238 (USA)
--- 1,3 ----
       Our address is: Numerical Recipes Software
!                      P.O. Box 380243
                       Cambridge, MA 02238 (USA)
diff -r -c1 recipes_c-kr/doc/DISTRIBUTION ../recipes_c-kr/doc/DISTRIBUTION
*** recipes_c-kr/doc/DISTRIBUTION	Thu Nov 19 03:55:02 1992
--- ../recipes_c-kr/doc/DISTRIBUTION	Fri Nov 30 18:01:22 2001
***************
*** 10,22 ****
       in your program executable.  Contact Numerical Recipes Software at
!      P.O. Box 243, Cambridge, MA 02238 (USA) [fax 617-863-1739] for
       details.  In distributing a program containing Numerical Recipes
       routines, you acknowledge acceptance of the DISCLAIMER OF WARRANTY
!      below, and of the fact that no business relationship is created
!      between your program's users and Numerical Recipes Software, the
!      authors of the Numerical Recipes books, or Cambridge University Press.
  
!      If you want to distribute software that has Numerical Recipes in the
!      form of source code or individually callable object modules, then you
!      must contact Numerical Recipes Software for further information.  A
!      fee (per Recipe) is charged, and we normally limit the total number
!      of Recipes distributed to 20.
--- 10,21 ----
       in your program executable.  Contact Numerical Recipes Software at
!      P.O. Box 380243, Cambridge, MA 02238 (USA) [fax 617-863-1739] for
       details.  In distributing a program containing Numerical Recipes
       routines, you acknowledge acceptance of the DISCLAIMER OF WARRANTY
!      distributed herewith, and of the fact that no business relationship
!      is created between your program's users and Numerical Recipes
!      Software, the authors of the Numerical Recipes books, or Cambridge
!      University Press.
  
!      For complete information on permissions, please go to
!           http://www.nr.com/infotop.html
!      or contact Numerical Recipes Software at the address given above.
diff -r -c1 recipes_c-kr/doc/ORDERING ../recipes_c-kr/doc/ORDERING
*** recipes_c-kr/doc/ORDERING	Wed Oct 21 01:17:34 1992
--- ../recipes_c-kr/doc/ORDERING	Sun Nov 25 19:07:18 2001
***************
*** 2,4 ****
  
!      Numerical Recipes books, and diskettes for IBM/PC-compatible and
       Apple Macintosh computers, are published exclusively by Cambridge
--- 2,4 ----
  
!      Numerical Recipes books, and CDROMs for IBM/PC-compatible and
       Apple Macintosh computers, are published exclusively by Cambridge
***************
*** 14,15 ****
--- 14,19 ----
  
+          e-mail: orders@cup.org
+ 
+          Web: http://www.cambridge.org
+ 
       Outside of North America, orders may be placed with:
***************
*** 20,33 ****
           Cambridge, CB2 2RU
!          U.K.                                         more...
  
!          fax: (U.K. country code+) 223 315052
           telex: 851 817256 CUPCAM G
           
       The following books are available (prices subject to change):
  
!      Numerical Recipes in C, 2nd Edition (hardcover) $49.95
!      Numerical Recipes C Example Book, 2nd Ed. (paper) $29.95
       
-      Numerical Recipes in FORTRAN, 2nd Edition (hardcover) $49.95
-      Numerical Recipes FORTRAN Example Book, 2nd Ed. (paper) $29.95
       
--- 24,48 ----
           Cambridge, CB2 2RU
!          U.K.
  
!          fax: (U.K. country code+) 223 325891
           telex: 851 817256 CUPCAM G
  
+          e-mail: directcustserv@cambridge.org
+ 
+          Web: http://www.cambridge.org
+          
       The following books are available (prices subject to change):
  
!      Numerical Recipes in Fortran 77, 2nd Edition (hardcover) $64.95
!         (Volume 1 of Fortran Numerical Recipes)
!      Numerical Recipes in Fortran 90  (hardcover) $47.95
!         (Volume 2 of Fortran Numerical Recipes)
!      Numerical Recipes FORTRAN Example Book, 2nd Ed. (paper) $32.95
!      
!      Numerical Recipes in C, 2nd Edition (hardcover) $64.95
!      Numerical Recipes C Example Book, 2nd Ed. (paper) $32.95
! 
!      Numerical Recipes in C++, 2nd Edition (hardcover) $69.95
!      Numerical Recipes Example Book (C++) (paper) $34.95
  
       
***************
*** 35,43 ****
  
!      Numerical Recipes in Pascal, Revised Edition (hardcover) $49.95
       Numerical Recipes Example Book (Pascal), Revised Ed. (paper) $29.95
  
!      Numerical Recipes Routines and Examples in BASIC (paper) $32.50
  
!      Diskettes are available in C, FORTRAN, Pascal, and BASIC (Microsoft
!      QuickBASIC), for IBM/PC-compatibles and Apple Macintosh.  Contact
       Cambridge University Press for pricing information.  A TrueBasic
--- 50,58 ----
  
!      Numerical Recipes in Pascal, Revised Edition (hardcover) $59.95
       Numerical Recipes Example Book (Pascal), Revised Ed. (paper) $29.95
  
!      Numerical Recipes Routines and Examples in BASIC (paper) $37.95
  
!      CDROMs are available for all languages for Windows, Apple Macintosh,
!      and Unix/Linux operating systems. Contact
       Cambridge University Press for pricing information.  A TrueBasic
diff -r -c1 recipes_c-kr/doc/VERSION ../recipes_c-kr/doc/VERSION
*** recipes_c-kr/doc/VERSION	Sat Jul 19 07:26:12 1997
--- ../recipes_c-kr/doc/VERSION	Fri Nov 30 18:46:40 2001
***************
*** 1,2 ****
  Numerical Recipes in C
! Unix Version 2.08
--- 1,2 ----
  Numerical Recipes in C
! UNIX/Linux Version 2.10
diff -r -c1 recipes_c-kr/recipes/banbks.c ../recipes_c-kr/recipes/banbks.c
*** recipes_c-kr/recipes/banbks.c	Mon Jul 14 18:34:06 1997
--- ../recipes_c-kr/recipes/banbks.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define SWAP(a,b) {dum=(a);(a)=(b);(b)=dum;}
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/beschb.c ../recipes_c-kr/recipes/beschb.c
*** recipes_c-kr/recipes/beschb.c	Mon Jul 14 18:34:06 1997
--- ../recipes_c-kr/recipes/beschb.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define NUSE1 5
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/bessj0.c ../recipes_c-kr/recipes/bessj0.c
*** recipes_c-kr/recipes/bessj0.c	Mon Jul 14 18:34:06 1997
--- ../recipes_c-kr/recipes/bessj0.c	Sun Nov 25 15:12:24 2001
***************
*** 23,25 ****
  			+y*(-0.6911147651e-5+y*(0.7621095161e-6
! 			-y*0.934935152e-7)));
  		ans=sqrt(0.636619772/ax)*(cos(xx)*ans1-z*sin(xx)*ans2);
--- 23,25 ----
  			+y*(-0.6911147651e-5+y*(0.7621095161e-6
! 			-y*0.934945152e-7)));
  		ans=sqrt(0.636619772/ax)*(cos(xx)*ans1-z*sin(xx)*ans2);
diff -r -c1 recipes_c-kr/recipes/broydn.c ../recipes_c-kr/recipes/broydn.c
*** recipes_c-kr/recipes/broydn.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/broydn.c	Sun Nov 25 15:12:24 2001
***************
*** 108,113 ****
  			for (sum=0.0,j=1;j<=n;j++) sum += qt[i][j]*fvec[j];
! 			g[i]=sum;
  		}
  		for (i=n;i>=1;i--) {
! 			for (sum=0.0,j=1;j<=i;j++) sum += r[j][i]*g[j];
  			g[i]=sum;
--- 108,113 ----
  			for (sum=0.0,j=1;j<=n;j++) sum += qt[i][j]*fvec[j];
! 			p[i] = -sum;
  		}
  		for (i=n;i>=1;i--) {
! 			for (sum=0.0,j=1;j<=i;j++) sum -= r[j][i]*p[j];
  			g[i]=sum;
***************
*** 119,124 ****
  		fold=f;
- 		for (i=1;i<=n;i++) {
- 			for (sum=0.0,j=1;j<=n;j++) sum += qt[i][j]*fvec[j];
- 			p[i] = -sum;
- 		}
  		rsolv(r,n,d,p);
--- 119,120 ----
diff -r -c1 recipes_c-kr/recipes/caldat.c ../recipes_c-kr/recipes/caldat.c
*** recipes_c-kr/recipes/caldat.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/caldat.c	Sun Nov 25 15:12:24 2001
***************
*** 10,12 ****
  	if (julian >= IGREG) {
! 		jalpha=(long)(((float) (julian-1867216)-0.25)/36524.25);
  		ja=julian+1+jalpha-(long) (0.25*jalpha);
--- 10,12 ----
  	if (julian >= IGREG) {
! 		jalpha=(long)(((double) (julian-1867216)-0.25)/36524.25);
  		ja=julian+1+jalpha-(long) (0.25*jalpha);
***************
*** 17,19 ****
  	jb=ja+1524;
! 	jc=(long)(6680.0+((float) (jb-2439870)-122.1)/365.25);
  	jd=(long)(365*jc+(0.25*jc));
--- 17,19 ----
  	jb=ja+1524;
! 	jc=(long)(6680.0+((double) (jb-2439870)-122.1)/365.25);
  	jd=(long)(365*jc+(0.25*jc));
***************
*** 26,28 ****
  	if (*iyyy <= 0) --(*iyyy);
! 	if (julian < 0) iyyy -= 100*(1-julian/36525);
  }
--- 26,28 ----
  	if (*iyyy <= 0) --(*iyyy);
! 	if (julian < 0) *iyyy -= 100*(1-julian/36525);
  }
diff -r -c1 recipes_c-kr/recipes/covsrt.c ../recipes_c-kr/recipes/covsrt.c
*** recipes_c-kr/recipes/covsrt.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/covsrt.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define SWAP(a,b) {swap=(a);(a)=(b);(b)=swap;}
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/dftcor.c ../recipes_c-kr/recipes/dftcor.c
*** recipes_c-kr/recipes/dftcor.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/dftcor.c	Sun Nov 25 15:12:24 2001
***************
*** 11,13 ****
  	th=w*delta;
! 	if (a >= b || th < 0.0e0 || th > 3.01416e0) nrerror("bad arguments to dftcor");
  	if (fabs(th) < 5.0e-2) {
--- 11,13 ----
  	th=w*delta;
! 	if (a >= b || th < 0.0e0 || th > 3.1416e0) nrerror("bad arguments to dftcor");
  	if (fabs(th) < 5.0e-2) {
diff -r -c1 recipes_c-kr/recipes/fasper.c ../recipes_c-kr/recipes/fasper.c
*** recipes_c-kr/recipes/fasper.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/fasper.c	Sun Nov 25 15:12:24 2001
***************
*** 21,22 ****
--- 21,23 ----
  	avevar(y,n,&ave,&var);
+ 	if (var == 0.0) nrerror("zero variance in fasper");
  	xmin=x[1];
diff -r -c1 recipes_c-kr/recipes/four1.c ../recipes_c-kr/recipes/four1.c
*** recipes_c-kr/recipes/four1.c	Mon Jul 14 18:34:08 1997
--- ../recipes_c-kr/recipes/four1.c	Sun Nov 25 15:12:24 2001
***************
*** 19,21 ****
  		}
! 		m=n >> 1;
  		while (m >= 2 && j > m) {
--- 19,21 ----
  		}
! 		m=nn;
  		while (m >= 2 && j > m) {
diff -r -c1 recipes_c-kr/recipes/gaussj.c ../recipes_c-kr/recipes/gaussj.c
*** recipes_c-kr/recipes/gaussj.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/gaussj.c	Sun Nov 25 15:12:24 2001
***************
*** 27,29 ****
  						}
! 					} else if (ipiv[k] > 1) nrerror("gaussj: Singular Matrix-1");
  				}
--- 27,29 ----
  						}
! 					}
  				}
***************
*** 36,38 ****
  		indxc[i]=icol;
! 		if (a[icol][icol] == 0.0) nrerror("gaussj: Singular Matrix-2");
  		pivinv=1.0/a[icol][icol];
--- 36,38 ----
  		indxc[i]=icol;
! 		if (a[icol][icol] == 0.0) nrerror("gaussj: Singular Matrix");
  		pivinv=1.0/a[icol][icol];
diff -r -c1 recipes_c-kr/recipes/hqr.c ../recipes_c-kr/recipes/hqr.c
*** recipes_c-kr/recipes/hqr.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/hqr.c	Sun Nov 25 15:12:24 2001
***************
*** 22,24 ****
  				if (s == 0.0) s=anorm;
! 				if ((float)(fabs(a[l][l-1]) + s) == s) break;
  			}
--- 22,27 ----
  				if (s == 0.0) s=anorm;
! 				if ((float)(fabs(a[l][l-1]) + s) == s) {
! 					a[l][l-1]=0.0;
! 					break;
! 				}
  			}
diff -r -c1 recipes_c-kr/recipes/iindexx.c ../recipes_c-kr/recipes/iindexx.c
*** recipes_c-kr/recipes/iindexx.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/iindexx.c	Sun Nov 25 15:12:24 2001
***************
*** 20,22 ****
  				a=arr[indxt];
! 				for (i=j-1;i>=1;i--) {
  					if (arr[indx[i]] <= a) break;
--- 20,22 ----
  				a=arr[indxt];
! 				for (i=j-1;i>=l;i--) {
  					if (arr[indx[i]] <= a) break;
***************
*** 32,36 ****
  			SWAP(indx[k],indx[l+1]);
- 			if (arr[indx[l+1]] > arr[indx[ir]]) {
- 				SWAP(indx[l+1],indx[ir])
- 			}
  			if (arr[indx[l]] > arr[indx[ir]]) {
--- 32,33 ----
***************
*** 38,41 ****
  			}
! 			if (arr[indx[l+1]] > arr[indx[l]]) {
! 				SWAP(indx[l+1],indx[l])
  			}
--- 35,41 ----
  			}
! 			if (arr[indx[l+1]] > arr[indx[ir]]) {
! 				SWAP(indx[l+1],indx[ir])
! 			}
! 			if (arr[indx[l]] > arr[indx[l+1]]) {
! 				SWAP(indx[l],indx[l+1])
  			}
***************
*** 43,45 ****
  			j=ir;
! 			indxt=indx[l];
  			a=arr[indxt];
--- 43,45 ----
  			j=ir;
! 			indxt=indx[l+1];
  			a=arr[indxt];
***************
*** 51,53 ****
  			}
! 			indx[l]=indx[j];
  			indx[j]=indxt;
--- 51,53 ----
  			}
! 			indx[l+1]=indx[j];
  			indx[j]=indxt;
diff -r -c1 recipes_c-kr/recipes/irbit2.c ../recipes_c-kr/recipes/irbit2.c
*** recipes_c-kr/recipes/irbit2.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/irbit2.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define IB1 1
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/ludcmp.c ../recipes_c-kr/recipes/ludcmp.c
*** recipes_c-kr/recipes/ludcmp.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/ludcmp.c	Sun Nov 25 15:12:24 2001
***************
*** 2,4 ****
  #include "nrutil.h"
! #define TINY 1.0e-20;
  
--- 2,4 ----
  #include "nrutil.h"
! #define TINY 1.0e-20
  
diff -r -c1 recipes_c-kr/recipes/medfit.c ../recipes_c-kr/recipes/medfit.c
*** recipes_c-kr/recipes/medfit.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/medfit.c	Sun Nov 25 15:12:24 2001
***************
*** 31,32 ****
--- 31,33 ----
  	f1=rofunc(b1);
+ 	if (sigb > 0.0) {
  		b2=bb+SIGN(3.0*sigb,f1);
***************
*** 57,58 ****
--- 58,60 ----
  				b2=bb;
+ 			}
  		}
diff -r -c1 recipes_c-kr/recipes/midinf.c ../recipes_c-kr/recipes/midinf.c
*** recipes_c-kr/recipes/midinf.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/midinf.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define FUNC(x) ((*funk)(1.0/(x))/((x)*(x)))
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/midpnt.c ../recipes_c-kr/recipes/midpnt.c
*** recipes_c-kr/recipes/midpnt.c	Mon Jul 14 18:34:10 1997
--- ../recipes_c-kr/recipes/midpnt.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define FUNC(x) ((*func)(x))
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/mp2dfr.c ../recipes_c-kr/recipes/mp2dfr.c
*** recipes_c-kr/recipes/mp2dfr.c	Mon Jul 14 18:34:12 1997
--- ../recipes_c-kr/recipes/mp2dfr.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define IAZ 48
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/mpdiv.c ../recipes_c-kr/recipes/mpdiv.c
*** recipes_c-kr/recipes/mpdiv.c	Mon Jul 14 18:34:12 1997
--- ../recipes_c-kr/recipes/mpdiv.c	Sun Nov 25 15:12:24 2001
***************
*** 15,17 ****
  	mpmul(rr,s,u,n+MACC,n);
! 	mpsad(s,rr,n+n+MACC/2,1);
  	mpmov(q,&s[2],n-m+1);
--- 15,17 ----
  	mpmul(rr,s,u,n+MACC,n);
! 	mpsad(s,rr,n+MACC-1,1);
  	mpmov(q,&s[2],n-m+1);
diff -r -c1 recipes_c-kr/recipes/mpops.c ../recipes_c-kr/recipes/mpops.c
*** recipes_c-kr/recipes/mpops.c	Mon Jul 14 18:34:12 1997
--- ../recipes_c-kr/recipes/mpops.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define LOBYTE(x) ((unsigned char) ((x) & 0xff))
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/pade.c ../recipes_c-kr/recipes/pade.c
*** recipes_c-kr/recipes/pade.c	Mon Jul 14 18:34:12 1997
--- ../recipes_c-kr/recipes/pade.c	Sun Nov 25 15:12:24 2001
***************
*** 36,40 ****
  	} while (rr < rrold);
! 	*resid=sqrt(rr);
  	for (k=1;k<=n;k++) {
! 		for (sum=cof[k],j=1;j<=k;j++) sum -= x[j]*cof[k-j];
  		y[k]=sum;
--- 36,40 ----
  	} while (rr < rrold);
! 	*resid=sqrt(rrold);
  	for (k=1;k<=n;k++) {
! 		for (sum=cof[k],j=1;j<=k;j++) sum -= z[j]*cof[k-j];
  		y[k]=sum;
***************
*** 43,45 ****
  		cof[j]=y[j];
! 		cof[j+n] = -x[j];
  	}
--- 43,45 ----
  		cof[j]=y[j];
! 		cof[j+n] = -z[j];
  	}
diff -r -c1 recipes_c-kr/recipes/period.c ../recipes_c-kr/recipes/period.c
*** recipes_c-kr/recipes/period.c	Mon Jul 14 18:34:12 1997
--- ../recipes_c-kr/recipes/period.c	Sun Nov 25 15:12:24 2001
***************
*** 21,22 ****
--- 21,23 ----
  	avevar(y,n,&ave,&var);
+ 	if (var == 0.0) nrerror("zero variance in period");
  	xmax=xmin=x[1];
diff -r -c1 recipes_c-kr/recipes/psdes.c ../recipes_c-kr/recipes/psdes.c
*** recipes_c-kr/recipes/psdes.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/psdes.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define NITER 4
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/qsimp.c ../recipes_c-kr/recipes/qsimp.c
*** recipes_c-kr/recipes/qsimp.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/qsimp.c	Sun Nov 25 15:12:24 2001
***************
*** 10,14 ****
  	int j;
! 	float s,st,ost,os;
  
- 	ost = os = -1.0e30;
  	for (j=1;j<=JMAX;j++) {
--- 10,13 ----
  	int j;
! 	float s,st,ost=0.0,os=0.0;
  
  	for (j=1;j<=JMAX;j++) {
diff -r -c1 recipes_c-kr/recipes/qtrap.c ../recipes_c-kr/recipes/qtrap.c
*** recipes_c-kr/recipes/qtrap.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/qtrap.c	Sun Nov 25 15:12:24 2001
***************
*** 10,14 ****
  	int j;
! 	float s,olds;
  
- 	olds = -1.0e30;
  	for (j=1;j<=JMAX;j++) {
--- 10,13 ----
  	int j;
! 	float s,olds=0.0;
  
  	for (j=1;j<=JMAX;j++) {
diff -r -c1 recipes_c-kr/recipes/ran0.c ../recipes_c-kr/recipes/ran0.c
*** recipes_c-kr/recipes/ran0.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/ran0.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define IA 16807
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/ran1.c ../recipes_c-kr/recipes/ran1.c
*** recipes_c-kr/recipes/ran1.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/ran1.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define IA 16807
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/ran2.c ../recipes_c-kr/recipes/ran2.c
*** recipes_c-kr/recipes/ran2.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/ran2.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define IM1 2147483563
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/ran3.c ../recipes_c-kr/recipes/ran3.c
*** recipes_c-kr/recipes/ran3.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/ran3.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- #include <math.h>
  #define MBIG 1000000000
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/rj.c ../recipes_c-kr/recipes/rj.c
*** recipes_c-kr/recipes/rj.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/rj.c	Sun Nov 25 15:12:24 2001
***************
*** 21,24 ****
  
! 	if (FMIN(FMIN(x,y),z) < 0.0 || FMIN(FMIN(x+y,x+z),FMIN(y+z,fabs(p))) < TINY
! 		|| FMAX(FMAX(x,y),FMAX(z,fabs(p))) > BIG)
  			nrerror("invalid arguments in rj");
--- 21,24 ----
  
! 	if (FMIN(FMIN(x,y),z) < 0.0 || FMIN(FMIN(FMIN(x+y,x+z),y+z),fabs(p)) < TINY
! 		|| FMAX(FMAX(FMAX(x,y),z),fabs(p)) > BIG)
  			nrerror("invalid arguments in rj");
***************
*** 60,63 ****
  		delp=(ave-pt)/ave;
! 	} while (FMAX(FMAX(fabs(delx),fabs(dely)),
! 		FMAX(fabs(delz),fabs(delp))) > ERRTOL);
  	ea=delx*(dely+delz)+dely*delz;
--- 60,63 ----
  		delp=(ave-pt)/ave;
! 	} while (FMAX(FMAX(FMAX(fabs(delx),fabs(dely)),
! 		fabs(delz)),fabs(delp)) > ERRTOL);
  	ea=delx*(dely+delz)+dely*delz;
diff -r -c1 recipes_c-kr/recipes/select.c ../recipes_c-kr/recipes/select.c
*** recipes_c-kr/recipes/select.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/select.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp;
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/simp2.c ../recipes_c-kr/recipes/simp2.c
*** recipes_c-kr/recipes/simp2.c	Mon Jul 14 18:34:14 1997
--- ../recipes_c-kr/recipes/simp2.c	Sun Nov 25 15:12:24 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define EPS 1.0e-6
--- 1 ----
diff -r -c1 recipes_c-kr/recipes/toeplz.c ../recipes_c-kr/recipes/toeplz.c
*** recipes_c-kr/recipes/toeplz.c	Mon Jul 14 18:34:16 1997
--- ../recipes_c-kr/recipes/toeplz.c	Sun Nov 25 15:12:26 2001
***************
*** 38,40 ****
  		}
! 		if (sd == 0.0 || sgd == 0.0) nrerror("toeplz-3 singular principal minor");
  		g[m1]=sgn/sgd;
--- 38,40 ----
  		}
! 		if (sgd == 0.0) nrerror("toeplz-3 singular principal minor");
  		g[m1]=sgn/sgd;
diff -r -c1 recipes_c-kr/recipes/trapzd.c ../recipes_c-kr/recipes/trapzd.c
*** recipes_c-kr/recipes/trapzd.c	Mon Jul 14 18:34:16 1997
--- ../recipes_c-kr/recipes/trapzd.c	Sun Nov 25 15:12:26 2001
***************
*** 1,2 ****
- /* note #undef's at end of file */
  #define FUNC(x) ((*func)(x))
--- 1 ----