на Delphii 7. Если алгоритм правильно считает, то выложу исходник. Кто захочет - загонит в компонент или в IC.
p.s. Мне лень