Добрый день.
Мне кажется в данном скрипте есть ошибка.
1. разместил квадратики так как было удобно мне для восприятия (связи между ними не менял естественно).
2. см рис. В самой формуле fs1 идет расчет по srs, а связь стоит от szs (на рис отмечена красным), то есть берется только сумма значений на интервале "к", а не их средняя величина.

Оставшуюся часть реализации формулы не смотрел, если не затруднит проверьте ее.

добавлено:
странную вещь обнаружил - значение fs1 остается неизменным как при наличии связи от srs так и при отсутствии её.