Программы на с++ для диплома

Мне нужно на c++ реализовать метод конечных объемов в цилиндрических координатах (через матрицу в профильном формате и решение слау). Есть написанный метод конечных элементов, могу скинуть для примера, нужно по аналогии реализовать МКО (сборку матрицы и правой части + учет краевых). Помимо кода нужно описание теории, как эти формулы получились (3-5 стр)