diff -r -b -B -c1 ansi\examples\xrlft3.c ..\ansi\examples\xrlft3.c
*** ansi\examples\xrlft3.c	Sat Dec 22 12:45:30 2001
--- ..\ansi\examples\xrlft3.c	Sat Dec 01 18:04:14 2001
***************
*** 39,41 ****
  	long idum=(-3);
! 	unsigned long err,i,j,k,l,nn1=NX,nn2=NY,nn3=NZ;
  	float fnorm,***data1,***data2,**speq1;
--- 40,42 ----
  	long idum=(-3);
! 	unsigned long err,i,j,k,nn1=NX,nn2=NY,nn3=NZ;
  	float fnorm,***data1,***data2,**speq1;
diff -r -b -B -c1 ansi\recipes\amotry.c ..\ansi\recipes\amotry.c
*** ansi\recipes\amotry.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\amotry.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\arcmak.c ..\ansi\recipes\arcmak.c
*** ansi\recipes\arcmak.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\arcmak.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\banbks.c ..\ansi\recipes\banbks.c
*** ansi\recipes\banbks.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\banbks.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define SWAP(a,b) {dum=(a);(a)=(b);(b)=dum;}
--- 1,2 ----
! 
  #define SWAP(a,b) {dum=(a);(a)=(b);(b)=dum;}
diff -r -b -B -c1 ansi\recipes\banmul.c ..\ansi\recipes\banmul.c
*** ansi\recipes\banmul.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\banmul.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\bcuint.c ..\ansi\recipes\bcuint.c
*** ansi\recipes\bcuint.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\bcuint.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\beschb.c ..\ansi\recipes\beschb.c
*** ansi\recipes\beschb.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\beschb.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NUSE1 5
--- 1,2 ----
! 
  #define NUSE1 5
diff -r -b -B -c1 ansi\recipes\bessj0.c ..\ansi\recipes\bessj0.c
*** ansi\recipes\bessj0.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\bessj0.c	Sat Dec 01 18:03:02 2001
***************
*** 22,24 ****
  			+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 -b -B -c1 ansi\recipes\broydn.c ..\ansi\recipes\broydn.c
*** ansi\recipes\broydn.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\broydn.c	Sat Dec 01 18:03:02 2001
***************
*** 113,118 ****
  			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;
--- 114,119 ----
  			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;
***************
*** 124,129 ****
  		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);
--- 125,126 ----
diff -r -b -B -c1 ansi\recipes\caldat.c ..\ansi\recipes\caldat.c
*** ansi\recipes\caldat.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\caldat.c	Sat Dec 01 18:03:02 2001
***************
*** 8,10 ****
  	if (julian >= IGREG) {
! 		jalpha=(long)(((float) (julian-1867216)-0.25)/36524.25);
  		ja=julian+1+jalpha-(long) (0.25*jalpha);
--- 9,11 ----
  	if (julian >= IGREG) {
! 		jalpha=(long)(((double) (julian-1867216)-0.25)/36524.25);
  		ja=julian+1+jalpha-(long) (0.25*jalpha);
***************
*** 15,17 ****
  	jb=ja+1524;
! 	jc=(long)(6680.0+((float) (jb-2439870)-122.1)/365.25);
  	jd=(long)(365*jc+(0.25*jc));
--- 16,18 ----
  	jb=ja+1524;
! 	jc=(long)(6680.0+((double) (jb-2439870)-122.1)/365.25);
  	jd=(long)(365*jc+(0.25*jc));
***************
*** 24,26 ****
  	if (*iyyy <= 0) --(*iyyy);
! 	if (julian < 0) iyyy -= 100*(1-julian/36525);
  }
--- 25,27 ----
  	if (*iyyy <= 0) --(*iyyy);
! 	if (julian < 0) *iyyy -= 100*(1-julian/36525);
  }
