def model (par,x):
n = len(x)
res = []
for i in range(1,n):
    A0 = par[3] + (par[4]-par[1])*par[6] + (par[5]-par[2])*par[6]^2
    if(x[i] == par[6]):
        res[i] = A0 + par[1]*x[i] + par[2]*x[i]^2
    else:
        res[i] = par[3] + par[4]*x[i] + par[5]*x[i]^2
return res