蒙特卡洛法求椭圆面积的MATLAB实现

syms x y
fimplicit((x.^2/4)+(y.^2/9)-1)
hold on
rectangle('Position',[-2,-3,4,6])
hold on
n=30000;
x=4*rand(1,n)-2;
y=6*rand(1,n)-3;
f=(x'.^2/4)+(y'.^2/9);
k=find(f<=1);
g=find(f>1);
num=numel(k);
X=x(k);
Y=y(k);
scatter(X,Y,1,'red')
num_2=numel(g);
X_1=x(g);
Y_1=y(g);
scatter(X_1,Y_1,1,'blue')
S=(num/n)*(4*6)