% ÀîÈÊÃñclc;clear all;close allmax_j=193；max_i=81; d11=0.005 ; d13=0.005； d14=0。001；d12=0。001;d15=0。0001； d16=0。0001 ; d17=0.0001 ； d18=0.0001；n=0。7; % ׿϶¶Ècdd=0。001； ％ Pcdà ú ¤¹ÌÏɳºË¸ÅÂÊnumtotal_need=（1-n）＊max_i*max_j;numsoild=0；arrgrid=zeros（max_i,max_j)； while numsoild2+1 i=index_i；j=index_j—1; if arrgrid(i,j)==0 &＆ rand(）〈d13 numsoild=numsoild+1；arrgrid(i,j)=1; soild（numsoild,1）=i；soild（numsoild，2)=j; end end ％¹ àòú ¤ 4·òÌÏÄÚºËÏÓÒɳ½Ï if index_i>2+1 i=index_i-1；j=index_j； if arrgrid（i,j)==0 && rand（)〈d14 numsoild=numsoild+1；arrgrid（i,j）=1； soild(numsoild,1)=i；soild（numsoild，2）=j; end end %àòú ¤ 5·ò¹ÌÏÄÚºËÏÓÒɳ½Ï if index_j〈max_j—1 && index_i
