function OUT = var_spher_old(y,b,M)

%VAR_SPHER estimates a spherical variogram for a given data set
%using nonlinear least squares.
%
%INPUTS: (i)   y = [r,s,n] = initial estimates of [range,sill,nug].
%        (ii)  b = number of bins to use in empirical variogram
%        (iii) M = (nx3)-matrix of locations and data values,
%                  M(i) = (Xi,Yi,Zi), i = 1:n

%
%OUTPUTS: OUT = cell(4,1) with
%         (i)   OUT{1}= [R S N] (= estimated [range,sill,nugget])
%         (ii)  OUT{2}= distance values, (Di, i=1:b) 
%         (iii) OUT{3}= empirical variogram values Gi at Di
%         (iv)  OUT{4}= spherical values Vi at Di

DAT = variogram(M,b) ;

D = DAT(:,1) ;

G = DAT(:,2) ;

p = var_fit(y,D,G) ;

OUT = cell(4,1) ;

OUT{1} = p ;

OUT{2} = D ;

OUT{3} = G ;

OUT{4} = spher(p,D) ;