diff -r -b -B -c1 ansi\recipes\chebpc.c ..\ansi\recipes\chebpc.c
*** ansi\recipes\chebpc.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\chebpc.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\convlv.c ..\ansi\recipes\convlv.c
*** ansi\recipes\convlv.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\convlv.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\correl.c ..\ansi\recipes\correl.c
*** ansi\recipes\correl.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\correl.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\covsrt.c ..\ansi\recipes\covsrt.c
*** ansi\recipes\covsrt.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\covsrt.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define SWAP(a,b) {swap=(a);(a)=(b);(b)=swap;}
--- 1,2 ----
! 
  #define SWAP(a,b) {swap=(a);(a)=(b);(b)=swap;}
diff -r -b -B -c1 ansi\recipes\cyclic.c ..\ansi\recipes\cyclic.c
*** ansi\recipes\cyclic.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\cyclic.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\daub4.c ..\ansi\recipes\daub4.c
*** ansi\recipes\daub4.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\daub4.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\df1dim.c ..\ansi\recipes\df1dim.c
*** ansi\recipes\df1dim.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\df1dim.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\dftcor.c ..\ansi\recipes\dftcor.c
*** ansi\recipes\dftcor.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\dftcor.c	Sat Dec 01 18:03:02 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) {
--- 12,14 ----
  	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 -b -B -c1 ansi\recipes\dlinmin.c ..\ansi\recipes\dlinmin.c
*** ansi\recipes\dlinmin.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\dlinmin.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\dsvbksb.c ..\ansi\recipes\dsvbksb.c
*** ansi\recipes\dsvbksb.c	Sat Dec 22 12:44:44 2001
--- ..\ansi\recipes\dsvbksb.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\f1dim.c ..\ansi\recipes\f1dim.c
*** ansi\recipes\f1dim.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\f1dim.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\fasper.c ..\ansi\recipes\fasper.c
*** ansi\recipes\fasper.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\fasper.c	Sat Dec 01 18:03:02 2001
***************
*** 24,25 ****
--- 25,27 ----
  	avevar(y,n,&ave,&var);
+ 	if (var == 0.0) nrerror("zero variance in fasper");
  	xmin=x[1];
diff -r -b -B -c1 ansi\recipes\fmin.c ..\ansi\recipes\fmin.c
*** ansi\recipes\fmin.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\fmin.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\four1.c ..\ansi\recipes\four1.c
*** ansi\recipes\four1.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\four1.c	Sat Dec 01 18:03:02 2001
***************
*** 16,18 ****
  		}
! 		m=n >> 1;
  		while (m >= 2 && j > m) {
--- 17,19 ----
  		}
! 		m=nn;
  		while (m >= 2 && j > m) {
diff -r -b -B -c1 ansi\recipes\fred2.c ..\ansi\recipes\fred2.c
*** ansi\recipes\fred2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\fred2.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\gaussj.c ..\ansi\recipes\gaussj.c
*** ansi\recipes\gaussj.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\gaussj.c	Sat Dec 01 18:03:02 2001
***************
*** 26,28 ****
  						}
! 					} else if (ipiv[k] > 1) nrerror("gaussj: Singular Matrix-1");
  				}
--- 27,29 ----
  						}
! 					}
  				}
***************
*** 35,37 ****
  		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 -b -B -c1 ansi\recipes\hqr.c ..\ansi\recipes\hqr.c
*** ansi\recipes\hqr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\hqr.c	Sat Dec 01 18:03:02 2001
***************
*** 21,23 ****
  				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 -b -B -c1 ansi\recipes\hufmak.c ..\ansi\recipes\hufmak.c
*** ansi\recipes\hufmak.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\hufmak.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\iindexx.c ..\ansi\recipes\iindexx.c
*** ansi\recipes\iindexx.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\iindexx.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
***************
*** 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 -b -B -c1 ansi\recipes\indexx.c ..\ansi\recipes\indexx.c
*** ansi\recipes\indexx.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\indexx.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\irbit2.c ..\ansi\recipes\irbit2.c
*** ansi\recipes\irbit2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\irbit2.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define IB1 1
--- 1,2 ----
! 
  #define IB1 1
