如何在matlab中写“圆上任意两点的中点”的程序代码？？求指点呀！！！

知道怎么用蒙特卡洛仿真求圆的面积吗？下面有一个pudn上的代码，使用蒙特卡洛方法求圆心在原点，半径为1的圆的面积。你可以用类似的方法解决你的问题，只要把条件改成落点同时在三个圆里面就行了sita=0:0.01:2*pi;x=sin(sita);y=cos(sita);%计算半径为1的圆周上的点，以便作出圆周观察m=0;%在圆内在落点计数器x1=2*rand(1000,1)-1;%产生均匀分布于[-1,+1]直接的两个独立随机数x1,y1y1=2*rand(1000,1)-1;N=1000;%设置试验次数forn=1:N%循环进行重复试验并统计p1=x1(1:n);q1=y1(1:n);if(x1(n)*x1(n)+y1(n)*y1(n))<1%计算落点到坐标原点的距离,判别落点是否在圆内m=m+1;%如果落入圆中,计数器加1endplot(p1,q1,'.',x,y,'-k',[-1-111-1],[-111-1-1],'-k');axisequal;%坐标纵横比例相同axis([-22-22]);%固定坐标范围text(-1,-1.2,['试验总次数n=',num2str(n)]);%显示试验结果text(-1,-1.4,['落入圆中数m=',num2str(m)]);text(-1,-1.6,['近似圆面积S_c=',num2str(m/n*4)]);set(gcf,'DoubleBuffer','on');%双缓冲避免作图闪烁drawnow;%显示结果end