У меня нет CoDeSys 3.5, посмотреть не смогу.
Было бы хорошо, если бы исходники были опубликованы в тексте сообщения, было кратко приведено описание применения - ведь нет смысла в каждом цикле выполнять пересчёт полиномов, значит он происходит или в отдельном ФБ или по отдельной команде.
Также интересно видеть формулы - как расчёта полиномов, так и вычислений функции, ограничения входных и выходных значений.

Без этого - применение чужого кода - "прыжок веры" Индианы Джонса.

В своё время хотел применить интерполяционные формулы Лагранжа и Ньютона для корнеизвлекающей зависимости расхода от перепада давления. Хорошо, что сначала построил график на компьютере - в узлах значения совпадали, а между ними формировались "вертикальные волны", т.е. между узлами результат интерполяции был со значительной ошибкой. Думал о кубических сплайнах, но не стал тратить время и сделал кусочно-линейной интерполяцией.
Понятно, что со сплайнами будет лучше, но всё же есть смысл проверить.
Была сетка значений:
dP, кПа 0.610 2.450 5.540 9.860 14.100 22.400 30.500 39.800 50.300
F, т/ч 5 10 15 20 25 30 35 40 45