“…Based on a Harnack type theorem and linearization, this technique avoids the use of supconvolutions, that in presence of distributed sources produces several complicacies. The method can be very well adapted to nonhomogeneous two-phase problems to prove that flat (see below) or Lipschitz free boundaries of (1) are C 1,γ , when the governing equation is the same in both phases (see [16], [17], [18]). Throughout this section, L 1 = L 2 and this common operator will be denoted by L. Also, f 1 = f 2 = f.…”