function p = phi_tot(x, a, b, e, u0)
    p = e*tanh(-(x+a)/b)+u0;
end