y.append(dfun(fxder, xi))