clear all;

% constants
f = linspace(1.25,4000,3200);
roh = 2.3e3; % Density                           
E = 2.6e10;  % Young's modulus                         
eta = 0.01;  % loss factor                     
h = 0.10;    % thickness of the plate                        
nue = 0.3;   % Poisson's ratio
lx = 2.8;                                
ly = 2;   
x1 = 1.99;                            
y1 = 1.27;                              
x2 = x1;                              
y2 = y1;
x3 = 1;
y3 = 1;
x4 = 1;
y4 = 1;
x5 = 1;
y5 = 1;
x6 = 1;
y6 = 1;
x7 = 1;
y7 = 1;
x8 = 1;
y8 = 1;
w=(2*pi).*f;                        

% max. frequency
fmax=5000;               

%zeros of the gamma-functions
yi = [7.8532 14.13716];
yj = [4.73004 10.9956];

groessef=size(f);
Yf1z1zsummeges15 = zeros(1,groessef(2));
Yf1z1zsummeges16 = zeros(1,groessef(2));
Yf1z1zsummeges17 = zeros(1,groessef(2));
Yf1z1zsummeges18 = zeros(1,groessef(2));
Yf1z1zsummeges25 = zeros(1,groessef(2));
Yf1z1zsummeges26 = zeros(1,groessef(2));
Yf1z1zsummeges27 = zeros(1,groessef(2));
Yf1z1zsummeges28 = zeros(1,groessef(2));
Yf1z1zsummeges35 = zeros(1,groessef(2));
Yf1z1zsummeges36 = zeros(1,groessef(2));
Yf1z1zsummeges37 = zeros(1,groessef(2));
Yf1z1zsummeges38 = zeros(1,groessef(2));
Yf1z1zsummeges45 = zeros(1,groessef(2));
Yf1z1zsummeges46 = zeros(1,groessef(2));
Yf1z1zsummeges47 = zeros(1,groessef(2));
Yf1z1zsummeges48 = zeros(1,groessef(2));

YMv1x1zsummeges15 = zeros(1,groessef(2));
YMv1x1zsummeges16 = zeros(1,groessef(2));
YMv1x1zsummeges17 = zeros(1,groessef(2));
YMv1x1zsummeges18 = zeros(1,groessef(2));
YMv1x1zsummeges25 = zeros(1,groessef(2));
YMv1x1zsummeges26 = zeros(1,groessef(2));
YMv1x1zsummeges27 = zeros(1,groessef(2));
YMv1x1zsummeges28 = zeros(1,groessef(2));
YMv1x1zsummeges35 = zeros(1,groessef(2));
YMv1x1zsummeges36 = zeros(1,groessef(2));
YMv1x1zsummeges37 = zeros(1,groessef(2));
YMv1x1zsummeges38 = zeros(1,groessef(2));
YMv1x1zsummeges45 = zeros(1,groessef(2));
YMv1x1zsummeges46 = zeros(1,groessef(2));
YMv1x1zsummeges47 = zeros(1,groessef(2));
YMv1x1zsummeges48 = zeros(1,groessef(2));

YMv1y1zsummeges15 = zeros(1,groessef(2));
YMv1y1zsummeges16 = zeros(1,groessef(2));
YMv1y1zsummeges17 = zeros(1,groessef(2));
YMv1y1zsummeges18 = zeros(1,groessef(2));
YMv1y1zsummeges25 = zeros(1,groessef(2));
YMv1y1zsummeges26 = zeros(1,groessef(2));
YMv1y1zsummeges27 = zeros(1,groessef(2));
YMv1y1zsummeges28 = zeros(1,groessef(2));
YMv1y1zsummeges35 = zeros(1,groessef(2));
YMv1y1zsummeges36 = zeros(1,groessef(2));
YMv1y1zsummeges37 = zeros(1,groessef(2));
YMv1y1zsummeges38 = zeros(1,groessef(2));
YMv1y1zsummeges45 = zeros(1,groessef(2));
YMv1y1zsummeges46 = zeros(1,groessef(2));
YMv1y1zsummeges47 = zeros(1,groessef(2));
YMv1y1zsummeges48 = zeros(1,groessef(2));

YfW1z1xsummeges15 = zeros(1,groessef(2));
YfW1z1xsummeges16 = zeros(1,groessef(2));
YfW1z1xsummeges17 = zeros(1,groessef(2));
YfW1z1xsummeges18 = zeros(1,groessef(2));
YfW1z1xsummeges25 = zeros(1,groessef(2));
YfW1z1xsummeges26 = zeros(1,groessef(2));
YfW1z1xsummeges27 = zeros(1,groessef(2));
YfW1z1xsummeges28 = zeros(1,groessef(2));
YfW1z1xsummeges35 = zeros(1,groessef(2));
YfW1z1xsummeges36 = zeros(1,groessef(2));
YfW1z1xsummeges37 = zeros(1,groessef(2));
YfW1z1xsummeges38 = zeros(1,groessef(2));
YfW1z1xsummeges45 = zeros(1,groessef(2));
YfW1z1xsummeges46 = zeros(1,groessef(2));
YfW1z1xsummeges47 = zeros(1,groessef(2));
YfW1z1xsummeges48 = zeros(1,groessef(2));

YM1x1xsummeges15 = zeros(1,groessef(2));
YM1x1xsummeges16 = zeros(1,groessef(2));
YM1x1xsummeges17 = zeros(1,groessef(2));
YM1x1xsummeges18 = zeros(1,groessef(2));
YM1x1xsummeges25 = zeros(1,groessef(2));
YM1x1xsummeges26 = zeros(1,groessef(2));
YM1x1xsummeges27 = zeros(1,groessef(2));
YM1x1xsummeges28 = zeros(1,groessef(2));
YM1x1xsummeges35 = zeros(1,groessef(2));
YM1x1xsummeges36 = zeros(1,groessef(2));
YM1x1xsummeges37 = zeros(1,groessef(2));
YM1x1xsummeges38 = zeros(1,groessef(2));
YM1x1xsummeges45 = zeros(1,groessef(2));
YM1x1xsummeges46 = zeros(1,groessef(2));
YM1x1xsummeges47 = zeros(1,groessef(2));
YM1x1xsummeges48 = zeros(1,groessef(2));

YM1y1xsummeges15 = zeros(1,groessef(2));
YM1y1xsummeges16 = zeros(1,groessef(2));
YM1y1xsummeges17 = zeros(1,groessef(2));
YM1y1xsummeges18 = zeros(1,groessef(2));
YM1y1xsummeges25 = zeros(1,groessef(2));
YM1y1xsummeges26 = zeros(1,groessef(2));
YM1y1xsummeges27 = zeros(1,groessef(2));
YM1y1xsummeges28 = zeros(1,groessef(2));
YM1y1xsummeges35 = zeros(1,groessef(2));
YM1y1xsummeges36 = zeros(1,groessef(2));
YM1y1xsummeges37 = zeros(1,groessef(2));
YM1y1xsummeges38 = zeros(1,groessef(2));
YM1y1xsummeges45 = zeros(1,groessef(2));
YM1y1xsummeges46 = zeros(1,groessef(2));
YM1y1xsummeges47 = zeros(1,groessef(2));
YM1y1xsummeges48 = zeros(1,groessef(2));

YfW1z1ysummeges15 = zeros(1,groessef(2));
YfW1z1ysummeges16 = zeros(1,groessef(2));
YfW1z1ysummeges17 = zeros(1,groessef(2));
YfW1z1ysummeges18 = zeros(1,groessef(2));
YfW1z1ysummeges25 = zeros(1,groessef(2));
YfW1z1ysummeges26 = zeros(1,groessef(2));
YfW1z1ysummeges27 = zeros(1,groessef(2));
YfW1z1ysummeges28 = zeros(1,groessef(2));
YfW1z1ysummeges35 = zeros(1,groessef(2));
YfW1z1ysummeges36 = zeros(1,groessef(2));
YfW1z1ysummeges37 = zeros(1,groessef(2));
YfW1z1ysummeges38 = zeros(1,groessef(2));
YfW1z1ysummeges45 = zeros(1,groessef(2));
YfW1z1ysummeges46 = zeros(1,groessef(2));
YfW1z1ysummeges47 = zeros(1,groessef(2));
YfW1z1ysummeges48 = zeros(1,groessef(2));

YM1x1ysummeges15 = zeros(1,groessef(2));
YM1x1ysummeges16 = zeros(1,groessef(2));
YM1x1ysummeges17 = zeros(1,groessef(2));
YM1x1ysummeges18 = zeros(1,groessef(2));
YM1x1ysummeges25 = zeros(1,groessef(2));
YM1x1ysummeges26 = zeros(1,groessef(2));
YM1x1ysummeges27 = zeros(1,groessef(2));
YM1x1ysummeges28 = zeros(1,groessef(2));
YM1x1ysummeges35 = zeros(1,groessef(2));
YM1x1ysummeges36 = zeros(1,groessef(2));
YM1x1ysummeges37 = zeros(1,groessef(2));
YM1x1ysummeges38 = zeros(1,groessef(2));
YM1x1ysummeges45 = zeros(1,groessef(2));
YM1x1ysummeges46 = zeros(1,groessef(2));
YM1x1ysummeges47 = zeros(1,groessef(2));
YM1x1ysummeges48 = zeros(1,groessef(2));

