power_equation_base function

use power equation to fit observed values