clear?clc
N=input（‘裂隙組數=‘）;
chang=input（‘生成域長=‘）;
kuan=input（‘生成域寬=‘）;
mianji=chang*kuan;?????????????????????????????????????????????????????????%面積
LXDD=[];
for?k=1:N
????midu=input（‘密度=‘）;
????zouxjz=input（‘走向均值=‘）;
????zouxbzc=input（‘走向標準差=‘）;
????jicjz=input（‘跡長均值=‘）;
????jicbzc=input（‘跡長標準差=‘）;
????xikanjz=input（‘隙寬均值=‘）;
????xikanbzc=input（‘隙寬標準差=‘）;
????n=mianji*midu;?????????????????????????????????????????????????????????%條數
????zhongxidian=unifrnd（-chang/2chang/2n2）;?????????????????????????????%中心點
????zouxiang=normrnd（zouxjzzouxbzcn1）;??????????????????????????????????%走向
????m=jicjz;
????v=jicbzc^2;
????mu=log（（m^2）/sqrt（v+m^2））;
????sigma=sqrt（log（v/（m^2）+1））;
????jichang=lognrnd（musigman1）;?????????????????????????????????????????%跡長
????xikuang=normrnd（xikan