YM1y1ysummeges15 = zeros(1,groessef(2));
YM1y1ysummeges16 = zeros(1,groessef(2));
YM1y1ysummeges17 = zeros(1,groessef(2));
YM1y1ysummeges18 = zeros(1,groessef(2));
YM1y1ysummeges25 = zeros(1,groessef(2));
YM1y1ysummeges26 = zeros(1,groessef(2));
YM1y1ysummeges27 = zeros(1,groessef(2));
YM1y1ysummeges28 = zeros(1,groessef(2));
YM1y1ysummeges35 = zeros(1,groessef(2));
YM1y1ysummeges36 = zeros(1,groessef(2));
YM1y1ysummeges37 = zeros(1,groessef(2));
YM1y1ysummeges38 = zeros(1,groessef(2));
YM1y1ysummeges45 = zeros(1,groessef(2));
YM1y1ysummeges46 = zeros(1,groessef(2));
YM1y1ysummeges47 = zeros(1,groessef(2));
YM1y1ysummeges48 = zeros(1,groessef(2));

Yf1z1zsumme15 = zeros(1,groessef(2));
Yf1z1zsumme16 = zeros(1,groessef(2));
Yf1z1zsumme17 = zeros(1,groessef(2));
Yf1z1zsumme18 = zeros(1,groessef(2));
Yf1z1zsumme25 = zeros(1,groessef(2));
Yf1z1zsumme26 = zeros(1,groessef(2));
Yf1z1zsumme27 = zeros(1,groessef(2));
Yf1z1zsumme28 = zeros(1,groessef(2));
Yf1z1zsumme35 = zeros(1,groessef(2));
Yf1z1zsumme36 = zeros(1,groessef(2));
Yf1z1zsumme37 = zeros(1,groessef(2));
Yf1z1zsumme38 = zeros(1,groessef(2));
Yf1z1zsumme45 = zeros(1,groessef(2));
Yf1z1zsumme46 = zeros(1,groessef(2));
Yf1z1zsumme47 = zeros(1,groessef(2));
Yf1z1zsumme48 = zeros(1,groessef(2));

YMv1x1zsumme15 = zeros(1,groessef(2));
YMv1x1zsumme16 = zeros(1,groessef(2));
YMv1x1zsumme17 = zeros(1,groessef(2));
YMv1x1zsumme18 = zeros(1,groessef(2));
YMv1x1zsumme25 = zeros(1,groessef(2));
YMv1x1zsumme26 = zeros(1,groessef(2));
YMv1x1zsumme27 = zeros(1,groessef(2));
YMv1x1zsumme28 = zeros(1,groessef(2));
YMv1x1zsumme35 = zeros(1,groessef(2));
YMv1x1zsumme36 = zeros(1,groessef(2));
YMv1x1zsumme37 = zeros(1,groessef(2));
YMv1x1zsumme38 = zeros(1,groessef(2));
YMv1x1zsumme45 = zeros(1,groessef(2));
YMv1x1zsumme46 = zeros(1,groessef(2));
YMv1x1zsumme47 = zeros(1,groessef(2));
YMv1x1zsumme48 = zeros(1,groessef(2));
        
YMv1y1zsumme15 = zeros(1,groessef(2));
YMv1y1zsumme16 = zeros(1,groessef(2));
YMv1y1zsumme17 = zeros(1,groessef(2));
YMv1y1zsumme18 = zeros(1,groessef(2));
YMv1y1zsumme25 = zeros(1,groessef(2));
YMv1y1zsumme26 = zeros(1,groessef(2));
YMv1y1zsumme27 = zeros(1,groessef(2));
YMv1y1zsumme28 = zeros(1,groessef(2));
YMv1y1zsumme35 = zeros(1,groessef(2));
YMv1y1zsumme36 = zeros(1,groessef(2));
YMv1y1zsumme37 = zeros(1,groessef(2));
YMv1y1zsumme38 = zeros(1,groessef(2));
YMv1y1zsumme45 = zeros(1,groessef(2));
YMv1y1zsumme46 = zeros(1,groessef(2));
YMv1y1zsumme47 = zeros(1,groessef(2));
YMv1y1zsumme48 = zeros(1,groessef(2));
             
YfW1z1xsumme15 = zeros(1,groessef(2));
YfW1z1xsumme16 = zeros(1,groessef(2));
YfW1z1xsumme17 = zeros(1,groessef(2));
YfW1z1xsumme18 = zeros(1,groessef(2));
YfW1z1xsumme25 = zeros(1,groessef(2));
YfW1z1xsumme26 = zeros(1,groessef(2));
YfW1z1xsumme27 = zeros(1,groessef(2));
YfW1z1xsumme28 = zeros(1,groessef(2));
YfW1z1xsumme35 = zeros(1,groessef(2));
YfW1z1xsumme36 = zeros(1,groessef(2));  
YfW1z1xsumme37 = zeros(1,groessef(2));
YfW1z1xsumme38 = zeros(1,groessef(2)); 
YfW1z1xsumme45 = zeros(1,groessef(2));
YfW1z1xsumme46 = zeros(1,groessef(2));      
YfW1z1xsumme47 = zeros(1,groessef(2));
YfW1z1xsumme48 = zeros(1,groessef(2));
             
YM1x1xsumme15 = zeros(1,groessef(2));
YM1x1xsumme16 = zeros(1,groessef(2));
YM1x1xsumme17 = zeros(1,groessef(2));
YM1x1xsumme18 = zeros(1,groessef(2));
YM1x1xsumme25 = zeros(1,groessef(2));
YM1x1xsumme26 = zeros(1,groessef(2));
YM1x1xsumme27 = zeros(1,groessef(2));
YM1x1xsumme28 = zeros(1,groessef(2));
YM1x1xsumme35 = zeros(1,groessef(2));
YM1x1xsumme36 = zeros(1,groessef(2));
YM1x1xsumme37 = zeros(1,groessef(2));
YM1x1xsumme38 = zeros(1,groessef(2));
YM1x1xsumme45 = zeros(1,groessef(2));
YM1x1xsumme46 = zeros(1,groessef(2));
YM1x1xsumme47 = zeros(1,groessef(2));
YM1x1xsumme48 = zeros(1,groessef(2));
       
YM1y1xsumme15 = zeros(1,groessef(2));
YM1y1xsumme16 = zeros(1,groessef(2));
YM1y1xsumme17 = zeros(1,groessef(2));
YM1y1xsumme18 = zeros(1,groessef(2));
YM1y1xsumme25 = zeros(1,groessef(2));
YM1y1xsumme26 = zeros(1,groessef(2));
YM1y1xsumme27 = zeros(1,groessef(2));
YM1y1xsumme28 = zeros(1,groessef(2));
YM1y1xsumme35 = zeros(1,groessef(2));
YM1y1xsumme36 = zeros(1,groessef(2));
YM1y1xsumme37 = zeros(1,groessef(2));
YM1y1xsumme38 = zeros(1,groessef(2));
YM1y1xsumme45 = zeros(1,groessef(2));
YM1y1xsumme46 = zeros(1,groessef(2));
YM1y1xsumme47 = zeros(1,groessef(2));
YM1y1xsumme48 = zeros(1,groessef(2));
             
YfW1z1ysumme15 = zeros(1,groessef(2));
YfW1z1ysumme16 = zeros(1,groessef(2)); 
YfW1z1ysumme17 = zeros(1,groessef(2)); 
YfW1z1ysumme18 = zeros(1,groessef(2));
YfW1z1ysumme25 = zeros(1,groessef(2));
YfW1z1ysumme26 = zeros(1,groessef(2));
YfW1z1ysumme27 = zeros(1,groessef(2));
YfW1z1ysumme28 = zeros(1,groessef(2)); 
YfW1z1ysumme35 = zeros(1,groessef(2));
YfW1z1ysumme36 = zeros(1,groessef(2));
YfW1z1ysumme37 = zeros(1,groessef(2)); 
YfW1z1ysumme38 = zeros(1,groessef(2)); 
YfW1z1ysumme45 = zeros(1,groessef(2));
YfW1z1ysumme46 = zeros(1,groessef(2));  
YfW1z1ysumme47 = zeros(1,groessef(2));
YfW1z1ysumme48 = zeros(1,groessef(2));
       
YM1x1ysumme15 = zeros(1,groessef(2));
YM1x1ysumme16 = zeros(1,groessef(2));
YM1x1ysumme17 = zeros(1,groessef(2));
YM1x1ysumme18 = zeros(1,groessef(2));
YM1x1ysumme25 = zeros(1,groessef(2));
YM1x1ysumme26 = zeros(1,groessef(2));
YM1x1ysumme27 = zeros(1,groessef(2));
YM1x1ysumme28 = zeros(1,groessef(2));
YM1x1ysumme35 = zeros(1,groessef(2));
YM1x1ysumme36 = zeros(1,groessef(2));
YM1x1ysumme37 = zeros(1,groessef(2));
YM1x1ysumme38 = zeros(1,groessef(2));
YM1x1ysumme45 = zeros(1,groessef(2));
YM1x1ysumme46 = zeros(1,groessef(2));
YM1x1ysumme47 = zeros(1,groessef(2));
YM1x1ysumme48 = zeros(1,groessef(2));
             
