function [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta)  
  
%[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta)  
%该函数用有限差分法求解有两种介质的正方形区域的二维拉普拉斯方程的数值解  
%函数返回迭代因子、迭代次数 ...
