eps0=8.854e-12;
k=1/(2*pi*eps0);
l=0.1; % dolzina palice
r=[0.1,1,2,5,10,100]*1e-2;
Q=2e-6;
q=Q/l;

Epremice=q./(2*pi*eps0*r)

cos_alfa=(l/2)./sqrt(r.^2+(l/2)^2)
Edaljice=q./(4*pi*eps0.*r)*2.*cos_alfa

Etocke=Q./(4*pi*eps0*r.^2)


loglog(r,Epremice,'--',r,Edaljice, r,Etocke,':', 'Linewidth',2)
xlabel('Razdalja od naelektrene daljice   /  cm ')
ylabel('Elektricna poljska jakost   /  V/m')