liti6.m

来自「此文件包括蒙特卡罗方法的基本思想及应用」· M 代码 · 共 7 行

Randnum=unifrnd(0,2*pi,1,10000);
xRandnum=cos(Randnum)
[Y,II] =sort(xRandnum)

yRandnum=sin(Randnum)
plot(xRandnum(II),yRandnum(II),'.')