vector a = (0,0,0);
vector b = (1,1,1);
vector n = (20,20,20);

mesh M = structured(n,a,b);

scene S = pov("plane.pov");
domain O = domain(S,inside(<1,0,0>));

mesh s = surface(O,M);

save(medit,"u", s);