fdivround
((y1) + fdivround(((x) - (x1)) * ((y2) - (y1)), (x2) - (x1), n))
((y1) + fdivround(((int)(x) - (x1)) * ((y2) - (y1)), (x2) - (x1), n))