function [TE,TM,Z_kr,A] = Eigenmoden(n,v,theta,theta1,f,r,phi)
%Eigenmode berechnet die Werte der Eigenmoden von TE,TM und TEM Wellen bei
%der Frequenz f und dem Abstand r vom Septum

% TE: Entspricht den Eigenmoden nach der Formel A6 aus der Publikation 

%Konstanten
mu0=4*pi*10^-7;
omega=2*pi*f;
mv=(v*pi)/(2*0.2618);
% phi=0.2618;

%Berechnung der Besselfunktion Z(kr);

ny=n+0.5;
kr=sqrt(ny*(ny+1));
Z_kr = besselh(ny,1,kr);

%Funktion U(phi) zur Beschreibung der phi-Abhängigkeit des Wellenfeldes

U_phi=cos(mv*phi);

%Berechnung der Legendrefunktionen P_n^mv(cos(theta)) und
%Q_n^mv(cos(theta)):


[LP3,nLP3,LP3ab,ALP3ab,nLP3ab] = legfonction(n,mv,cos(theta));
[LP4,nLP4,LP4ab,nLP4ab]        = my_legendre(n,mv,cos(theta1));
[LP5,nLP5,LP5ab,ALP5ab,nLP5ab] = legfonction(n,mv,cos(-theta));
P_costheta1    = LP4;     %my_legendre(n,mv,cos(theta1));
Pn_costheta1   = nLP4;    %my_legendre(n,mv,cos(-theta1));
Pab_costheta1  = LP4ab;  %derivest(@(theta) my_legendre(n,mv,cos(theta)),theta1);
Pnab_costheta1 = nLP4ab;  %derivest(@(theta) my_legendre(n,mv,cos(-theta)),-theta1);
Pn_costheta    = LP5;
Pnab_costheta  = ALP5ab;
P_costheta     = LP3;
Pab_costheta   = ALP3ab;

% Berechnet die Eigenmoden der TE-Wellen mit den Indizes xy
TE=((1i*(mu0*omega)./(r*sin(theta))* sqrt(kr)*Z_kr).*(Pnab_costheta1*P_costheta-Pab_costheta1*Pn_costheta)'/(Pnab_costheta1)*U_phi)';

%Berechnet die Eigenmoden der TM-Wellen mit den Indizes xy
TM=1/kr*(derivest(@(kr) sqrt(kr)*Z_kr,kr)*(Pn_costheta1*Pab_costheta-P_costheta1*Pnab_costheta))/(Pn_costheta1)*U_phi;
%TM=1/kr*derivest(@(kr) sqrt(kr)*Z_kr,kr)*(nLP2*LP1ab-LP2*nLP1)/(nLP2)*U_phi;
end