% Zonenplatte

lambda_he_ne = 632.8*10^(-6);

% Verwendete Blenden
b_durchmesser(1)=60.68-59.24;
b_durchmesser(2)=60.38-59.24;
b_durchmesser(3)=60.09-59.64;
d_b_durchmesser = 0.01;

%Stellung der Blende
x_blende=679;
%brennebene_mikroskop=45;
brennebene_mikroskop=0;

% Blende1 vermessen
blende1(1)=767;
blende1(2)=750;
blende1(3)=743;
blende1(4)=736;
blende1(5)=732.5;
blende1(6)=729.5;
blende1(7)=726.5;

% Blende2 vermessen
blende2(1)=915;
blende2(2)=851;
blende2(3)=814;
blende2(4)=795;
blende2(5)=780;
blende2(6)=770;
blende2(7)=761;
blende2(8)=757;
blende2(9)=752;
blende2(10)=748;
blende2(11)=746;

% Blende3 vermessen
blende3(1)=1165;
blende3(2)=1037;
blende3(3)=943;
blende3(4)=897;
blende3(5)=865;
blende3(6)=843;
blende3(7)=827;
blende3(8)=812;
blende3(9)=806;
blende3(10)=796;
blende3(11)=791;

% genauigkeit beim ablesen
dx(1)=10;
dx(2)=10;
dx(3)=7;
dx(4)=5;
dx(5)=3;
dx(6)=2;
dx(7)=1;
dx(8)=1;
dx(9)=1;
dx(10)=1;
dx(11)=1;

imax=11;
i=1:imax;
x=0:0.1:imax;
z(i) = 1./(blende3(i)-x_blende-brennebene_mikroskop);

%fit
k=i'\z';

figure(1);
plot(i,z(i),'r*',x,k*x);
axis([0,11,0,0.03]);

%k=0.001;
r=sqrt(lambda_he_ne/k);
disp('gemessen:');
disp(b_durchmesser(3));
disp('errechnet');
disp(r*2);



f_alt=100;
z0=1000000;
d=-35;
%for z0=-10000:100:-200
for d=40:-1:-110
	r(i)=2*sqrt((2*i)./(1./(z0)+1./(blende3(i)-x_blende+d)).*lambda_he_ne);
   dr(i) = (i .* z0^2 .* lambda_he_ne)./(2*(blende3(i)+d+z0).^2.*sqrt((blende3(i)+d).*i.*z0.*lambda_he_ne)/(blende3(i)+d+z0)).*dx(i);
   
   rd=sum(r)/imax;
	f=sqrt(sum(1./dr(i).*(r(i)-rd).^2))/imax;
   if (f/rd)<f_alt
      disp([num2str(rd),' +/-',num2str(f/rd),' bei z0=',num2str(z0),' und d=',num2str(d)]);
      f_alt=(f/rd);
      r_i=r;
   end
   
%end
end