Dear Yves,
First of all, thank you very much for your reply!
You're right: the placing of polynomials was not normal, I am really grateful for that notice. However, the overall picture is basically the same: I understand how to code one distinct term, but for now I do not know how to combine that terms under different regions of integration in a single equation. For example:
http://image.ibb.co/mBSUfb/La_Te_X_Example_2.png.
Here we have basically an equation, in which a sum of an integral over region and an integral over region's boundary is equal to other integral over the same region.
(Actually, this situation emerges from the need to reformulate weak form for the general linear elliptic PDE, in which there is a Hessian in an integral over region; I know that getfem++ can operate with Hessians, so, maybe, I should try to leave it in "as is" condition)
As I understand, I have to use the generic assembly bricks for that, as I do not see anything comparable in the section about high-level generic assembly procedures.
May I ask you one more question: is there any example of generic matrix computation, when I can directly calculate the entry of the model matrix and the right-hand side iteratively? I use the Python interface, but any example will do. Thank you again for your time!
(p.s. Sorry, I'm sending that again since I forgot to check "reply to all" option at the first time)