YM1y1ysumme15 = zeros(1,groessef(2));
YM1y1ysumme16 = zeros(1,groessef(2));
YM1y1ysumme17 = zeros(1,groessef(2));
YM1y1ysumme18 = zeros(1,groessef(2));
YM1y1ysumme25 = zeros(1,groessef(2));
YM1y1ysumme26 = zeros(1,groessef(2));
YM1y1ysumme27 = zeros(1,groessef(2));
YM1y1ysumme28 = zeros(1,groessef(2));
YM1y1ysumme35 = zeros(1,groessef(2));
YM1y1ysumme36 = zeros(1,groessef(2));
YM1y1ysumme37 = zeros(1,groessef(2));
YM1y1ysumme38 = zeros(1,groessef(2));
YM1y1ysumme45 = zeros(1,groessef(2));
YM1y1ysumme46 = zeros(1,groessef(2));
YM1y1ysumme47 = zeros(1,groessef(2));
YM1y1ysumme48 = zeros(1,groessef(2));

m = 1;
n = 0;
z = 0; %Zähler 
while (1)
    z = z + 1;
    if m == 0
        Gx = 0;
        Hx = 0;
        Jx = 0;
    elseif m == 1
        Gx = 1.506;
        Hx = 1.248;
        Jx = 5.017;
    else
        Gx = m+0.5;
        Hx = (m+0.5)^2*(1-4/((2*m+1)*pi));
        Jx = (m+0.5)^2*(1+12/((2*m+1)*pi));
    end
    if n == 0
        Gy = 0;
        Hy = 0;
        Jy = 0;
    elseif n == 1
        Gy = 1.506;
        Hy = 1.248;
        Jy = 5.017;
    else
        Gy = n+0.5;
        Hy = (n+0.5)^2*(1-4/((2*n+1)*pi));
        Jy = (n+0.5)^2*(1+12/((2*n+1)*pi));
    end
    qmn = sqrt(Gx^4+Gy^4*(lx/ly)^4+2*(lx/ly)^2*(nue*Hx*Hy+(1-nue)*Jx*Jy)); 
    wmnr = sqrt(E*h^2/(12*roh*(1-nue^2)))*(pi/lx)^2*qmn; 
    wmn = wmnr*sqrt(1+j*eta); 
    fmn = wmn/(2*pi);                                     
     
       if m==0
           phix1 = 1;
           phix1st = 0; 
           phix2 = 1;
           phix2st = 0;
           phix3 = 1;
           phix3st = 0; 
           phix4 = 1;
           phix4st = 0;
           phix5 = 1;
           phix5st = 0; 
           phix6 = 1;
           phix6st = 0;
           phix7 = 1;
           phix7st = 0; 
           phix8 = 1;
           phix8st = 0;
           
       else
           if m/2 == round(m/2); %Abfrage: m gerade Zahl
                ii = m/2;
                if ii > 2;
                   yi(ii) = (4*ii+1)*pi/2;
                end                
                kn = sin(0.5*yi(ii))/sinh(0.5*yi(ii));
                phix1 = sqrt(2)*(sin(yi(ii)*(x1/lx-0.5))+kn*sinh(yi(ii)*(x1/lx-0.5))); %phi an der Stelle x1
                phix1st = sqrt(2)*(cos(yi(ii)*(x1/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x1/lx-0.5))*yi(ii)/lx);
                phix2 = sqrt(2)*(sin(yi(ii)*(x2/lx-0.5))+kn*sinh(yi(ii)*(x2/lx-0.5))); %phi an der Stelle x2
                phix2st = sqrt(2)*(cos(yi(ii)*(x2/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x2/lx-0.5))*yi(ii)/lx);
                phix3 = sqrt(2)*(sin(yi(ii)*(x3/lx-0.5))+kn*sinh(yi(ii)*(x3/lx-0.5))); %phi an der Stelle x1
                phix3st = sqrt(2)*(cos(yi(ii)*(x3/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x3/lx-0.5))*yi(ii)/lx);
                phix4 = sqrt(2)*(sin(yi(ii)*(x4/lx-0.5))+kn*sinh(yi(ii)*(x4/lx-0.5))); %phi an der Stelle x2
                phix4st = sqrt(2)*(cos(yi(ii)*(x4/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x4/lx-0.5))*yi(ii)/lx);
                phix5 = sqrt(2)*(sin(yi(ii)*(x5/lx-0.5))+kn*sinh(yi(ii)*(x5/lx-0.5))); %phi an der Stelle x1
                phix5st = sqrt(2)*(cos(yi(ii)*(x5/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x5/lx-0.5))*yi(ii)/lx);
                phix6 = sqrt(2)*(sin(yi(ii)*(x6/lx-0.5))+kn*sinh(yi(ii)*(x6/lx-0.5))); %phi an der Stelle x2
                phix6st = sqrt(2)*(cos(yi(ii)*(x6/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x6/lx-0.5))*yi(ii)/lx);
                phix7 = sqrt(2)*(sin(yi(ii)*(x7/lx-0.5))+kn*sinh(yi(ii)*(x7/lx-0.5))); %phi an der Stelle x1
                phix7st = sqrt(2)*(cos(yi(ii)*(x7/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x7/lx-0.5))*yi(ii)/lx);
                phix8 = sqrt(2)*(sin(yi(ii)*(x8/lx-0.5))+kn*sinh(yi(ii)*(x8/lx-0.5))); %phi an der Stelle x2
                phix8st = sqrt(2)*(cos(yi(ii)*(x8/lx-0.5))*yi(ii)/lx+kn*cosh(yi(ii)*(x8/lx-0.5))*yi(ii)/lx);
                
           else % m ungerade Zahl
                jj = (m+1)/2;
                if jj > 2;
                   yj(jj) = (4*jj-1)*pi/2;
                end
                kn = -sin(0.5*yj(jj))/sinh(0.5*yj(jj));
                phix1 = sqrt(2)*(cos(yj(jj)*(x1/lx-0.5))+kn*cosh(yj(jj)*(x1/lx-0.5))); 
                phix1st = sqrt(2)*(-sin(yj(jj)*(x1/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x1/lx-0.5))*yj(jj)/lx);
                phix2 = sqrt(2)*(cos(yj(jj)*(x2/lx-0.5))+kn*cosh(yj(jj)*(x2/lx-0.5)));
                phix2st = sqrt(2)*(-sin(yj(jj)*(x2/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x2/lx-0.5))*yj(jj)/lx);
                phix3 = sqrt(2)*(cos(yj(jj)*(x3/lx-0.5))+kn*cosh(yj(jj)*(x3/lx-0.5))); 
                phix3st = sqrt(2)*(-sin(yj(jj)*(x3/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x3/lx-0.5))*yj(jj)/lx);
                phix4 = sqrt(2)*(cos(yj(jj)*(x4/lx-0.5))+kn*cosh(yj(jj)*(x4/lx-0.5)));
                phix4st = sqrt(2)*(-sin(yj(jj)*(x4/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x4/lx-0.5))*yj(jj)/lx);
                phix5 = sqrt(2)*(cos(yj(jj)*(x5/lx-0.5))+kn*cosh(yj(jj)*(x5/lx-0.5))); 
                phix5st = sqrt(2)*(-sin(yj(jj)*(x5/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x5/lx-0.5))*yj(jj)/lx);
                phix6 = sqrt(2)*(cos(yj(jj)*(x6/lx-0.5))+kn*cosh(yj(jj)*(x6/lx-0.5)));
                phix6st = sqrt(2)*(-sin(yj(jj)*(x6/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x6/lx-0.5))*yj(jj)/lx);
                phix7 = sqrt(2)*(cos(yj(jj)*(x7/lx-0.5))+kn*cosh(yj(jj)*(x7/lx-0.5))); 
                phix7st = sqrt(2)*(-sin(yj(jj)*(x7/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x7/lx-0.5))*yj(jj)/lx);
                phix8 = sqrt(2)*(cos(yj(jj)*(x8/lx-0.5))+kn*cosh(yj(jj)*(x8/lx-0.5)));
                phix8st = sqrt(2)*(-sin(yj(jj)*(x8/lx-0.5))*yj(jj)/lx+kn*sinh(yj(jj)*(x8/lx-0.5))*yj(jj)/lx);
                
           end % der if Abfrage m gerade Zahl
        end
        if n == 0
            phiy1 = 1;
            phiy1st = 0;
            phiy2 = 1;
            phiy2st = 0;
            phiy3 = 1;
            phiy3st = 0;
            phiy4 = 1;
            phiy4st = 0;
            phiy5 = 1;
            phiy5st = 0;
            phiy6 = 1;
            phiy6st = 0;
            phiy7 = 1;
            phiy7st = 0;
            phiy8 = 1;
            phiy8st = 0;
        else
            if n/2 == round(n/2); %Abfrage: n gerade Zahl
                ii = n/2;
                if ii > 2;
                   yi(ii) = (4*ii+1)*pi/2;
                end                
                kn = sin(0.5*yi(ii))/sinh(0.5*yi(ii));
                phiy1 = sqrt(2)*(sin(yi(ii)*(y1/ly-0.5))+kn*sinh(yi(ii)*(y1/ly-0.5))); 
                phiy1st = sqrt(2)*(cos(yi(ii)*(y1/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y1/ly-0.5))*yi(ii)/ly);
                phiy2 = sqrt(2)*(sin(yi(ii)*(y2/ly-0.5))+kn*sinh(yi(ii)*(y2/ly-0.5))); 
                phiy2st = sqrt(2)*(cos(yi(ii)*(y2/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y2/ly-0.5))*yi(ii)/ly);
                phiy3 = sqrt(2)*(sin(yi(ii)*(y3/ly-0.5))+kn*sinh(yi(ii)*(y3/ly-0.5))); 
                phiy3st = sqrt(2)*(cos(yi(ii)*(y3/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y3/ly-0.5))*yi(ii)/ly);
                phiy4 = sqrt(2)*(sin(yi(ii)*(y4/ly-0.5))+kn*sinh(yi(ii)*(y4/ly-0.5))); 
                phiy4st = sqrt(2)*(cos(yi(ii)*(y4/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y4/ly-0.5))*yi(ii)/ly);
                phiy5 = sqrt(2)*(sin(yi(ii)*(y5/ly-0.5))+kn*sinh(yi(ii)*(y5/ly-0.5))); 
                phiy5st = sqrt(2)*(cos(yi(ii)*(y5/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y5/ly-0.5))*yi(ii)/ly);
                phiy6 = sqrt(2)*(sin(yi(ii)*(y6/ly-0.5))+kn*sinh(yi(ii)*(y6/ly-0.5))); 
                phiy6st = sqrt(2)*(cos(yi(ii)*(y6/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y6/ly-0.5))*yi(ii)/ly);
                phiy7 = sqrt(2)*(sin(yi(ii)*(y7/ly-0.5))+kn*sinh(yi(ii)*(y7/ly-0.5))); 
                phiy7st = sqrt(2)*(cos(yi(ii)*(y7/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y7/ly-0.5))*yi(ii)/ly);
                phiy8 = sqrt(2)*(sin(yi(ii)*(y8/ly-0.5))+kn*sinh(yi(ii)*(y8/ly-0.5))); 
                phiy8st = sqrt(2)*(cos(yi(ii)*(y8/ly-0.5))*yi(ii)/ly+kn*cosh(yi(ii)*(y8/ly-0.5))*yi(ii)/ly);
                
            else % n ungerade Zahl
                jj = (n+1)/2;
                if jj > 2;
                   yj(jj) = (4*jj-1)*pi/2;
                end
                kn = -sin(0.5*yj(jj))/sinh(0.5*yj(jj));
                phiy1 = sqrt(2)*(cos(yj(jj)*(y1/ly-0.5))+kn*cosh(yj(jj)*(y1/ly-0.5))); 
                phiy1st = sqrt(2)*(-sin(yj(jj)*(y1/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y1/ly-0.5))*yj(jj)/ly);
                phiy2 = sqrt(2)*(cos(yj(jj)*(y2/ly-0.5))+kn*cosh(yj(jj)*(y2/ly-0.5))); 
                phiy2st = sqrt(2)*(-sin(yj(jj)*(y2/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y2/ly-0.5))*yj(jj)/ly);
                phiy3 = sqrt(2)*(cos(yj(jj)*(y3/ly-0.5))+kn*cosh(yj(jj)*(y3/ly-0.5))); 
                phiy3st = sqrt(2)*(-sin(yj(jj)*(y3/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y3/ly-0.5))*yj(jj)/ly);
                phiy4 = sqrt(2)*(cos(yj(jj)*(y4/ly-0.5))+kn*cosh(yj(jj)*(y4/ly-0.5))); 
                phiy4st = sqrt(2)*(-sin(yj(jj)*(y4/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y4/ly-0.5))*yj(jj)/ly);
                phiy5 = sqrt(2)*(cos(yj(jj)*(y5/ly-0.5))+kn*cosh(yj(jj)*(y5/ly-0.5))); 
                phiy5st = sqrt(2)*(-sin(yj(jj)*(y5/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y5/ly-0.5))*yj(jj)/ly);
                phiy6 = sqrt(2)*(cos(yj(jj)*(y6/ly-0.5))+kn*cosh(yj(jj)*(y6/ly-0.5))); 
                phiy6st = sqrt(2)*(-sin(yj(jj)*(y6/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y6/ly-0.5))*yj(jj)/ly);
                phiy7 = sqrt(2)*(cos(yj(jj)*(y7/ly-0.5))+kn*cosh(yj(jj)*(y7/ly-0.5))); 
                phiy7st = sqrt(2)*(-sin(yj(jj)*(y7/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y7/ly-0.5))*yj(jj)/ly);
                phiy8 = sqrt(2)*(cos(yj(jj)*(y8/ly-0.5))+kn*cosh(yj(jj)*(y8/ly-0.5))); 
                phiy8st = sqrt(2)*(-sin(yj(jj)*(y8/ly-0.5))*yj(jj)/ly+kn*sinh(yj(jj)*(y8/ly-0.5))*yj(jj)/ly);
                
             end % der if Abfrage n gerade Zahl
        end

       psi1 = phix1*phiy1; % psi an der Stelle 1
       psi2 = phix2*phiy2; % psi an der Stelle 2
       psi3 = phix3*phiy3; % psi an der Stelle 1
       psi4 = phix4*phiy4; % psi an der Stelle 2
       psi5 = phix5*phiy5; % psi an der Stelle 1
       psi6 = phix6*phiy6; % psi an der Stelle 2
       psi7 = phix7*phiy7; % psi an der Stelle 1
       psi8 = phix8*phiy8; % psi an der Stelle 2
       
       psiX1 = phix1*phiy1st; % psi x (m,n) an der Stelle 1 
       psiX2 = phix2*phiy2st; % psi x (m,n) an der Stelle 2
       psiX3 = phix3*phiy3st; % psi x (m,n) an der Stelle 1 
       psiX4 = phix4*phiy4st; % psi x (m,n) an der Stelle 2
       psiX5 = phix5*phiy5st; % psi x (m,n) an der Stelle 1 
       psiX6 = phix6*phiy6st; % psi x (m,n) an der Stelle 2
       psiX7 = phix7*phiy7st; % psi x (m,n) an der Stelle 1 
       psiX8 = phix8*phiy8st; % psi x (m,n) an der Stelle 2
              
       psiY1 = -phix1st*phiy1; % psi y (m,n) an der Stelle 1 
       psiY2 = -phix2st*phiy2; % psi y (m,n) an der Stelle 2
       psiY3 = -phix3st*phiy3; % psi y (m,n) an der Stelle 1 
       psiY4 = -phix4st*phiy4; % psi y (m,n) an der Stelle 2
       psiY5 = -phix5st*phiy5; % psi y (m,n) an der Stelle 1 
       psiY6 = -phix6st*phiy6; % psi y (m,n) an der Stelle 2
       psiY7 = -phix7st*phiy7; % psi y (m,n) an der Stelle 1 
       psiY8 = -phix8st*phiy8; % psi y (m,n) an der Stelle 2
                   
       Yf1z1zsumme15(z,:) = psi1.*psi5./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme16(z,:) = psi1.*psi6./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme17(z,:) = psi1.*psi7./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme18(z,:) = psi1.*psi8./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme25(z,:) = psi2.*psi5./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme26(z,:) = psi2.*psi6./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme27(z,:) = psi2.*psi7./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme28(z,:) = psi2.*psi8./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme35(z,:) = psi3.*psi5./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme36(z,:) = psi3.*psi6./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme37(z,:) = psi3.*psi7./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme38(z,:) = psi3.*psi8./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme45(z,:) = psi4.*psi5./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme46(z,:) = psi4.*psi6./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme47(z,:) = psi4.*psi7./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       Yf1z1zsumme48(z,:) = psi4.*psi8./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
                   
       YMv1x1zsumme15(z,:) = psi5.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme16(z,:) = psi6.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme17(z,:) = psi7.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme18(z,:) = psi8.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme25(z,:) = psi5.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme26(z,:) = psi6.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme27(z,:) = psi7.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme28(z,:) = psi8.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme35(z,:) = psi5.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme36(z,:) = psi6.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme37(z,:) = psi7.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme38(z,:) = psi8.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme45(z,:) = psi5.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme46(z,:) = psi6.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme47(z,:) = psi7.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1x1zsumme48(z,:) = psi8.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
        
       YMv1y1zsumme15(z,:) = psi5.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme16(z,:) = psi6.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme17(z,:) = psi7.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme18(z,:) = psi8.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme25(z,:) = psi5.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme26(z,:) = psi6.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme27(z,:) = psi7.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme28(z,:) = psi8.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme35(z,:) = psi5.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme36(z,:) = psi6.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme37(z,:) = psi7.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme38(z,:) = psi8.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme45(z,:) = psi5.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme46(z,:) = psi6.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme47(z,:) = psi7.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YMv1y1zsumme48(z,:) = psi8.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
             
       YfW1z1xsumme15(z,:) = psiX5.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme16(z,:) = psiX6.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));      
       YfW1z1xsumme17(z,:) = psiX7.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme18(z,:) = psiX8.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));    
       YfW1z1xsumme25(z,:) = psiX5.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme26(z,:) = psiX6.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));      
       YfW1z1xsumme27(z,:) = psiX7.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme28(z,:) = psiX8.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));    
       YfW1z1xsumme35(z,:) = psiX5.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme36(z,:) = psiX6.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));      
       YfW1z1xsumme37(z,:) = psiX7.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme38(z,:) = psiX8.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));    
       YfW1z1xsumme45(z,:) = psiX5.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme46(z,:) = psiX6.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));      
       YfW1z1xsumme47(z,:) = psiX7.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1xsumme48(z,:) = psiX8.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
             
       YM1x1xsumme15(z,:) = psiX5.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme16(z,:) = psiX6.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme17(z,:) = psiX7.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme18(z,:) = psiX8.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme25(z,:) = psiX5.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme26(z,:) = psiX6.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme27(z,:) = psiX7.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme28(z,:) = psiX8.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme35(z,:) = psiX5.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme36(z,:) = psiX6.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme37(z,:) = psiX7.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme38(z,:) = psiX8.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme45(z,:) = psiX5.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme46(z,:) = psiX6.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme47(z,:) = psiX7.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1xsumme48(z,:) = psiX8.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       
       YM1y1xsumme15(z,:) = psiX5.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme16(z,:) = psiX6.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme17(z,:) = psiX7.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme18(z,:) = psiX8.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme25(z,:) = psiX5.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme26(z,:) = psiX6.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme27(z,:) = psiX7.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme28(z,:) = psiX8.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme35(z,:) = psiX5.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme36(z,:) = psiX6.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme37(z,:) = psiX7.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme38(z,:) = psiX8.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme45(z,:) = psiX5.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme46(z,:) = psiX6.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme47(z,:) = psiX7.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1xsumme48(z,:) = psiX8.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
             
       YfW1z1ysumme15(z,:) = psiY5.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme16(z,:) = psiY6.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme17(z,:) = psiY7.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme18(z,:) = psiY8.*psi1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme25(z,:) = psiY5.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme26(z,:) = psiY6.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme27(z,:) = psiY7.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme28(z,:) = psiY8.*psi2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme35(z,:) = psiY5.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme36(z,:) = psiY6.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme37(z,:) = psiY7.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme38(z,:) = psiY8.*psi3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme45(z,:) = psiY5.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme46(z,:) = psiY6.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme47(z,:) = psiY7.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YfW1z1ysumme48(z,:) = psiY8.*psi4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       
       YM1x1ysumme15(z,:) = psiY2.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme16(z,:) = psiY2.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme17(z,:) = psiY2.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme18(z,:) = psiY2.*psiX1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme25(z,:) = psiY2.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme26(z,:) = psiY2.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme27(z,:) = psiY2.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme28(z,:) = psiY2.*psiX2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme35(z,:) = psiY2.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme36(z,:) = psiY2.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme37(z,:) = psiY2.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme38(z,:) = psiY2.*psiX3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme45(z,:) = psiY2.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme46(z,:) = psiY2.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme47(z,:) = psiY2.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1x1ysumme48(z,:) = psiY2.*psiX4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
             
       YM1y1ysumme15(z,:) = psiY5.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme16(z,:) = psiY6.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme17(z,:) = psiY7.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme18(z,:) = psiY8.*psiY1./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme25(z,:) = psiY5.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme26(z,:) = psiY6.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme27(z,:) = psiY7.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme28(z,:) = psiY8.*psiY2./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme35(z,:) = psiY5.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme36(z,:) = psiY6.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme37(z,:) = psiY7.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme38(z,:) = psiY8.*psiY3./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme45(z,:) = psiY5.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme46(z,:) = psiY6.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme47(z,:) = psiY7.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       YM1y1ysumme48(z,:) = psiY8.*psiY4./(roh.*h.*lx.*ly.*(wmn.^2-w.^2));
       
       
       
       Yf1z1zsummeges15(1,:) = Yf1z1zsumme15(z,:)+Yf1z1zsummeges15(1,:);
       Yf1z1zsummeges16(1,:) = Yf1z1zsumme16(z,:)+Yf1z1zsummeges16(1,:);
       Yf1z1zsummeges17(1,:) = Yf1z1zsumme17(z,:)+Yf1z1zsummeges17(1,:);
       Yf1z1zsummeges18(1,:) = Yf1z1zsumme18(z,:)+Yf1z1zsummeges18(1,:);
       Yf1z1zsummeges25(1,:) = Yf1z1zsumme25(z,:)+Yf1z1zsummeges25(1,:);
       Yf1z1zsummeges26(1,:) = Yf1z1zsumme26(z,:)+Yf1z1zsummeges26(1,:);
       Yf1z1zsummeges27(1,:) = Yf1z1zsumme27(z,:)+Yf1z1zsummeges27(1,:);
       Yf1z1zsummeges28(1,:) = Yf1z1zsumme28(z,:)+Yf1z1zsummeges28(1,:);
       Yf1z1zsummeges35(1,:) = Yf1z1zsumme35(z,:)+Yf1z1zsummeges35(1,:);
       Yf1z1zsummeges36(1,:) = Yf1z1zsumme36(z,:)+Yf1z1zsummeges36(1,:);
       Yf1z1zsummeges37(1,:) = Yf1z1zsumme37(z,:)+Yf1z1zsummeges37(1,:);
       Yf1z1zsummeges38(1,:) = Yf1z1zsumme38(z,:)+Yf1z1zsummeges38(1,:);
       Yf1z1zsummeges45(1,:) = Yf1z1zsumme45(z,:)+Yf1z1zsummeges45(1,:);
       Yf1z1zsummeges46(1,:) = Yf1z1zsumme46(z,:)+Yf1z1zsummeges46(1,:);
       Yf1z1zsummeges47(1,:) = Yf1z1zsumme47(z,:)+Yf1z1zsummeges47(1,:);
       Yf1z1zsummeges48(1,:) = Yf1z1zsumme48(z,:)+Yf1z1zsummeges48(1,:);
             
       YMv1x1zsummeges15(1,:) = YMv1x1zsumme15(z,:)+YMv1x1zsummeges15(1,:);
       YMv1x1zsummeges16(1,:) = YMv1x1zsumme16(z,:)+YMv1x1zsummeges16(1,:);
       YMv1x1zsummeges17(1,:) = YMv1x1zsumme17(z,:)+YMv1x1zsummeges17(1,:);
       YMv1x1zsummeges18(1,:) = YMv1x1zsumme18(z,:)+YMv1x1zsummeges18(1,:);
       YMv1x1zsummeges25(1,:) = YMv1x1zsumme25(z,:)+YMv1x1zsummeges25(1,:);
       YMv1x1zsummeges26(1,:) = YMv1x1zsumme26(z,:)+YMv1x1zsummeges26(1,:);
       YMv1x1zsummeges27(1,:) = YMv1x1zsumme27(z,:)+YMv1x1zsummeges27(1,:);
       YMv1x1zsummeges28(1,:) = YMv1x1zsumme28(z,:)+YMv1x1zsummeges28(1,:);
       YMv1x1zsummeges35(1,:) = YMv1x1zsumme35(z,:)+YMv1x1zsummeges35(1,:);
       YMv1x1zsummeges36(1,:) = YMv1x1zsumme36(z,:)+YMv1x1zsummeges36(1,:);
       YMv1x1zsummeges37(1,:) = YMv1x1zsumme37(z,:)+YMv1x1zsummeges37(1,:);
       YMv1x1zsummeges38(1,:) = YMv1x1zsumme38(z,:)+YMv1x1zsummeges38(1,:);
       YMv1x1zsummeges45(1,:) = YMv1x1zsumme45(z,:)+YMv1x1zsummeges45(1,:);
       YMv1x1zsummeges46(1,:) = YMv1x1zsumme46(z,:)+YMv1x1zsummeges46(1,:);
       YMv1x1zsummeges47(1,:) = YMv1x1zsumme47(z,:)+YMv1x1zsummeges47(1,:);
       YMv1x1zsummeges48(1,:) = YMv1x1zsumme48(z,:)+YMv1x1zsummeges48(1,:);
           
       YMv1y1zsummeges15(1,:) = YMv1y1zsumme15(z,:)+YMv1y1zsummeges15(1,:);
       YMv1y1zsummeges16(1,:) = YMv1y1zsumme16(z,:)+YMv1y1zsummeges16(1,:);
       YMv1y1zsummeges17(1,:) = YMv1y1zsumme17(z,:)+YMv1y1zsummeges17(1,:);
       YMv1y1zsummeges18(1,:) = YMv1y1zsumme18(z,:)+YMv1y1zsummeges18(1,:);
       YMv1y1zsummeges25(1,:) = YMv1y1zsumme25(z,:)+YMv1y1zsummeges25(1,:);
       YMv1y1zsummeges26(1,:) = YMv1y1zsumme26(z,:)+YMv1y1zsummeges26(1,:);
       YMv1y1zsummeges27(1,:) = YMv1y1zsumme27(z,:)+YMv1y1zsummeges27(1,:);
       YMv1y1zsummeges28(1,:) = YMv1y1zsumme28(z,:)+YMv1y1zsummeges28(1,:);
       YMv1y1zsummeges35(1,:) = YMv1y1zsumme35(z,:)+YMv1y1zsummeges35(1,:);
       YMv1y1zsummeges36(1,:) = YMv1y1zsumme36(z,:)+YMv1y1zsummeges36(1,:);
       YMv1y1zsummeges37(1,:) = YMv1y1zsumme37(z,:)+YMv1y1zsummeges37(1,:);
       YMv1y1zsummeges38(1,:) = YMv1y1zsumme38(z,:)+YMv1y1zsummeges38(1,:);
       YMv1y1zsummeges45(1,:) = YMv1y1zsumme45(z,:)+YMv1y1zsummeges45(1,:);
       YMv1y1zsummeges46(1,:) = YMv1y1zsumme46(z,:)+YMv1y1zsummeges46(1,:);
       YMv1y1zsummeges47(1,:) = YMv1y1zsumme47(z,:)+YMv1y1zsummeges47(1,:);
       YMv1y1zsummeges48(1,:) = YMv1y1zsumme48(z,:)+YMv1y1zsummeges48(1,:);
       
       YfW1z1xsummeges15(1,:) = YfW1z1xsumme15(z,:)+YfW1z1xsummeges15(1,:);
       YfW1z1xsummeges16(1,:) = YfW1z1xsumme16(z,:)+YfW1z1xsummeges16(1,:);
       YfW1z1xsummeges17(1,:) = YfW1z1xsumme17(z,:)+YfW1z1xsummeges17(1,:);
       YfW1z1xsummeges18(1,:) = YfW1z1xsumme18(z,:)+YfW1z1xsummeges18(1,:);
       YfW1z1xsummeges25(1,:) = YfW1z1xsumme25(z,:)+YfW1z1xsummeges25(1,:);
       YfW1z1xsummeges26(1,:) = YfW1z1xsumme26(z,:)+YfW1z1xsummeges26(1,:);
       YfW1z1xsummeges27(1,:) = YfW1z1xsumme27(z,:)+YfW1z1xsummeges27(1,:);
       YfW1z1xsummeges28(1,:) = YfW1z1xsumme28(z,:)+YfW1z1xsummeges28(1,:);
       YfW1z1xsummeges35(1,:) = YfW1z1xsumme35(z,:)+YfW1z1xsummeges35(1,:);
       YfW1z1xsummeges36(1,:) = YfW1z1xsumme36(z,:)+YfW1z1xsummeges36(1,:);
       YfW1z1xsummeges37(1,:) = YfW1z1xsumme37(z,:)+YfW1z1xsummeges37(1,:);
       YfW1z1xsummeges38(1,:) = YfW1z1xsumme38(z,:)+YfW1z1xsummeges38(1,:);
       YfW1z1xsummeges45(1,:) = YfW1z1xsumme45(z,:)+YfW1z1xsummeges45(1,:);
       YfW1z1xsummeges46(1,:) = YfW1z1xsumme46(z,:)+YfW1z1xsummeges46(1,:);
       YfW1z1xsummeges47(1,:) = YfW1z1xsumme47(z,:)+YfW1z1xsummeges47(1,:);
       YfW1z1xsummeges48(1,:) = YfW1z1xsumme48(z,:)+YfW1z1xsummeges48(1,:);
             
       YM1x1xsummeges15(1,:) = YM1x1xsumme15(z,:)+YM1x1xsummeges15(1,:);
       YM1x1xsummeges16(1,:) = YM1x1xsumme16(z,:)+YM1x1xsummeges16(1,:);
       YM1x1xsummeges17(1,:) = YM1x1xsumme17(z,:)+YM1x1xsummeges17(1,:);
       YM1x1xsummeges18(1,:) = YM1x1xsumme18(z,:)+YM1x1xsummeges18(1,:);
       YM1x1xsummeges25(1,:) = YM1x1xsumme25(z,:)+YM1x1xsummeges25(1,:);
       YM1x1xsummeges26(1,:) = YM1x1xsumme26(z,:)+YM1x1xsummeges26(1,:);
       YM1x1xsummeges27(1,:) = YM1x1xsumme27(z,:)+YM1x1xsummeges27(1,:);
       YM1x1xsummeges28(1,:) = YM1x1xsumme28(z,:)+YM1x1xsummeges28(1,:);
       YM1x1xsummeges35(1,:) = YM1x1xsumme35(z,:)+YM1x1xsummeges35(1,:);
       YM1x1xsummeges36(1,:) = YM1x1xsumme36(z,:)+YM1x1xsummeges36(1,:);
       YM1x1xsummeges37(1,:) = YM1x1xsumme37(z,:)+YM1x1xsummeges37(1,:);
       YM1x1xsummeges38(1,:) = YM1x1xsumme38(z,:)+YM1x1xsummeges38(1,:);
       YM1x1xsummeges45(1,:) = YM1x1xsumme45(z,:)+YM1x1xsummeges45(1,:);
       YM1x1xsummeges46(1,:) = YM1x1xsumme46(z,:)+YM1x1xsummeges46(1,:);
       YM1x1xsummeges47(1,:) = YM1x1xsumme47(z,:)+YM1x1xsummeges47(1,:);
       YM1x1xsummeges48(1,:) = YM1x1xsumme48(z,:)+YM1x1xsummeges48(1,:);
              
       YM1y1xsummeges15(1,:) = YM1y1xsumme15(z,:)+YM1y1xsummeges15(1,:);
       YM1y1xsummeges16(1,:) = YM1y1xsumme16(z,:)+YM1y1xsummeges16(1,:);
       YM1y1xsummeges17(1,:) = YM1y1xsumme17(z,:)+YM1y1xsummeges17(1,:);
       YM1y1xsummeges18(1,:) = YM1y1xsumme18(z,:)+YM1y1xsummeges18(1,:);
       YM1y1xsummeges25(1,:) = YM1y1xsumme25(z,:)+YM1y1xsummeges25(1,:);
       YM1y1xsummeges26(1,:) = YM1y1xsumme26(z,:)+YM1y1xsummeges26(1,:);
       YM1y1xsummeges27(1,:) = YM1y1xsumme27(z,:)+YM1y1xsummeges27(1,:);
       YM1y1xsummeges28(1,:) = YM1y1xsumme28(z,:)+YM1y1xsummeges28(1,:);
       YM1y1xsummeges35(1,:) = YM1y1xsumme35(z,:)+YM1y1xsummeges35(1,:);
       YM1y1xsummeges36(1,:) = YM1y1xsumme36(z,:)+YM1y1xsummeges36(1,:);
       YM1y1xsummeges37(1,:) = YM1y1xsumme37(z,:)+YM1y1xsummeges37(1,:);
       YM1y1xsummeges38(1,:) = YM1y1xsumme38(z,:)+YM1y1xsummeges38(1,:);
       YM1y1xsummeges45(1,:) = YM1y1xsumme45(z,:)+YM1y1xsummeges45(1,:);
       YM1y1xsummeges46(1,:) = YM1y1xsumme46(z,:)+YM1y1xsummeges46(1,:);
       YM1y1xsummeges47(1,:) = YM1y1xsumme47(z,:)+YM1y1xsummeges47(1,:);
       YM1y1xsummeges48(1,:) = YM1y1xsumme48(z,:)+YM1y1xsummeges48(1,:);
       
       YfW1z1ysummeges15(1,:) = YfW1z1ysumme15(z,:)+YfW1z1ysummeges15(1,:);
       YfW1z1ysummeges16(1,:) = YfW1z1ysumme16(z,:)+YfW1z1ysummeges16(1,:);
       YfW1z1ysummeges17(1,:) = YfW1z1ysumme17(z,:)+YfW1z1ysummeges17(1,:);
       YfW1z1ysummeges18(1,:) = YfW1z1ysumme18(z,:)+YfW1z1ysummeges18(1,:);
       YfW1z1ysummeges25(1,:) = YfW1z1ysumme25(z,:)+YfW1z1ysummeges25(1,:);
       YfW1z1ysummeges26(1,:) = YfW1z1ysumme26(z,:)+YfW1z1ysummeges26(1,:);
       YfW1z1ysummeges27(1,:) = YfW1z1ysumme27(z,:)+YfW1z1ysummeges27(1,:);
       YfW1z1ysummeges28(1,:) = YfW1z1ysumme28(z,:)+YfW1z1ysummeges28(1,:);
       YfW1z1ysummeges35(1,:) = YfW1z1ysumme35(z,:)+YfW1z1ysummeges35(1,:);
       YfW1z1ysummeges36(1,:) = YfW1z1ysumme36(z,:)+YfW1z1ysummeges36(1,:);
       YfW1z1ysummeges37(1,:) = YfW1z1ysumme37(z,:)+YfW1z1ysummeges37(1,:);
       YfW1z1ysummeges38(1,:) = YfW1z1ysumme38(z,:)+YfW1z1ysummeges38(1,:);
       YfW1z1ysummeges45(1,:) = YfW1z1ysumme45(z,:)+YfW1z1ysummeges45(1,:);
       YfW1z1ysummeges46(1,:) = YfW1z1ysumme46(z,:)+YfW1z1ysummeges46(1,:);
       YfW1z1ysummeges47(1,:) = YfW1z1ysumme47(z,:)+YfW1z1ysummeges47(1,:);
       YfW1z1ysummeges48(1,:) = YfW1z1ysumme48(z,:)+YfW1z1ysummeges48(1,:);
       
       YM1x1ysummeges15(1,:)  = YM1x1ysumme15(z,:)+YM1x1ysummeges15(1,:);
       YM1x1ysummeges16(1,:)  = YM1x1ysumme16(z,:)+YM1x1ysummeges16(1,:);
       YM1x1ysummeges17(1,:)  = YM1x1ysumme17(z,:)+YM1x1ysummeges17(1,:);
       YM1x1ysummeges18(1,:)  = YM1x1ysumme18(z,:)+YM1x1ysummeges18(1,:);
       YM1x1ysummeges25(1,:)  = YM1x1ysumme25(z,:)+YM1x1ysummeges25(1,:);
       YM1x1ysummeges26(1,:)  = YM1x1ysumme26(z,:)+YM1x1ysummeges26(1,:);
       YM1x1ysummeges27(1,:)  = YM1x1ysumme27(z,:)+YM1x1ysummeges27(1,:);
       YM1x1ysummeges28(1,:)  = YM1x1ysumme28(z,:)+YM1x1ysummeges28(1,:);
       YM1x1ysummeges35(1,:)  = YM1x1ysumme35(z,:)+YM1x1ysummeges35(1,:);
       YM1x1ysummeges36(1,:)  = YM1x1ysumme36(z,:)+YM1x1ysummeges36(1,:);
       YM1x1ysummeges37(1,:)  = YM1x1ysumme37(z,:)+YM1x1ysummeges37(1,:);
       YM1x1ysummeges38(1,:)  = YM1x1ysumme38(z,:)+YM1x1ysummeges38(1,:);
       YM1x1ysummeges45(1,:)  = YM1x1ysumme45(z,:)+YM1x1ysummeges45(1,:);
       YM1x1ysummeges46(1,:)  = YM1x1ysumme46(z,:)+YM1x1ysummeges46(1,:);
       YM1x1ysummeges47(1,:)  = YM1x1ysumme47(z,:)+YM1x1ysummeges47(1,:);
       YM1x1ysummeges48(1,:)  = YM1x1ysumme48(z,:)+YM1x1ysummeges48(1,:);
           
       YM1y1ysummeges15(1,:)  = YM1y1ysumme15(z,:)+YM1y1ysummeges15(1,:);
       YM1y1ysummeges16(1,:)  = YM1y1ysumme16(z,:)+YM1y1ysummeges16(1,:);
       YM1y1ysummeges17(1,:)  = YM1y1ysumme17(z,:)+YM1y1ysummeges17(1,:);
       YM1y1ysummeges18(1,:)  = YM1y1ysumme18(z,:)+YM1y1ysummeges18(1,:);
       YM1y1ysummeges25(1,:)  = YM1y1ysumme25(z,:)+YM1y1ysummeges25(1,:);
       YM1y1ysummeges26(1,:)  = YM1y1ysumme26(z,:)+YM1y1ysummeges26(1,:);
       YM1y1ysummeges27(1,:)  = YM1y1ysumme27(z,:)+YM1y1ysummeges27(1,:);
       YM1y1ysummeges28(1,:)  = YM1y1ysumme28(z,:)+YM1y1ysummeges28(1,:);
       YM1y1ysummeges35(1,:)  = YM1y1ysumme35(z,:)+YM1y1ysummeges35(1,:);
       YM1y1ysummeges36(1,:)  = YM1y1ysumme36(z,:)+YM1y1ysummeges36(1,:);
       YM1y1ysummeges37(1,:)  = YM1y1ysumme37(z,:)+YM1y1ysummeges37(1,:);
       YM1y1ysummeges38(1,:)  = YM1y1ysumme38(z,:)+YM1y1ysummeges38(1,:);
       YM1y1ysummeges45(1,:)  = YM1y1ysumme45(z,:)+YM1y1ysummeges45(1,:);
       YM1y1ysummeges46(1,:)  = YM1y1ysumme46(z,:)+YM1y1ysummeges46(1,:);
       YM1y1ysummeges47(1,:)  = YM1y1ysumme47(z,:)+YM1y1ysummeges47(1,:);
       YM1y1ysummeges48(1,:)  = YM1y1ysumme48(z,:)+YM1y1ysummeges48(1,:);

    if (fmn < fmax)
        m = m + 1;
    elseif (m > 0)
        n = n + 1;
        m = 0;
        z = z - 1;
    else
        break;
    end
