雷诺方程的数值解法

				function [Delta]=all_zone_m()
				ps_location=evalin('base','ps_location');
				beta=evalin('base','beta');
				n=evalin('base','n');
				M=evalin('base','M');
				N=evalin('base','N');
				ps=[ps_location,beta];
				num=size(ps);
				for i=1:num(1)
				    m(i)=m_in(ps(i,1),ps(i,2),ps(i,3));
				end
				min1=sum(m);
				min2=mz(1,1,M+1);
				min3=mz(N+1,1,M+1);
				Dm=min1+min2-min3;
				D_ratio=Dm/min1;
				Delta=[min1,min2-min3,D_ratio,Dm];
				end