Sa fonction optimisation (après développement de Taylor) est du type:
E ~= sum [Vx Ix + Vy Iy + D]², où D est ta disparité (shift).
Tu optimises cette fonction par rapport à Vx et Vy pour trouver...