Rao:Neutrophil:Dmatrix: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
No edit summary
 
No edit summary
 
Line 1: Line 1:
function A = Dmatrix(N)
<nowiki>function A = Dmatrix(N)
%Diffusion Matrix
%Diffusion Matrix
D = 1;
D = 1;
Line 11: Line 11:
end
end
A = (N-1)*(N-1)*A;
A = (N-1)*(N-1)*A;
A = sparse(A);
A = sparse(A);</nowiki>

Latest revision as of 14:02, 18 August 2006

function A = Dmatrix(N) %Diffusion Matrix D = 1; A = zeros(N,N); A(1,1) = -2*D; A(1,2) = D; A(1,N) = D; A(N,N) = -2*D; A(N,N-1) = D; A(N,1) = D; for i=2:N-1 A(i,i) = -2*D; A(i,i+1) = D; A(i,i-1) = D; end A = (N-1)*(N-1)*A; A = sparse(A);