Fitting experimental data with differential equations : Optimization

[Gribouillis]

You could define a functional that measures the distance between the two curves such as

X, Y = ...

def distance(a, B, G):
    x, y = Solver(a, B, G)
    return np.linalg.norm(y - Y)

Then minimize this function of 3 variables with scipy.optimize.minimize(distance, **args).