вот по пунктам как я понимаю:
1. вычисляем dХ и dY
2. определяем которая из них (дельт) больше - например dy
3. переменная err - на входе цикла равна 0
4. делаем цикл на столько итераций сколько в бОльшей дельте for n=0 to (y-1)
5. ставим точку и увеличиваем ту координату, дельта которой больше - y=y+1
6. err = err + dx - прибавляем меньшую дельту к переменной
7. if err > dy then err = err-dy:x=x+1: end if;