Когда-то давно я писал топик про REDUCE — open source систему для всяческих вычислений.

Тут нужно было накидать вот такой расчет эффективной позиции линзы ELP после имплантации ИОЛ:

doctor-hill.com/iol-main/elp.html

Если ссылка не работает, то я сохранил тут.

Вещь конечно специфическая, но может кому пригодится:



Непосредственно программулька:

Обозначения я упростил, POST = APostRx, IOL = IOLe, L = ALo.
Данные K, L, IOL, V — взяты из первого примера статьи.

X := 1336/(1000/(1000/POST-V)+K);
A := IOL;
B := -1*IOL*(L+X);
C := 1336*(L-X)+IOL*X*L;
IOL := 22.5;
L := 22.48;
K := 46.685*0.98765431;
V := 12;
POST := -0.5;
ELP := (-1*B-SQRT(B*B-4*A*C))/(2*A);
on rounded; X;
on rounded; A;
on rounded; B;
on rounded; C;
on rounded; ELP;
end;

Просто сохраняем в файл и открываем в REDUCE.