%Franziska Minolts, Stand 30.09.2011
%--------------------------------------------------------------------------
function us=test1()                   
df=0.9997;
dsp=2.51;
f=100*10^6;                                                                  
w=2*pi*f;                                                                   
k1=w/1497;                                                                    
k2l=w/5523;                                                              
k2t=w/3382;                                                             
nmax=40;                                                                    %Anzahl der Schleifendurchläufe

m=1;
rad=0;
for z=0.8:0.001:1.5
    rad=rad+exp(z)-1;
    r=rad*10^-6;
    ka=k1*r;                                                                    
    kla=k2l*r;
    kta=k2t*r;
    an=test2(df,dsp,kla,kta,ka,nmax);
    sum1=0;
   
    for n=1:1:nmax
        l_sum1=(2*n-1)*real(an(n,1));                                       %Berechnung der Summe 
        sum1=sum1+l_sum1;
    end;
     streu(m,1)=-sum1*(4*pi/(k1^2));                                        
     durchmesser(m,1)=2*r*10^6;                                             %Durchmesser in µm für graphische Ausgabe

    m=m+1;
end;


loglog(durchmesser,streu,'b');
us=[durchmesser,streu];