end

YWFz15 = abs(Yf1z1zsummeges15(1,:).*j.*w);  % T = theta!
YWFz16 = abs(Yf1z1zsummeges16(1,:).*j.*w);
YWFz17 = abs(Yf1z1zsummeges17(1,:).*j.*w);
YWFz18 = abs(Yf1z1zsummeges18(1,:).*j.*w);
YWFz25 = abs(Yf1z1zsummeges25(1,:).*j.*w);
YWFz26 = abs(Yf1z1zsummeges26(1,:).*j.*w);
YWFz27 = abs(Yf1z1zsummeges27(1,:).*j.*w);
YWFz28 = abs(Yf1z1zsummeges28(1,:).*j.*w);
YWFz35 = abs(Yf1z1zsummeges35(1,:).*j.*w);
YWFz36 = abs(Yf1z1zsummeges36(1,:).*j.*w);
YWFz37 = abs(Yf1z1zsummeges37(1,:).*j.*w);
YWFz38 = abs(Yf1z1zsummeges38(1,:).*j.*w);
YWFz45 = abs(Yf1z1zsummeges45(1,:).*j.*w);
YWFz46 = abs(Yf1z1zsummeges46(1,:).*j.*w);
YWFz47 = abs(Yf1z1zsummeges47(1,:).*j.*w);
YWFz48 = abs(Yf1z1zsummeges48(1,:).*j.*w);

