One way is this:
power = 1/2; shift = 0; add = 0; genericEquation = (y == m (x + shift)^power + b + add) instances = equation /. {{x -> 1, y -> 5}, {x -> 5, y -> 10}} solutions = NSolve[instances] finalEquation = genericEquation /. First[solutions]
