Есть решение ДУсЧП, представленное в виде ряда. Хочу оценить погрешность решения, полученного МКЭ. Это можно сделать с помощью кода
Код:
errornorm(u, uh, norm_type='l2')
где uh -- МКЭ-решение, а u -- точное решение. В качестве u может быть подставлено либо Function, либо Expression. Не понятно,
как создать Expression, отвечающий конечной сумме с произвольным количеством слагаемых.
Были такие идеи:
1)Создать для каждого слагаемого по Expression'у и суммировать их. Однако, в результате получается объект класса Sum, который не принимается функцией errornorm.
2)Создать строку, в которой будет записана искомая сумма в синтаксисе C++, и создать по ней Expression. Однако, такой способ оказывается неоправданно затратным по памяти.
Неужели в таком мощном инструменте нет возможности задавать функции, представимые суммой?
