Energy integral method for gravity field determination from satellite orbit coordinates