for ii=1:1:100
 
a(ii)=ii*0.00001;
%,00.0014

    
element{1} = struct('hs',P,'hp',G,'b',0.0015,'L',0.011 ,'rhos',11340,'E',complex(0.16E11,0),'poisson',0.26,'k_schub',sqrt(0.87),'s33E',20.5E-12,'s11E',15.8E-12,'epsT',7000*8.85E-12,'d31',-385E-12,'Qm',45,'rhop',8000,'tandelta',0.025,'s55E',57.8E-12,'N',0,'l_klemm',0);


for kk=1:100;

% Frequenz
fall =10*kk ;

% Lasten   
load{1} =[1    1 ;
         0    0 ];
load{2} =[1    2 ; 
         0    0 ];   
load{3} = [1   5 ; 
          130  0 ];
load{4} = [2   3 ; 
          0   0 ]; 
% load{5} = [2   4 ; 
%           0   0 ]; 

%Diskrete Elemente
diskretEle{1}  = struct('type','maxwell+','iRand',2,'c',C3(kk),'d',D(kk));



% Zusammenfassung zu Aufgaben Struct
problem.load        = load;
problem.kopplungen  = kopplungen;
problem.element     = element;
problem.linEqn      = linEqn;
problem.diskretEle  = diskretEle;

problem = complete_Input(problem);


% L?sung
for i=1:length(fall)
    problem.w  = fall(i)*2*pi;
    resAll(i,:,:)  = make_Calculation(problem);                                                             
end
Freq(kk)=fall;
geschw(kk)=resAll(:,2,1);

end
geschwindigkeit(:,ii,hh)=geschw;
weg(:,ii,hh)=geschwindigkeit(:,ii,hh)./fall'/2/pi/1i;
end





s_a=size(a);
s_b=size(b);

        figure 
        hold on
        for gg=1:s_b(2)

                surf(a,Freq,abs(geschwindigkeit(:,:,gg)))
               title('Auslenkung ','FontSize',24);xlabel('L?nge [m]');ylabel('Frequenz[Hz]');zlabel('Auslenkung[m]');
    
        end