Вычисление параметров сложного теплообмена с использованием FreeFEM++

Уравнение Пуассона
или
dx(u)*dx(v)+dy(u)*dy(v) - (f*v)
Для каждого элемента сетки:
mesh Sh= square(4,4);
plot(Sh);

mesh Sh= square(10,10); // сетка
fespace Vh(Sh,P1); // пространство имен Vh в 2d
Vh u,v; // переменные
func f=cos(x)*y; // функция источника
problem Poisson(u,v)=
int2d(Sh)(dx(u)*dx(v)+dy(u)*dy(v))
-int2d(Sh)(f*v)
+on(1,2,3,4,u=0); // ГУ I типа
Poisson; // определяем problem для каждого FE
plot(u);