eps0=8.854e-12;

Q=2e-9;
R=5e-2; A=pi*R^2;
sigma=Q/A;

z=0:0.005:0.1;

cos_d=z./sqrt(R^2+z.^2);
Ediska=sigma./(2*eps0)*(1-cos_d);

q=Q/(2*pi*R);
Eobroca=q*R*z./(2*eps0.*(R^2+z.^2).^(3/2))

plot(z,Ediska,  z, Eobroca, '--','Linewidth',2)
legend('DISK','OBROC')

xlabel('Razdalja vzdolz Z osi   /  m ')
ylabel('Elektricna poljska jakost   /  V/m')

Eravnine=sigma/(2*eps0)*ones(length(z));
figure
plot(z,Ediska,  z, Eravnine,'--', 'Linewidth',2)
legend('DISK','RAVNINA')

xlabel('Razdalja vzdolz Z osi   /  m ')
ylabel('Elektricna poljska jakost   /  V/m')