请教一个问题关于蒙特卡罗随机数模拟 求参数方程确定的爱心面积

问题遇到的现象和发生背景 蒙特卡罗随机数模拟 求参数方程确定的爱心面积

问题相关代码，请勿粘贴截图 %核心算法：取矩形框套住爱心，随机点落在爱心区域内的概率乘以矩形面积

sj_point=rand(100,2); %生成1000个随机二维随机数组 1002
%确定矩形大小
x1=32sj_point(:,1)+14; %确定随机点的x坐标范围[13,45)，即矩形宽为32.
y1=29sj_point(:,2)+13; %确定随机点的y坐标范围[12,42)，即矩形长为30.
love=find(); %找到落在心形图案里的点
point_number=length(love) %计算落在爱心里的个数
S=928point_number/100; %矩形面积乘以落在心里的概率即是爱心面积

运行结果及报错内容 love=find(); %找到落在心形图案里的点

我的解答思路和尝试过的方法 我没有办法确定点落在爱心里面的边界条件

我想要达到的结果