Решение уравнения Пуассона с нулевыми граничными условиями
Листинг 11.7.Решение уравнения Пуассона с нулевыми граничными условиями
В первой строке листинга задается значение м=32, в двух следующих строках создается матрица правой части уравнения Пуассона, состоящая из всех нулевых элементов, за исключением одного, задающего расположение источника. В последней строке матрице с присваивается результат действия функции multigrid. Обратите внимание, первый ее аргумент сопровождается знаком "минус", что соответствует записи правой части уравнения Пуассона (11.11). Графики решения показаны на Рисунок 11.17 и 11.18 в виде трехмерной поверхности и линий уровня соответственно.
