function J = diabetes_err(p,td,gd)
% Least squares error
y = GTT(td,p);
err = y - gd;
J = err*err';
end

