ydata = np.array([0.699369, 0.700462, 0.695354, 1.03905, 1.97389, 2.41143, 
    1.91091, 0.919576, -0.730975, -1.42001])

def model0(x, p1, p2):
    pass
How to calculate the likelihood of curve-fitting in scipy?