Метод матричной прогонки для уравнения Пуассона

[Юрикк]

Помогите с методом матричной прогонки.
Я нашел готовый код матричной прогонки, но не могу его разобрать, чтоб переделать под себя. Код на с++ билдере в консоли.

Этот код написан по книге Самарский, Гулин - Чисельные методы, находится начиная с 411 страницы

Здесь решение выводится в ексель - сама матрица и решение уравнения.
Я не могу найти, где здесь менять и как вводить начальные данные, чтоб эта матрица и уравнение небыли статическими.

То-есть мы имеем такое уравнение(подчеркнуто красным):

http://s019.radikal.ru/i605/1505/79/df716657b908.png

и начальные условия, мне нужно сделать это на форме и чтоб можно было задавать параметры и функции некоторые начальные параметры, Ф, мью, для начала хотя-бы числами. Потом уже можно будет писать с символами и обрабатывать их.
В этом же коде не могу понять, где и как можно добавить такой функционал сюда, даже в консоли.
Если кто-то сможет с этим управиться, уточню задание в личку.