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 % 求椭圆面积