numerical recipes in C source codes

				void rsolv(float **a, int n, float d[], float b[])
				{
					int i,j;
					float sum;
				
					b[n] /= d[n];
					for (i=n-1;i>=1;i--) {
						for (sum=0.0,j=i+1;j						b[i]=(b[i]-sum)/d[i];
					}
				}
				/* (C) Copr. 1986-92 Numerical Recipes Software "U,6VV'. */