% PRIMER IZRISA POLJA IN POTENCIALA V SFERI?NEM KONDEZATORJU S PROGRAMOM MATLAB
e0=8.854e-12;
U=2000; rn=2e-3; rz=5e-3;
Q=U*4*pi*e0/(1/rz-1/rn);
R=0:1e-5:rz;
E=abs(Q./(4*pi*e0.*R));
V=Q/(4*pi*e0).*(1./R-1/rn);
for i=1:1:length(R)
    if R(i)<rn 
    V(i)=0;
    E(i)=0;
    end
end
[ax ax1 ax2]=plotyy(R,V,R,E,'plot'); 
axes(ax(1)); ylabel(' Potencial [V]'); xlabel('Razdalja    [m]')
axes(ax(2)); ylabel(' El. poljska jakost [V/m]');
set(ax1,'LineStyle',':')
set(ax1,'Linewidth',3)
set(ax2,'Linewidth',2)