power_equation_all function

use power equation to fit observed values