function f=pij_cyl(r,sig,ngauss) % compute the reduced collision probability matrix in 1D cylindrical geometry % function f=pij_cyl(r,sig,ngauss) % (c) 2005 Alain Hebert, Polytechnique Montreal npij=size(r,2) ; vol=zeros(1,npij) ; rik1=0 ; for i=1:npij rik2=pi*r(i)^2 ; vol(i)=rik2-rik1 ; rik1=rik2 ; end f=zeros(npij) ; track=sybt1d(r,false,ngauss) ; % for ix=1:npij for i=1:ngauss w=track{ix,i}{1} ; seg=track{ix,i}{3} ; % complete Matlab statements after ... f(ix,ix)=f(ix,ix)+ ... ; tau=2*seg(1)*sig(ix) ; tau1j=0 ; for ip=ix+1:npij f(ip,ip)=f(ip,ip)+ ... f(ix,ip)=f(ix,ip)+ ... tauij=0 ; for jp=ip+1:npij f(ip,jp)=f(ip,jp)+ ... tauij=tauij+seg(jp-ix+1)*sig(jp) ; end tau=tau+2*seg(ip-ix+1)*sig(ip) ; tau1j=tau1j+seg(ip-ix+1)*sig(ip) ; end end end % for i=1:npij for j=i:npij val=f(i,j) ; f(i,j)=val/vol(i) ; f(j,i)=val/vol(j) ; end end