planetary modeling with power functions