求解答，有没有大佬帮忙看看这段代码（用来生成某种符合特定分布的样本点的），我看不懂
function output1 = UND(input1,input2) %按照randn函数的使用方法使用本函数
e=0;
sigma=0.81;
a=-4;
b=4;
jd=0.01;
x1=a:jd:b;
x2=a:jd:b-jd;
y1=1./(1+exp(pi*(e-x1)/sigma/sqrt(3)));
for i=1:1:length(x1)-1;
y2(i)=y1(i+1)-y1(i);
end
yy=min(y2);
y3=round(y2/yy);
a=sum(y3);
y4=[];
for i=1:1:length(y3)
z1=ones(1,y3(i))*x2(i);
y4=[y4,z1];
end
for i=1:1:input1
for j=1:1:input2
z=round(rand()*a);
output1(i,j)=y4(1,z);
end
end
end