c数值算法的源码

				#include "nrutil.h"
				#define MACC 3
				
				void mpdiv(q,r,u,v,n,m)
				int m,n;
				unsigned char q[],r[],u[],v[];
				{
					void mpinv(),mpmov(),mpmul(),mpsub();
					int is;
					unsigned char *rr,*s;
				
					rr=cvector(1,(n+MACC)					s=cvector(1,n+MACC);
					mpinv(s,v,n-m+MACC,m);
					mpmul(rr,s,u,n-m+MACC,n);
					mpmov(q,&rr[1],n-m+1);
					mpmul(rr,q,v,n-m+1,m);
					mpsub(&is,&rr[1],u,&rr[1],n);
					if (is) nrerror("MACC too small in mpdiv");
					mpmov(r,&rr[n-m+1],m);
					free_cvector(s,1,n+MACC);
					free_cvector(rr,1,(n+MACC)				}
				#undef MACC