YWMx15 = abs(YMv1x1zsummeges15(1,:).*j.*w); 
YWMx16 = abs(YMv1x1zsummeges16(1,:).*j.*w); 
YWMx17 = abs(YMv1x1zsummeges17(1,:).*j.*w); 
YWMx18 = abs(YMv1x1zsummeges18(1,:).*j.*w); 
YWMx25 = abs(YMv1x1zsummeges25(1,:).*j.*w); 
YWMx26 = abs(YMv1x1zsummeges26(1,:).*j.*w); 
YWMx27 = abs(YMv1x1zsummeges27(1,:).*j.*w); 
YWMx28 = abs(YMv1x1zsummeges28(1,:).*j.*w); 
YWMx35 = abs(YMv1x1zsummeges35(1,:).*j.*w); 
YWMx36 = abs(YMv1x1zsummeges36(1,:).*j.*w); 
YWMx37 = abs(YMv1x1zsummeges37(1,:).*j.*w); 
YWMx38 = abs(YMv1x1zsummeges38(1,:).*j.*w); 
YWMx45 = abs(YMv1x1zsummeges45(1,:).*j.*w); 
YWMx46 = abs(YMv1x1zsummeges46(1,:).*j.*w); 
YWMx47 = abs(YMv1x1zsummeges47(1,:).*j.*w); 
YWMx48 = abs(YMv1x1zsummeges48(1,:).*j.*w); 

