function PSet4Q3_vect_q

N = 1000;
tol = 1e-5;

% *** ENTER TRUE SOLUTIONS ***
r1 = ???
r2 = ???
r3 = ???

% generates N equally spaced points between -5 and 5
xv = linspace(-5,5,N);
yv = linspace(-5,5,N);

% matrix X stores the x-values of the initial guesses, Y the y-values
[X,Y] = meshgrid(xv,yv);
X0 = X;
Y0 = Y;

% *** WRITE CODE FOR NEWTON ITERATION HERE ***


% find the distance between the numerical solution and each exact solution
Nrm1 = sqrt((X-r1(1)).^2+(Y-r1(2)).^2);
Nrm2 = sqrt((X-r2(1)).^2+(Y-r2(2)).^2);
Nrm3 = sqrt((X-r3(1)).^2+(Y-r3(2)).^2);

% color the initial guess according to the solution it finds
C = zeros(N,N);
C(Nrm1 < tol) = 1;
C(Nrm2 < tol) = 2;
C(Nrm3 < tol) = 3;

figure(1)
surf(X0,Y0,rand(N,N),C,'edgecolor','none'); view([0,90])