“Hubbard” equation function