diff -r -b -B -c1 ansi\recipes\lfit.c ..\ansi\recipes\lfit.c
*** ansi\recipes\lfit.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\lfit.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\linmin.c ..\ansi\recipes\linmin.c
*** ansi\recipes\linmin.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\linmin.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\ludcmp.c ..\ansi\recipes\ludcmp.c
*** ansi\recipes\ludcmp.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\ludcmp.c	Sat Dec 01 18:03:02 2001
***************
*** 3,5 ****
  #include "nrutil.h"
! #define TINY 1.0e-20;
  
--- 4,6 ----
  #include "nrutil.h"
! #define TINY 1.0e-20
  
diff -r -b -B -c1 ansi\recipes\medfit.c ..\ansi\recipes\medfit.c
*** ansi\recipes\medfit.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\medfit.c	Sat Dec 01 18:03:02 2001
***************
*** 30,31 ****
--- 31,33 ----
  	f1=rofunc(b1);
+ 	if (sigb > 0.0) {
  		b2=bb+SIGN(3.0*sigb,f1);
***************
*** 56,57 ****
--- 58,60 ----
  				b2=bb;
+ 			}
  		}
diff -r -b -B -c1 ansi\recipes\mgfas.c ..\ansi\recipes\mgfas.c
*** ansi\recipes\mgfas.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mgfas.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mglin.c ..\ansi\recipes\mglin.c
*** ansi\recipes\mglin.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mglin.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\midinf.c ..\ansi\recipes\midinf.c
*** ansi\recipes\midinf.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\midinf.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define FUNC(x) ((*funk)(1.0/(x))/((x)*(x)))
--- 1,2 ----
! 
  #define FUNC(x) ((*funk)(1.0/(x))/((x)*(x)))
diff -r -b -B -c1 ansi\recipes\midpnt.c ..\ansi\recipes\midpnt.c
*** ansi\recipes\midpnt.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\midpnt.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define FUNC(x) ((*func)(x))
--- 1,2 ----
! 
  #define FUNC(x) ((*func)(x))
diff -r -b -B -c1 ansi\recipes\mmid.c ..\ansi\recipes\mmid.c
*** ansi\recipes\mmid.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mmid.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mp2dfr.c ..\ansi\recipes\mp2dfr.c
*** ansi\recipes\mp2dfr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mp2dfr.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define IAZ 48
--- 1,2 ----
! 
  #define IAZ 48
diff -r -b -B -c1 ansi\recipes\mpdiv.c ..\ansi\recipes\mpdiv.c
*** ansi\recipes\mpdiv.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mpdiv.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
***************
*** 22,24 ****
  	mpmul(rr,s,u,n+MACC,n);
! 	mpsad(s,rr,n+n+MACC/2,1);
  	mpmov(q,&s[2],n-m+1);
--- 22,24 ----
  	mpmul(rr,s,u,n+MACC,n);
! 	mpsad(s,rr,n+MACC-1,1);
  	mpmov(q,&s[2],n-m+1);
diff -r -b -B -c1 ansi\recipes\mpinv.c ..\ansi\recipes\mpinv.c
*** ansi\recipes\mpinv.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mpinv.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mpmul.c ..\ansi\recipes\mpmul.c
*** ansi\recipes\mpmul.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mpmul.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mpops.c ..\ansi\recipes\mpops.c
*** ansi\recipes\mpops.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mpops.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define LOBYTE(x) ((unsigned char) ((x) & 0xff))
--- 1,2 ----
! 
  #define LOBYTE(x) ((unsigned char) ((x) & 0xff))
diff -r -b -B -c1 ansi\recipes\mprove.c ..\ansi\recipes\mprove.c
*** ansi\recipes\mprove.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mprove.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mrqcof.c ..\ansi\recipes\mrqcof.c
*** ansi\recipes\mrqcof.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mrqcof.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\mrqmin.c ..\ansi\recipes\mrqmin.c
*** ansi\recipes\mrqmin.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\mrqmin.c	Sat Dec 01 18:03:02 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\orthog.c ..\ansi\recipes\orthog.c
*** ansi\recipes\orthog.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\orthog.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\pade.c ..\ansi\recipes\pade.c
*** ansi\recipes\pade.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\pade.c	Sat Dec 01 18:03:04 2001
***************
*** 37,41 ****
  	} 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;
