function [k1,k2] = dist_double_power_calib(s1,f1,s2,f2)

low = 1;
high = 6;

k0 = 1;

%k1 = fminbnd('dist_double_power_criterion_new',low,high,[],s1,f1,s2,f2) ;

k1 = fminsearch('dist_double_power_criterion',k0,[],s1,f1,s2,f2) ;

if f2^k1 > 10^(-16)
    
    k2 = log(s2)/log( 1 - f2^k1);
    
else
    disp('***************************************************');
    disp(' ');
    disp('VALUE OUT OF RANGE: k2 is too large to calculate')
    disp(' ');
    disp('***************************************************');
    error(' ');
end