Для решения дифференциальных уравнений Пуассона (в частных производных второго порядка) и уравнений Лапласа в систему введены следующие функции:
Фbvalfit(vl, v2, х1, х2, xi, — устанавливает начальные условия для F, LI, L2, score) краевой задачи, заданной в векторах F, vl и v2
на интервале от х1 до х2, где решение известно в некоторой промежуточной точке xi;
Фге1ах(М1, М2, МЗ, — возвращает квадратную матрицу решения М4, М5, A, U, г) уравнения Пуассона для спектрального радиуса г, Фmultigrid(M, п) — возвращает матрицу решения уравнения Пуассона, у которого решение равно нулю на границах;
sbval(y, х1, х2, F, L.,score) — дает установку начальных условий для краевой задачи, определенной в символьном векторе F; вектор у — начальные условия на интервале от х1 до х2.