La simulation ne sera pas réalisée par manque de temps. Nous en donnons l'idée derrière : Nous avons déjà écrit l'équation de l'indice optique dépendante de la position, du temps et de la température en prenant en compte la perturbation : Modélisation milieu perturbé (eau sucrée).

Cette équation donne un gradient en 2 dimensions. L'obtention d'une équation différentielle régissant la propagation du rayon n'est pas évidente et ne sera pas faite. À la place, nous avons décidé de résoudre la propagation à l'aide de la loi de Snell-Descartes et d'une subdivision de l'espace en petits carrés. Nous calculerons l'indice optique uniquement dans les carrés que le rayon rencontrera.