Q=1e-9;
r0=0.01;
d=0.02;
epsr=3;
eps0=8.854e-12;

% tvorimo vektor vrednosti radijev od r0 do r0+d:

r0d=linspace(r0,r0+d,100)