NUMERICAL RECIPES THIRD EDITION v3.02 to v3.04 FILE COMPARISON
   
Left base folder: NR302\code (previous version files) 
Right base folder: NR304\code (new version files) 
   
File: besselfrac.h  

L4	Doub xo,nuo;
R4	Doub xjy,nujy,xik,nuik,xai,xsph;
L11	Bessel() : xo(9.99e99), nuo(9.99e99), sphno(-9999) {}
R11	Bessel() : xjy(9.99e99), nujy(9.99e99), xik(9.99e99), nuik(9.99e99),
	xai(9.99e99), sphno(-9999) {}
L17	if (nu != nuo || x != xo) besseljy(nu,x);
R18	if (nu != nujy || x != xjy) besseljy(nu,x);
L21	if (nu != nuo || x != xo) besseljy(nu,x);
R22	if (nu != nujy || x != xjy) besseljy(nu,x);
L25	if (nu != nuo || x != xo) besselik(nu,x);
R26	if (nu != nuik || x != xik) besselik(nu,x);
L29	if (nu != nuo || x != xo) besselik(nu,x);
R30	if (nu != nuik || x != xik) besselik(nu,x);
L205	xo = x;
	nuo = nu;
R206	xjy = x;
	nujy = nu;
L326	xo = x;
	nuo = nu;
R327	xik = x;
	nuik = nu;
L359	if (x != xo) airy(x);
R360	if (x != xai) airy(x);
L363	if (x != xo) airy(x);
R364	if (x != xai) airy(x);
L381	if (n != sphno || x != xo) sphbes(n,x);
R382	if (n != sphno || x != xsph) sphbes(n,x);
L385	if (n != sphno || x != xo) sphbes(n,x);
R386	if (n != sphno || x != xsph) sphbes(n,x);

   
File: deviates.h  

L218	if (k < 0) continue;
R218	if (k < 0 || k > n) continue;

   
File: hmm.h  

L79	while (lhood < BIGI) {lhood *= BIG; lrnrm++;}
R79	if (lhood != 0.) while (lhood < BIGI) {lhood *= BIG; lrnrm++;}

   
File: incgammabeta.h  

L93	return (psig?(ans>0.0? 1.0-ans:-ans):(ans>=0.0? ans:1.0+ans));
R93	return (psig?(x>a1? 1.0-ans:-ans):(x>a1? ans:1.0+ans));

   
File: interp_1d.h  

L7	dj = MIN(1,(int)pow((Doub)n,0.25));
R7	dj = MAX(1,(int)pow((Doub)n,0.25));
L183	Int n=y2.size();

   
File: ksdist.h  

L17	if (z == 0.) return 0.;
R17	if (z < 0.042) return 0.;

   
File: markovgen.h  

L16	r = ran.doub()/cum[j][m-1];
R16	r = ran.doub()*cum[j][m-1];

   
File: qroot.h  

L21	b += (delb=(r*sc-s*rc)*dv);
	c += (delc=(-r*sb+s*rb)*dv);
R21	b += delb;
	c += delc;

   
File: quad3d.h  

R1	template <class T>
L3	Doub (*func3d)(const Doub, const Doub, const Doub);
R4	T *func3d;
L6	return func3d(xsav,ysav,z);
R7	return (*func3d)(xsav,ysav,z);
R10	template <class T, class Z1, class Z2>
L10	NRf3 f3;
	Doub (*z1)(Doub, Doub);
	Doub (*z2)(Doub, Doub);
	NRf2(Doub zz1(Doub, Doub), Doub zz2(Doub, Doub)) : z1(zz1), z2(zz2) {}
R12	NRf3<T> f3;
	Z1 &z1;
	Z2 &z2;
	NRf2(Z1 &zz1, Z2 &zz2) : z1(zz1), z2(zz2) {}
R22	template <class T, class Y1, class Y2, class Z1, class Z2>
L21	Doub (*y1)(Doub);
	Doub (*y2)(Doub);
	NRf2 f2;
	NRf1(Doub yy1(Doub), Doub yy2(Doub), Doub z1(Doub, Doub),
	Doub z2(Doub, Doub)) : y1(yy1),y2(yy2), f2(z1,z2) {}
R24	Y1 &y1;
	Y2 &y2;
	NRf2<T,Z1,Z2> f2;
	NRf1(Y1 &yy1, Y2 &yy2, Z1 &z1, Z2 &z2) : y1(yy1),y2(yy2), f2(z1,z2) {}
L32
	template <class T>
	Doub quad3d(T &func, const Doub x1, const Doub x2, Doub y1(Doub), Doub y2(Doub),
	Doub z1(Doub, Doub), Doub z2(Doub, Doub))
R34	template <class T, class Y1, class Y2, class Z1, class Z2>
	Doub quad3d(T &func, const Doub x1, const Doub x2, Y1 &y1, Y2 &y2, Z1 &z1,
	Z2 &z2)
L37	NRf1 f1(y1,y2,z1,z2);
	f1.f2.f3.func3d=func;
R38	NRf1<T,Y1,Y2,Z1,Z2> f1(y1,y2,z1,z2);
	f1.f2.f3.func3d=&func;