Тут пообщался с VTO. Он подсказал про очень хороший пакет для решения систем линейных уравнений. Я пользовался самопальной куда как менее эффективной версией, но она сейчас неактуальна. На Западе используется очень широко. Его перевод на PL может помочь популяризации языка. Я сам проверил — комплексная система 5000-го порядка колется за 3 сек. (Гауссу нужно 500) и это безо всяких ухищрения, которые тоже возможны. Т.к. мое нынешнее знание PL на низком уровне, я могу поучаствовать в переводе комментариев (собственно, перевести их все и сделать две версии: английскую и русскую). Оптимально написать еще и драйвер, т.к. количество параметров велико и не всем они нужны, да не все их и поймут. Объём — 1700 строк кода с комментариями (которых сотни три). Кто возьмется?
Там четыре версии: действительная и комплексная одинарной и двойной точности, но они практически идентичны, я даже понимаю, как их свести в фортране в одну (в PL, думаю, так не удастся).