pkg load msh
pkg load fem-fenics

fem_init_env ();
problem = 'Poisson';
fem_create_all (problem);

 

x = y = linspace (0, 1, 32);
msho = msh2m_structured_mesh (x, y, 1, 1:4);
mshd = fem_init_mesh (msho);

V  = fem_fs_Poisson (mshd);

bc = fem_bc (V, @(x,y) 0, [2, 4]);

f = fem_coeff ('f', @(x,y) 10*exp(-((x - 0.5)^2 + (y - 0.5)^2) / 0.02));
g = fem_coeff ('g', @(x,y) sin (5.0 * x));

A = fem_rhs_Poisson (V, bc);
b = fem_lhs_Poisson (V, f, g, bc);

u = A \ b;


func = fem_func ('u', V, u);

fem_plot (func);

fem_save (func, problem);