I have
l = IndexedBase(\'l\') i = symbols("i", cls=Idx) omega = symbols("omega", cls=Idx) var_l = l[i, omega] equation = Eq(-(b[omega] -