%Franziska Minolts, Stand 06.10.11
%--------------------------------------------------------------------------
function us=test2(df,dsp,kla,kta,ka,nmax)                               

jn_ka(1,1)=sin(ka)/(ka);                                                    %Besselfunktion 1. Ordnung für den Wert n=0 für ka  
jn_kla(1,1)=sin(kla)/(kla);                                                 %Besselfunktion 1. Ordnung für den Wert n=0 für kla
jn_kta(1,1)=sin(kta)/(kta);                                                 %Besselfunktion 1. Ordnung für den Wert n=0 für kta
yn_ka(1,1)=-cos(ka)/(ka);                                                   %Besselfunktion 2. Ordnung für den Wert n=0 für ka
sqrt_ka=sqrt(pi/(2*ka));                                                    %Hilfsvariable zur Berechnung der Besselfunktion für ka
sqrt_kla=sqrt(pi/(2*kla));                                                  %Hilfsvariable zur Berechnung der Besselfunktion für kla
sqrt_kta=sqrt(pi/(2*kta));                                                  %Hilfsvariable zur Berechnung der Besselfunktion für kta
ajn_ka(1,1)=(cos(ka)*ka-sin(ka))/(ka^2);                                    %1. Ableitung der Bessel-Funktion 1. Ordnung nach ka für n=0
ajn_kla(1,1)=(cos(kla)*kla-sin(kla))/(kla^2);                               %1. Ableitung der Bessel-Funktion 1. Ordnung nach kla für n=0
ajn_kta(1,1)=(cos(kta)*kta-sin(kta))/(kta^2);                               %1. Ableitung der Bessel-Funktion 1. Ordnung nach kta für n=0
ayn_ka(1,1)=(sin(ka)*ka+cos(ka))/(ka^2);                                    %1. Ableitung der Bessel-Funktion 2. Ordnung nach ka für n=0

for n=1:1:nmax-1
    m=n+1;                                                                  %Hilfsparameter zum Speichern der Berechneten Wert in Matrix
    k=n+0.5;                                                                %Hilfsparameter zur Berechnung der Besselfunktion 
    jn_ka(m,1)=sqrt_ka*besselj(k,ka);                                       %Berechnung der Bessel-Funktion 1. Ordnung für ka
    jn_kla(m,1)=sqrt_kla*besselj(k,kla);                                    %Berechnung der Bessel-Funktion 1. Ordnung für kla
    jn_kta(m,1)=sqrt_kta*besselj(k,kta);                                    %Berechnung der Bessel-Funktion 1. Ordnung für kta
    yn_ka(m,1)=sqrt_ka*bessely(k,ka);                                       %Berechnung der Bessel-Funktion 2. Ordnung für ka
    
    ajn_ka(m,1)=abesj(n,ka);                                                %1. Ableitung der Bessel-Funktion 1. Ordnung nach ka
    ajn_kla(m,1)=abesj(n,kla);                                              %1. Ableitung der Bessel-Funktion 1. Ordnung nach kla
    ajn_kta(m,1)=abesj(n,kta);                                              %1. Ableitung der Bessel-Funktion 1. Ordnung nach kta 
    ayn_ka(m,1)=abesy(n,ka);                                                %1. Ableitung der Bessel-Funktion 2. Ordnung nach ka
end;
   
for n=1:1:nmax
    m=n-1;                                                                  %Hilfsparameter für die Berechung der einzelnen Koeffizienten
    dn=-jn_ka(n,1)/yn_ka(n,1);                                              
    an_ka=-ka*ajn_ka(n,1)/jn_ka(n,1);                                       
    an_kla=-kla*ajn_kla(n,1)/jn_kla(n,1);                                   
    an_kta=-kta*ajn_kta(n,1)/jn_kta(n,1);                                   
    bn=-ka*ayn_ka(n,1)/yn_ka(n,1);                                          
   
    zaehler=-kta^2*((an_kla/(an_kla+1))-((m^2+m)/((m^2+m)-1-0.5*kta^2+an_kta)));                                
    nenner=2*((((m^2+m)-0.5*kta^2+2*an_kla)/(an_kla+1))-(((m^2+m)*(an_kta+1))/((m^2+m)-1-0.5*kta^2+an_kta)));   
    en=zaehler/nenner;                                                                                          
    
    nn=dn*((an_ka-df/dsp*en)/(bn-df/dsp*en));                               
    an(n,1)=-i*nn/(1+i*nn);                                                 
end;
us=an;