function xn = newton_pop(x0,t,pop,tol,N)
% 2D Newton's Method
xn = x0; err = 1; k = 1;
while ((err > tol)&&(k < N))
    x1 = xn;
    xn = x1 - Jpop(x1,t,pop)\gpop(x1,t,pop)
    err = abs(norm(x1)-norm(xn))
    k = k+1
end
end