YWMy15 = abs(YMv1y1zsummeges15(1,:).*j.*w);
YWMy16 = abs(YMv1y1zsummeges16(1,:).*j.*w);
YWMy17 = abs(YMv1y1zsummeges17(1,:).*j.*w);
YWMy18 = abs(YMv1y1zsummeges18(1,:).*j.*w);
YWMy25 = abs(YMv1y1zsummeges25(1,:).*j.*w);
YWMy26 = abs(YMv1y1zsummeges26(1,:).*j.*w);
YWMy27 = abs(YMv1y1zsummeges27(1,:).*j.*w);
YWMy28 = abs(YMv1y1zsummeges28(1,:).*j.*w);
YWMy35 = abs(YMv1y1zsummeges35(1,:).*j.*w);
YWMy36 = abs(YMv1y1zsummeges36(1,:).*j.*w);
YWMy37 = abs(YMv1y1zsummeges37(1,:).*j.*w);
YWMy38 = abs(YMv1y1zsummeges38(1,:).*j.*w);
YWMy45 = abs(YMv1y1zsummeges45(1,:).*j.*w);
YWMy46 = abs(YMv1y1zsummeges46(1,:).*j.*w);
YWMy47 = abs(YMv1y1zsummeges47(1,:).*j.*w);
YWMy48 = abs(YMv1y1zsummeges48(1,:).*j.*w);

