figure();
axis([-100, 100 -100 100]);
fontsize = 12;
xlabel('X', 'Fontsize', fontsize);
ylabel('Y', 'Fontsize', fontsize);
hold on
x_c = 0;
y_c = 0;
radii_plane = 80;
radii_vein = 4;
radii_inf_range = 12;
center_plane = [x_c, y_c];
viscircles(center_plane, radii_plane, 'color', 'b');
hold on
VNA = [0,80;69.2820323027551,40.0000000000000;-69.2820323027551,40.0000000000000;0,-80;-69.2820323027551,-40.0000000000000;69.2820323027551,-40.0000000000000];
SizeXY = size(VNA,1);
VNAx = VNA(:,1);
VNAy = VNA(:,2);
plot(VNAx, VNAy, '.', 'color', 'r', 'markersize', pi*6^2);
for i = 1:SizeXY
centers_nodeXY = [VNAx(i), VNAy(i)];
viscircles(centers_nodeXY, radii_inf_range, 'color', 'r');
end