To make it plain and simple, if I have a holonomic constraint, that I want to treat using a lagrange multiplier, in any textbook I concern, they are just expressed as "$\lambda$" (omitting possible arguments). What I'd like to know is whether a lagrange multiplier is something like "$\lambda(t)$", that would be $\lambda$ to be an additional degree of freedom, whose time dependence is not known yet. Or instead, is the lagrange multiplier a (yet to exactly determine) extension to the lagrange function, and as such "$\lambda(q, \dot{q}, t)$"?
I will expand my thoughts here on why I think that both versions do work out: Starting with "$\lambda(q, \dot{q}, t)$", The complete Lagrangian of the System would be $$ \mathcal L = L(q, \dot{q},t) - \lambda(q, \dot{q},t) f(q,t) $$ Variations in the paths $\delta q$ and the requirement that $f(q(t),t)=0$ holds, will yield the correct EOM: $$ \frac{d}{dt} \frac{\partial \mathcal L}{\partial \dot{q}} (q(t),\dot{q}(t),t) - \frac{\partial \mathcal L}{\partial q} (q(t),\dot{q}(t),t) = \\ \frac{d}{dt} \frac{\partial \lambda }{\partial \dot{q}} (q(t),\dot{q}(t),t) - \frac{\partial \lambda}{\partial q} (q(t),\dot{q}(t),t))f(q(t), t) + \frac{\partial \lambda}{\partial \dot{q}}(q(t),\dot{q}(t),t) \frac{d}{dt}f(q(t),t) + \lambda(q(t), \dot{q}(t), t) \frac{\partial f}{\partial q} $$
Applying the constraints, $f = 0$ and $\frac{d}{dt} f = 0$, then we get:
$$ \frac{d}{dt} \frac{\partial L}{\partial \dot{q}} (q(t),\dot{q}(t),t) - \frac{\partial L}{\partial q} (q(t),\dot{q}(t),t) = + \lambda(q(t), \dot{q}(t), t) \frac{\partial f}{\partial q} $$ And $$ f(q(t), t) = 0 $$ which are just the equations that describe the motion of the system.
Alternatively, (as said), I can treat the multipliers as additional degrees of freedom in the configuration-space. The total lagrangian, dependent on $q$, $\dot{q}$, $\lambda$, $\dot{\lambda}$ (this dependence is just listed here for completeness, the total lagrangian won't depend on $\dot{\lambda}$), and t, is: $$ L(q, \dot{q},t) - \lambda(t) f(q,t) $$ Lagrange's equations then will be: $$ \frac{d}{dt} \frac{\partial L}{\partial \dot{q}}(q(t), \dot{q}(t),t) - \frac{\partial L}{\partial q} (q(t), \dot{q}(t), t) = \lambda(t) \frac{\partial f}{\partial q}(q(t), t) $$ And $$ f(q(t), t) = 0 $$
Both methods, although using different assumptions, yield the same equations of motion. Which one is more feasible? Are there cases where my reasoning doesn't work, which favor one of the two options I have given?