count = 24;
dim_x = 1.2;
dim_y = dim_x;
x = linspace(-1. * dim_x, dim_x,count+1);
y = linspace(dim_y,-1. * dim_x,count+1);

[X Y] = meshgrid(x,y);

% X_vec = 0.2 * sin(2 * pi * X);
% Y_vec = 0.2 * sin(2 * pi * Y);

% Y_vec = -0.2 * X;
% X_vec = -0.2 * Y;

Radius = sqrt(X.^2+Y.^2);
% Angle

X_vec = 0.4 * 1. ./ Radius .* (Y./Radius);
Y_vec = -0.4 * 1. ./ Radius .* (X./Radius);


center = (count/2) + 1;

X_vec(center - 2:center +2, center - 1: center + 1) = zeros(5,3);
Y_vec(center - 2:center +2, center - 1: center + 1) = zeros(5,3);
X_vec(center - 1:center +1, center - 2: center + 2) = zeros(3,5);
Y_vec(center - 1:center +1, center - 2: center + 2) = zeros(3,5);

hold on;

plot (X , Y, "x", "color", "k", "markersize", 2. );

quiver(X, Y, X_vec, Y_vec, "linewidth", 1.5, "filled");




%plot (0, 0, "+", "markersize", 20, "linewidth", 5);


plot (0,0, "o", "markersize", 35, "linewidth", 5, "markerfacecolor", "r", "color", "red");

axis ([-1.05 * dim_x, +1.05 * dim_x, -1.05 * dim_y, +1.05 * dim_y], "square");
set (gca(), "linewidth", 2.);


xlabel ("x");
ylabel ("y");
title("Magnetische Feldstärke eines geraden stromdurchflossenen Leiters")

hold off;


print("magnetfeld.svg", "-dsvg", "-S1000,1000")

%printf("Weiter mit beliebiger Taste\n");
%pause;