Liczba modów w światłowodzie planarnym

Manipulate[Plot[{k Tan[a k],-k Cot[a k],(Sqrt[R2-(k^2)])/.{R2->((n2^2)-(n1^2))/(λ^2)}},{k,0,5}, PlotRange->{0,5}, AspectRatio->Automatic],{a,0,30},{n2,1,2},{n1,1,2},{λ,0.01,3}]

Rozwiązanie równań na liczby falowe k i κ

a=10;
k0=1;
n2=1.6;
n1=1.5;
solsym=NSolve[{k Tan[a k]==κ,(κ^2)+(k^2)==((n2^2)-(n1^2))(k0^2), κ≥0, k≥0},{k,κ},Reals]
solantysym=NSolve[{-k Cot[a k]==κ,(κ^2)+(k^2)==((n2^2)-(n1^2))(k0^2), κ≥0, k≥0},{k,κ},Reals]

Mody symetryczne w falowodzie

Manipulate[Plot[(Exp[κ a]Cos[k a] Exp[κ z] HeavisideTheta[-(z+a)]+HeavisideTheta[-(z-a)]HeavisideTheta[(z+a)] Cos[k z]+
Exp[-κ z]HeavisideTheta[(z-a)] Exp[κ a]Cos[k a])/.solsym[[i]], {z, -20,20}, PlotRange->All],{i,1, Length[solsym],1}]

Mody antysymetryczne w falowodzie

Manipulate[Plot[(-Exp[κ a]Sin[k a] Exp[κ z] HeavisideTheta[-(z+a)]+HeavisideTheta[-(z-a)]HeavisideTheta[(z+a)] Sin[k z]+
Exp[-κ z]HeavisideTheta[(z-a)] Exp[κ a]Sin[k a])/.solantysym[[i]], {z, -20,20}, PlotRange->All],{i,1, Length[solsym],1}]