![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
clc clear close format long syms u v x y u1 u2 u3 u3 u x1 x2 x3 x4 x5 d L %int(int(int(int(1,x5,x4,1),x4,x3,1),x3,x2,1),x2,x1,1) %int(int(int(6,x3,x2+d,L),x2,x1+d,L-d),x1,0,L-2*d) a=3 b=200 y_P=linspace(0,0,a*b); y_N=linspace(0,0,a*b); index=0; t=1000000 for i=1:t u1=rand(); u2=rand(); v=sqrt(-log(u1)*2)*cos(2*pi*u2); if v>0 if v <=a index=round(v *b); if index == 0 index=1; end y_P(index)=y_P(index) +1; end else if v>=-a index=round(abs(v) *b); if index == 0 index=1; end y_N(index)=y_N(index) +1; end end end % x=[0,0.1,0.2,0.3,0.4] grid on hold x=linspace(0,a,a*b) plot(x,y_P*b/t) x=linspace(0,-a,a*b) plot(x,y_N*b/t)