< C语言数值算法程序大全>>配套程序

				void rsolv(a,n,d,b)
				float **a,b[],d[];
				int n;
				{
					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];
					}
				}