%%%%%%%%%%%%%%%%%%
%Fuktion%%%%%%%%%%
   
    function f = myfun(x)
  
    [Nr, x01, y]   = textread('KontaktNodesStart.txt','%f %f %f', -1);   

    AnzVec     = size(x01);
    Anz        = AnzVec(1);
    
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  
       
    fidAnsysNmodif = fopen('AnsysNmodif.txt', 'wt');
    
    for i=1:Anz

        s1 = 'nmodif,';
        s2 = ',';
        
        fprintf(fidAnsysNmodif, '%s %f %s %f %s %f\n ',s1,Nr(i),s2,x(i),s2,y(i));

    end
    
    fclose(fidAnsysNmodif)

    fidX = fopen('X.txt', 'wt');
    
    for i=1:Anz

        
        fprintf(fidX, '%f\n ',x(i));

    end
    
    fclose(fidX)
    
    x
    
    pause
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    %Ansys Aufruf

     % !module add ansys

     % !ansys110 -p aa_t_a -b -j HomModell < auto.txt > HomModell.out 
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    
    [Nr, Stress] = textread('Stress.txt','%f %f', -1); 

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

    StressMean = mean(Stress); 
    %Xmean = mean(x);
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   
    for i = 1:Anz
        
        Resid(i) = abs(StressMean - Stress(i));
        %Resid(i) = abs(Xmean - x(i));
        
    end
    
    f = sum(Resid);
    %f
    %f = mean(x);
     
    %x
    %f
    %pause
    
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Funktion Ende
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%