clear all; clc;

n = 1:3;

figure()
hold on
z = 0;
for x = -1:.5:1
    for y = -1:.5:1
        E = (9e9*3e-6) / (sqrt(x^2+y^2))
        quiver3(x, y, z, E(1), E(2), E(3))
    end
end