Thank you but it seems it doesn't work. What i would like to do is this:
a[1][\[Omega]_] = \[Omega]
a[2][l_] = l
V[x_, a_[1], a_[2]] := 1/4*a[1]^2*x^2 + a[2]*(a[2] + 1)/x^2
V[x, a[1] + 1, a[2]]
This works, but after i found a k value (like k=1) i would be able to eliminate a[2] dependence from V, defining a new function.
k := 1
G[x_, a_[k]] := V[x, a[1], a[2]]
G[x, a[k] + 1]
This doesn't work