%Laden van gegevens

    XX=load('-ascii', 'GOP2.txt');
    omega_uexp=XX(:,1)*2*pi;		%1e kolom
    domega_uexp=XX(:,2)*2*pi;		%2e kolom
    Aexp=XX(:,3)/1000;			%3e kolom
    dAexp=XX(:,4)/1000;			%4e kolom
    phiexp=XX(:,5);         %5e kolom
    dphiexp=XX(:,6);        %6e kolom
    N=length(freqexp);

%fitparameters  
    
    iplot=2;
    n1=50;
    n2=50;
    n3=50;
    n4=50;
    n5=50;
    omega_0testvec=linspace(0,1e5,n1)';
    b_cvec=linspace(0,10000,n2)';
    F_0vec=linspace(0,10000,n3)';
    bvec=linspace(0,10000,n4)';
    mvec=linspace(1e-3,1,n5)';
    chi2mat=zeros(n2,n3);

    
%fitten van fase(omega_uexp) - fase(b,omega_0,m)

    for i1=1:n1,
        for i4=1:n4,
            for i5=1:n5, 
              omega_0test=omega_0testvec(i1);
              btest=bvec(i4);
              mtest=mvec(i5);
              faseTheoretischvec=fase_functie(btest,omega_uexp,mtest,omega_0test); %Theoretische functie
              chi2=sum((abs(faseTheoretischvec-phiexp)).^2);
              chi2mat(i1,i4,i5)=chi2;	
            end;%end for i5
        end;%end for i4
    end;%end for i1
  
%bepaling van beste fitparameters voor minimale chi2 

    [Y1,I1]=min(chi2mat);	%I1 = beste indices van omega_0, 1e coordinaat
    [Y2,I2]=min(Y1);        %I2 = beste indices van 2e coordinaat
    [Y3,I3]=min(Y2);        %I2 = beste indices van 3e coordinaat

    omega_0best=omega_0testvec(I1(I2(I3)))
    bbest=bvec(I2(I3))
    mbest=mvec(I3)


%fitten van Aexp(omega_uexp)-A(omega_0,b_c,F_0)  