How to estimate pi within a 1x1 square

I was asked to create a function that plots points randomly in a [0,1]x[0,1] square. Then to estimate Pi based on the number of points that fall within the 1st quadrant unit circle.

• So: points within the unit circle / points total.
• This should roughly equal pi.I'm stuck on how to create that unit circle boundary and how to look at the # number of points within boundary.This is what I have so far:

function estpi=estimate(npts)
x=zeros(1,npts);
y=zeros(1,npts);
for i=1:npts
    x(i)=rand(1);
    y(i)=rand(1);

end

plot(x,y,'o')

end
%number of points inside circle/ total points

Answer 1

[Walter Roberson]

Remember, inside or on a circle, x^2 + y^2 <= r^2

[Matthew Lozancich]

Thanks for the hint!! figured it out :)