% K_0 Asymptotic estimate for k_0, given the accuracy % level eps0 % function y=k_0(eps0,a,b) x=-((1-b)/a)*eps0^(1/a); if x<-exp(-1) error('invalid input data for k0') return end w=wofy1(x); y=ceil(-1-a*w/((1-b)*log(10)));