% POTENCIAL TOCKASTEGA NABOJA
Q=1e-8;
e0=8.854e-12;
r=-2e-2:1e-3:2e-2;
V=Q./(4*pi*e0.*abs(r));    %potencial
E=sign(r)*Q./(4*pi*e0.*r.^2); % funkcija sign() poskrbi za pravilen predznak polja
plot(r,V); xlabel('Razdalja  / m'); ylabel('Potencial   / V');
figure;
zero=zeros(length(r),1); % vektor ni?el potrebujemo za izris linije ni?le polja
plot(r,E,r,zero); xlabel('Razdalja  / m'); ylabel('Polje   / V/m');