clear 

[Nr, x0, y]   = textread('KontaktNodesStart.txt','%f %f %f', -1);   

AnzVec     = size(x0);
Anz        = AnzVec(1);

for i=1:31
    
    lb(i)  = 18;
    ub(i)  = 20;
    
end


b = [18; 20];

option=optimset('TolX',1e-10,'Display','iter','DiffMinChange',0.005,'TypicalX',19,'TolFun',1);
[x, fval,exitflag,output] = fmincon(@myfun,x0,[],[],[],[],lb,ub,[])
%x = fminunc(@myfun,x0,option);

fidAnsysXY = fopen('AnsysXYopti.txt', 'wt');
    
    for i=1:Anz

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

    end

fclose(fidAnsysXY)

fidResid = fopen('Resid.txt', 'wt');
fprintf(fidResid, '%f\n ',fval);
fclose(fidResid)