return (fxder[3] * xi[0] ** 4. + fxder[0] * xi[3] +