YTxFz15 = abs(YfW1z1xsummeges15(1,:).*j.*w); 
YTxFz16 = abs(YfW1z1xsummeges16(1,:).*j.*w);
YTxFz17 = abs(YfW1z1xsummeges17(1,:).*j.*w); 
YTxFz18 = abs(YfW1z1xsummeges18(1,:).*j.*w);
YTxFz25 = abs(YfW1z1xsummeges25(1,:).*j.*w); 
YTxFz26 = abs(YfW1z1xsummeges26(1,:).*j.*w);
YTxFz27 = abs(YfW1z1xsummeges27(1,:).*j.*w); 
YTxFz28 = abs(YfW1z1xsummeges28(1,:).*j.*w);
YTxFz35 = abs(YfW1z1xsummeges35(1,:).*j.*w); 
YTxFz36 = abs(YfW1z1xsummeges36(1,:).*j.*w);
YTxFz37 = abs(YfW1z1xsummeges37(1,:).*j.*w); 
YTxFz38 = abs(YfW1z1xsummeges38(1,:).*j.*w);
YTxFz45 = abs(YfW1z1xsummeges45(1,:).*j.*w); 
YTxFz46 = abs(YfW1z1xsummeges46(1,:).*j.*w);
YTxFz47 = abs(YfW1z1xsummeges47(1,:).*j.*w); 
YTxFz48 = abs(YfW1z1xsummeges48(1,:).*j.*w);

