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”);
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”);