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