include "ffmatlib.idp" real Rdomaine = 1.0; // rayon du domaine real rayon = 0.02; // rayon du fil // fonction source (mu0*J=1) func real source(real &x, real &y){ if (sqrt(x*x+y*y)