clear all
Q1=1e-6; 
Q2=-3e-6;	
e0=8.854e-12;
k=1/(4*pi*e0);  

X=(-5:0.5:20)*1e-2
%Q1 postavimo v koord. izhodišce, Q2 pa 12 cm stran v smeri poz. x osi
E1=k*Q1*X./X.^3;	% operator za deljenje vektorja je "./"
XT=12*1e-2*ones(1,length(X)); 	%vektor z vrednostmi pozicije lektrine Q2
Xraz=X-XT;
E2=k*Q2*Xraz./Xraz.^3
E=E1+E2;

plot(X,E1,X,E2,X,E,'o')
xlabel('X os   [cm]');
ylabel('E1, E2 in E1+E2     [V/m]');

% Če želimo izpisati točno lego:
a=find(E<0);
X(a(1))
x=int2str(1e3*X(a(1)))
answer=['Polje je nic pri X='x ' mm!'];
disp(answer)