# differential equations – ODE resolution system with additional parameter

I would like to solve a $$2 times 2$$ system of form
$$frac {d} {d theta} T = TA, quad T (0) = Id$$
or $$theta$$ is real and $$A$$ is of the form
$$A = begin {pmatrix} 0 & frac {e ^ {- i theta}} { lambda} \ frac {1} {36} e ^ {- i theta} left (9 lambda + 2 ( lambda-1) ^ 2 (6 cos { theta} + cos {2 theta} + 6) right) & 0 end {pmatrix},$$
with $$lambda$$ a free parameter in the unitary circle.

In particular, I'm interested in getting digital solutions to $$theta = 2 pi$$ depending on the additional parameter $$lambda$$. I'm pretty new using Mathematicaand that's what I've tried so far:

``````T[θ_] = {{T11[θ], T12[θ]}, {T21[θ], T22[θ]}};
A[θ_] = {{0, E ^ (- I θ) / λ}, {1/36 E ^ (- I θ) (9 λ + 2 (-1 + λ) ^ 2 (6 + 6 Cos[θ] + Cos[2 θ])), 0}};
sys = {T & # 39;[θ] == T[θ].A[θ]};
``````

The previous code defines the system I want to solve and I am now trying to solve it numerically. I first tried

``````NSol = NDSolve[{sys, T11[0] == 1, T12[0] == 0, T21[0] == 0, T22[0] == 1}, {T11[θ], T12[θ], T21[θ], T22[θ]}, {θ}, {θ, 0, 2 Pi}];
``````

which gives me the exit

``````NDSolve :: dupv: "Double variable θ found in NDSolve[<[<[<[<<1>>]. "
``````

I have also tried

``````Nsol2 = ParametricNDSolve[{sys, T11[0] == 1, T12[0] == 0, T21[0] == 0, T22[0] == 1}, {T11, T12, T21, T22}, {θ, 0, 2 Pi}, {λ}];
``````

which gives me exit $$T_ {11}, dots, T_ {22}$$ as ParametricFunctions dependent on each other and on $$lambda$$.

I do not know if this is the right approach and, if so, how to extract a numeric expression based on $$lambda$$ of the last output – all that I saw in the documentation are examples plotted for specific values ​​of the parameter. Any help is very appreciated.