Исправление ziemin, 01.12.13 00:46 (текущая версия) :
А тута:
(x2-x1)*(tx-Dx) / (y2-y1) + ty;
Будет
b*(tx-Dx)*sqrt(x1)/a*sqrt(tx) +ty
корни тоже вычислить заранее
Исходная версия ziemin, 01.12.13 00:44:
b*(tx-Dx)/a +ty