import?matplotlib.pyplot?as?plt
from?scipy.optimize?import?curve_fit
import?numpy?as?np
#?np.seterr（invalid=‘ignore‘）
x2?=?np.array（[0.0000?0.0100?0.0200?0.0400?0.0600?0.0900?0.1200?0.1600?0.2000?0.2400]）
y2?=?np.array（[0.721?0.676?0.659?0.619?0.587?0.545?0.496?0.471?0.459?0.431]）


def?func（x?a?b?c?d）:
????return?a*np.log（1+b*x）+c*x+d


pop