include “MUMPS”

// Parameters func f = 1.;

// Mesh int nn = 25; //Mesh quality mesh Th = square(nn, nn);

// Fespace func Pk = P2; fespace Uh(Th, Pk); Uh u;

// Plot plot(u, nbiso=30, fill=true, value=true, cmm=“A”);