What is wrong with the following code snippet for fitting a quadratic model using curve_fit?

What is wrong with the following code snippet for fitting a quadratic model using curve_fit?

import numpy as np
from scipy.optimize import curve_fit

def quad(x, a, b, c):
    return a * x**2 + b * x + c

xdata = np.array([1, 2, 3, 4])
ydata = np.array([3, 7, 13, 21])

popt, pcov = curve_fit(quad, ydata, xdata)
print(popt)