(1) calculate intensity at vertices (I1, I2, I3) from illumination model.
(2) interpolate vertex intensities along edges (Ia, Ib)