gusucode.com > 蒙特卡洛法求椭圆面积项目 > chengxu.m

    clear
n=500000;   % 随机生成n个点
x=rand(1,n);
y=rand(1,n);
x=2.*x;
y=3.*y;
r=(1/4).*x.*x+(1/9).*y.*y;  % 计算点到圆心的距离
m=find(r<=1);  % 找到落在椭圆里面的点
mm=length(m);  % 落在椭圆里面点的个数
S=(mm/n)*24    % 求椭圆面积