function err = err_fit_degenerate(s,r0,u,D,G)

%This form of err_fit forces nugget to zero

%INPUTS: s  = sill 
%        r0 = max(D)

p = [r0,s,0]; %form parameter vector

v = spher(p,D);

err = sqsum(u,G,v) ;