c = 1
f0 = a * pow(b, 2) * c
f1 = a * b / c
return np.array([f0, f1])

def fit(x, a0, a1):
    res = a0 * x[0] + a1 * x[1]
return np.array([res])
x = getXval()
y = np.array([0.15, 0.34])
popt, pcov = curve_fit(fit, x, y)
print(popt, pcov)
scipy curve fit on a matrix