YTxMx15 = abs(YM1x1xsummeges15(1,:).*j.*w);
YTxMx16 = abs(YM1x1xsummeges16(1,:).*j.*w);
YTxMx17 = abs(YM1x1xsummeges17(1,:).*j.*w);
YTxMx18 = abs(YM1x1xsummeges18(1,:).*j.*w);
YTxMx25 = abs(YM1x1xsummeges25(1,:).*j.*w);
YTxMx26 = abs(YM1x1xsummeges26(1,:).*j.*w);
YTxMx27 = abs(YM1x1xsummeges27(1,:).*j.*w);
YTxMx28 = abs(YM1x1xsummeges28(1,:).*j.*w);
YTxMx35 = abs(YM1x1xsummeges35(1,:).*j.*w);
YTxMx36 = abs(YM1x1xsummeges36(1,:).*j.*w);
YTxMx37 = abs(YM1x1xsummeges37(1,:).*j.*w);
YTxMx38 = abs(YM1x1xsummeges38(1,:).*j.*w);
YTxMx45 = abs(YM1x1xsummeges45(1,:).*j.*w);
YTxMx46 = abs(YM1x1xsummeges46(1,:).*j.*w);
YTxMx47 = abs(YM1x1xsummeges47(1,:).*j.*w);
YTxMx48 = abs(YM1x1xsummeges48(1,:).*j.*w);

YTxMy15 = abs(YM1y1xsummeges15(1,:) .*j.*w); 
YTxMy16 = abs(YM1y1xsummeges16(1,:) .*j.*w);
YTxMy17 = abs(YM1y1xsummeges17(1,:) .*j.*w);
YTxMy18 = abs(YM1y1xsummeges18(1,:) .*j.*w);
YTxMy25 = abs(YM1y1xsummeges25(1,:) .*j.*w); 
YTxMy26 = abs(YM1y1xsummeges26(1,:) .*j.*w);
YTxMy27 = abs(YM1y1xsummeges27(1,:) .*j.*w);
YTxMy28 = abs(YM1y1xsummeges28(1,:) .*j.*w);
YTxMy35 = abs(YM1y1xsummeges35(1,:) .*j.*w); 
YTxMy36 = abs(YM1y1xsummeges36(1,:) .*j.*w);
YTxMy37 = abs(YM1y1xsummeges37(1,:) .*j.*w);
YTxMy38 = abs(YM1y1xsummeges38(1,:) .*j.*w);
YTxMy45 = abs(YM1y1xsummeges45(1,:) .*j.*w); 
YTxMy46 = abs(YM1y1xsummeges46(1,:) .*j.*w);
YTxMy47 = abs(YM1y1xsummeges47(1,:) .*j.*w);
YTxMy48 = abs(YM1y1xsummeges48(1,:) .*j.*w);

YTyFz15 = abs(YfW1z1ysummeges15(1,:).*j.*w);
YTyFz16 = abs(YfW1z1ysummeges16(1,:).*j.*w);
YTyFz17 = abs(YfW1z1ysummeges17(1,:).*j.*w);
YTyFz18 = abs(YfW1z1ysummeges18(1,:).*j.*w);
YTyFz25 = abs(YfW1z1ysummeges15(1,:).*j.*w);
YTyFz26 = abs(YfW1z1ysummeges16(1,:).*j.*w);
YTyFz27 = abs(YfW1z1ysummeges17(1,:).*j.*w);
YTyFz28 = abs(YfW1z1ysummeges18(1,:).*j.*w);
YTyFz35 = abs(YfW1z1ysummeges35(1,:).*j.*w);
YTyFz36 = abs(YfW1z1ysummeges36(1,:).*j.*w);
YTyFz37 = abs(YfW1z1ysummeges37(1,:).*j.*w);
YTyFz38 = abs(YfW1z1ysummeges38(1,:).*j.*w);
YTyFz45 = abs(YfW1z1ysummeges45(1,:).*j.*w);
YTyFz46 = abs(YfW1z1ysummeges46(1,:).*j.*w);
YTyFz47 = abs(YfW1z1ysummeges47(1,:).*j.*w);
YTyFz48 = abs(YfW1z1ysummeges48(1,:).*j.*w);

YTyMx15 = abs(YM1x1ysummeges15(1,:) .*j.*w);
YTyMx16 = abs(YM1x1ysummeges16(1,:) .*j.*w);
YTyMx17 = abs(YM1x1ysummeges17(1,:) .*j.*w);
YTyMx18 = abs(YM1x1ysummeges18(1,:) .*j.*w);
YTyMx25 = abs(YM1x1ysummeges25(1,:) .*j.*w);
YTyMx26 = abs(YM1x1ysummeges26(1,:) .*j.*w);
YTyMx27 = abs(YM1x1ysummeges27(1,:) .*j.*w);
YTyMx28 = abs(YM1x1ysummeges28(1,:) .*j.*w);
YTyMx35 = abs(YM1x1ysummeges35(1,:) .*j.*w);
YTyMx36 = abs(YM1x1ysummeges36(1,:) .*j.*w);
YTyMx37 = abs(YM1x1ysummeges37(1,:) .*j.*w);
YTyMx38 = abs(YM1x1ysummeges38(1,:) .*j.*w);
YTyMx45 = abs(YM1x1ysummeges45(1,:) .*j.*w);
YTyMx46 = abs(YM1x1ysummeges46(1,:) .*j.*w);
YTyMx47 = abs(YM1x1ysummeges47(1,:) .*j.*w);
YTyMx48 = abs(YM1x1ysummeges48(1,:) .*j.*w);

YTyMy15 = abs(YM1y1ysummeges15(1,:).*j.*w);
YTyMy16 = abs(YM1y1ysummeges16(1,:).*j.*w);
YTyMy17 = abs(YM1y1ysummeges17(1,:).*j.*w);
YTyMy18 = abs(YM1y1ysummeges18(1,:).*j.*w);
YTyMy25 = abs(YM1y1ysummeges25(1,:).*j.*w);
YTyMy26 = abs(YM1y1ysummeges26(1,:).*j.*w);
YTyMy27 = abs(YM1y1ysummeges27(1,:).*j.*w);
YTyMy28 = abs(YM1y1ysummeges28(1,:).*j.*w);
YTyMy35 = abs(YM1y1ysummeges35(1,:).*j.*w);
YTyMy36 = abs(YM1y1ysummeges36(1,:).*j.*w);
YTyMy37 = abs(YM1y1ysummeges37(1,:).*j.*w);
YTyMy38 = abs(YM1y1ysummeges38(1,:).*j.*w);
YTyMy45 = abs(YM1y1ysummeges45(1,:).*j.*w);
YTyMy46 = abs(YM1y1ysummeges46(1,:).*j.*w);
YTyMy47 = abs(YM1y1ysummeges47(1,:).*j.*w);
YTyMy48 = abs(YM1y1ysummeges48(1,:).*j.*w);

Mob_Matrix = [YWFz15 YWFz16 YWFz17 YWFz18; 
              YWFz25 YWFz26 YWFz27 YWFz28;
              YWFz35 YWFz36 YWFz37 YWFz38;
              YWFz45 YWFz46 YWFz47 YWFz48
              YWMx15 YWMx16 YWMx17 YWMx18
              YWMx25 YWMx26 YWMx27 YWMx28
              YWMx35 YWMx36 YWMx37 YWMx38
              YWMx45 YWMx46 YWMx47 YWMx48
              YWMy15 YWMy16 YWMy17 YWMy18
              YWMy25 YWMy16 YWMy17 YWMy18
              YWMy35 YWMy36 YWMy37 YWMy38
              YWMy45 YWMy46 YWMy47 YWMy48];
          
          
v_matrix = [1;1;1;1;1;1;1;1;1;1;1;1];

F_matrix = Mob_Matrix*v_matrix;