num=20000??;%2000個數據
x0=abs（sin（randn））;y0=0;%取x0為[01]之間的隨機數
lamda=3.7
%for?lamda=3.1:0.1:4;
?%???a=a+1;
xn=x0;
for?n=1:num
???xn=lamda*xn*（1-xn）;%先進行循環迭代以達到穩定的定常解
???h（n）=xn;
end
%h=zeros（1num）
k=zeros（1num）;kk=zeros（1num）;
%h為序列的數值，k為對應的度，kk為各個度具有的點數
for?i=1:num-1??????????%這個大循環獲得了各點的度值
????k（i）=k（i）+1;
????k（i+1）=k（i+1）+1;%相鄰的節點度值加一
????for?j=2:num-i
????????f1=（h（i+j）-h（i））/j;
????????cont1=0;
????????for?u=1:j-1%判斷ij之間的點是否阻礙了二者的可視
????????????f2=（h（i+u）-h（i））/u;
?????????