--- 38,42 ----
  	} 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;
***************
*** 44,46 ****
  		cof[j]=y[j];
! 		cof[j+n] = -x[j];
  	}
--- 45,47 ----
  		cof[j]=y[j];
! 		cof[j+n] = -z[j];
  	}
diff -r -b -B -c1 ansi\recipes\period.c ..\ansi\recipes\period.c
*** ansi\recipes\period.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\period.c	Sat Dec 01 18:03:04 2001
***************
*** 21,22 ****
--- 22,24 ----
  	avevar(y,n,&ave,&var);
+ 	if (var == 0.0) nrerror("zero variance in period");
  	xmax=xmin=x[1];
diff -r -b -B -c1 ansi\recipes\polcoe.c ..\ansi\recipes\polcoe.c
*** ansi\recipes\polcoe.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\polcoe.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\polin2.c ..\ansi\recipes\polin2.c
*** ansi\recipes\polin2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\polin2.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\predic.c ..\ansi\recipes\predic.c
*** ansi\recipes\predic.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\predic.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\psdes.c ..\ansi\recipes\psdes.c
*** ansi\recipes\psdes.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\psdes.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NITER 4
--- 1,2 ----
! 
  #define NITER 4
diff -r -b -B -c1 ansi\recipes\pwt.c ..\ansi\recipes\pwt.c
*** ansi\recipes\pwt.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\pwt.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\pzextr.c ..\ansi\recipes\pzextr.c
*** ansi\recipes\pzextr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\pzextr.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\qsimp.c ..\ansi\recipes\qsimp.c
*** ansi\recipes\qsimp.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\qsimp.c	Sat Dec 01 18:03:04 2001
***************
*** 9,13 ****
  	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 -b -B -c1 ansi\recipes\qtrap.c ..\ansi\recipes\qtrap.c
*** ansi\recipes\qtrap.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\qtrap.c	Sat Dec 01 18:03:04 2001
***************
*** 9,13 ****
  	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 -b -B -c1 ansi\recipes\quadvl.c ..\ansi\recipes\quadvl.c
*** ansi\recipes\quadvl.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\quadvl.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\ran0.c ..\ansi\recipes\ran0.c
*** ansi\recipes\ran0.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\ran0.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define IA 16807
--- 1,2 ----
! 
  #define IA 16807
diff -r -b -B -c1 ansi\recipes\ran1.c ..\ansi\recipes\ran1.c
*** ansi\recipes\ran1.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\ran1.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define IA 16807
--- 1,2 ----
! 
  #define IA 16807
diff -r -b -B -c1 ansi\recipes\ran2.c ..\ansi\recipes\ran2.c
*** ansi\recipes\ran2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\ran2.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define IM1 2147483563
--- 1,2 ----
! 
  #define IM1 2147483563
diff -r -b -B -c1 ansi\recipes\rj.c ..\ansi\recipes\rj.c
*** ansi\recipes\rj.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\rj.c	Sat Dec 01 18:03:04 2001
***************
*** 22,25 ****
  
! 	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");
--- 23,26 ----
  
! 	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");
***************
*** 61,64 ****
  		delp=(ave-pt)/ave;
! 	} while (FMAX(FMAX(fabs(delx),fabs(dely)),
! 		FMAX(fabs(delz),fabs(delp))) > ERRTOL);
  	ea=delx*(dely+delz)+dely*delz;
--- 62,65 ----
  		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 -b -B -c1 ansi\recipes\rk4.c ..\ansi\recipes\rk4.c
