%% Polje tockastega naboja
eps0=8.854e-12;
k=1/(2*pi*eps0);
Q=2e-4;
x=(-1:0.01:1)*1;

E=Q./(4*pi*eps0*x.^2).*sign(x)



plot(x,E)

%% Polje dveh pozitivnih nabojev
x=(-0.5:0.01:1)*1;
E1=Q./(4*pi*eps0*x.^2).*sign(x)
E2=Q./(4*pi*eps0*(x-0.5).^2).*sign(x-0.5)
plot(x,E1,x,E2,'LineWidth',2)

%% Polje pozitivnega in negativnega naboja
x=(-0.5:0.01:1)*1;
E1=Q./(4*pi*eps0*x.^2).*sign(x)
E2=-Q./(4*pi*eps0*(x-0.5).^2).*sign(x-0.5)
plot(x,E1,x,E2,'LineWidth',2)