f=([z1,z2]f-[z0,z1]f)/(z2-z0),