*** ansi\recipes\rk4.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\rk4.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\rkck.c ..\ansi\recipes\rkck.c
*** ansi\recipes\rkck.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\rkck.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\rkdumb.c ..\ansi\recipes\rkdumb.c
*** ansi\recipes\rkdumb.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\rkdumb.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\rzextr.c ..\ansi\recipes\rzextr.c
*** ansi\recipes\rzextr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\rzextr.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\select.c ..\ansi\recipes\select.c
*** ansi\recipes\select.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\select.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp;
--- 1,2 ----
! 
  #define SWAP(a,b) temp=(a);(a)=(b);(b)=temp;
diff -r -b -B -c1 ansi\recipes\selip.c ..\ansi\recipes\selip.c
*** ansi\recipes\selip.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\selip.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\shoot.c ..\ansi\recipes\shoot.c
*** ansi\recipes\shoot.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\shoot.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\shootf.c ..\ansi\recipes\shootf.c
*** ansi\recipes\shootf.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\shootf.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\simp2.c ..\ansi\recipes\simp2.c
*** ansi\recipes\simp2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\simp2.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define EPS 1.0e-6
--- 1,2 ----
! 
  #define EPS 1.0e-6
diff -r -b -B -c1 ansi\recipes\simplx.c ..\ansi\recipes\simplx.c
*** ansi\recipes\simplx.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\simplx.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\simpr.c ..\ansi\recipes\simpr.c
*** ansi\recipes\simpr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\simpr.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\sobseq.c ..\ansi\recipes\sobseq.c
*** ansi\recipes\sobseq.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\sobseq.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\sort.c ..\ansi\recipes\sort.c
*** ansi\recipes\sort.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\sort.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\sort2.c ..\ansi\recipes\sort2.c
*** ansi\recipes\sort2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\sort2.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\sort3.c ..\ansi\recipes\sort3.c
*** ansi\recipes\sort3.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\sort3.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\splin2.c ..\ansi\recipes\splin2.c
*** ansi\recipes\splin2.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\splin2.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\spline.c ..\ansi\recipes\spline.c
*** ansi\recipes\spline.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\spline.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\spread.c ..\ansi\recipes\spread.c
*** ansi\recipes\spread.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\spread.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\stoerm.c ..\ansi\recipes\stoerm.c
*** ansi\recipes\stoerm.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\stoerm.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\svbksb.c ..\ansi\recipes\svbksb.c
*** ansi\recipes\svbksb.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\svbksb.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\svdfit.c ..\ansi\recipes\svdfit.c
*** ansi\recipes\svdfit.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\svdfit.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\svdvar.c ..\ansi\recipes\svdvar.c
*** ansi\recipes\svdvar.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\svdvar.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\toeplz.c ..\ansi\recipes\toeplz.c
*** ansi\recipes\toeplz.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\toeplz.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
***************
*** 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 -b -B -c1 ansi\recipes\trapzd.c ..\ansi\recipes\trapzd.c
*** ansi\recipes\trapzd.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\trapzd.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define FUNC(x) ((*func)(x))
--- 1,2 ----
! 
  #define FUNC(x) ((*func)(x))
diff -r -b -B -c1 ansi\recipes\tridag.c ..\ansi\recipes\tridag.c
*** ansi\recipes\tridag.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\tridag.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\trnspt.c ..\ansi\recipes\trnspt.c
*** ansi\recipes\trnspt.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\trnspt.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\vander.c ..\ansi\recipes\vander.c
*** ansi\recipes\vander.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\vander.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\voltra.c ..\ansi\recipes\voltra.c
*** ansi\recipes\voltra.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\voltra.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\wtn.c ..\ansi\recipes\wtn.c
*** ansi\recipes\wtn.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\wtn.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI
diff -r -b -B -c1 ansi\recipes\zrhqr.c ..\ansi\recipes\zrhqr.c
*** ansi\recipes\zrhqr.c	Sat Dec 22 12:44:46 2001
--- ..\ansi\recipes\zrhqr.c	Sat Dec 01 18:03:04 2001
***************
*** 1,2 ****
! /* note #undef's at end of file */
  #define NRANSI
--- 1,2 ----
! 
  #define NRANSI