setrval.m.txt

这是马尔可夫-蒙特卡罗算法的MATLAB源程序.

function DR = SetRVal(DR,inds1,inds2,r)
% DR = SetRVal(DR,inds1,inds2,r)
% set the value of the resistor from point inds1 (array-style indices) 
% to point inds2 to value r

% Colin Fox, 21 Jan 2003

nnps = round(sqrt(size(DR,1))); % number of nodes per side
DR(sub2ind([nnps nnps],inds1(2),inds1(1)),sub2ind([nnps nnps],inds2(2),inds2(1)))=r;