clear all;
clear globals;
close all;
clc;
tic;

%Nacteni identickych bodu
dir_ib  = '_data\_input\ib\';
names_ib = dir([dir_ib,'*_ib.txt']);
ib_files = size(names_ib,1);

% Matice A a vektor l
radek = 0;
for i=1:ib_files
    points_ib = (load([dir_ib,names_ib(i).name]))';
    ib_count = size(points_ib,2);
    for j=1:ib_count
       A(2*radek+j,1+12*(i-1)) = points_ib(1,j)*points_ib(1,j);
       A(2*radek+j,2+12*(i-1)) = points_ib(2,j)*points_ib(2,j);
       A(2*radek+j,3+12*(i-1)) = points_ib(1,j)*points_ib(2,j);
       A(2*radek+j,4+12*(i-1)) = points_ib(1,j);
       A(2*radek+j,5+12*(i-1)) = points_ib(2,j);
       A(2*radek+j,6+12*(i-1)) = 1;
       A(2*radek+j+ib_count,7+12*(i-1)) = points_ib(1,j)*points_ib(1,j);
       A(2*radek+j+ib_count,8+12*(i-1)) = points_ib(2,j)*points_ib(2,j);
       A(2*radek+j+ib_count,9+12*(i-1)) = points_ib(1,j)*points_ib(2,j);
       A(2*radek+j+ib_count,10+12*(i-1)) = points_ib(1,j);
       A(2*radek+j+ib_count,11+12*(i-1)) = points_ib(2,j);
       A(2*radek+j+ib_count,12+12*(i-1)) = 1;    
       vektor_l(2*radek+j) = points_ib(3,j);
       vektor_l(2*radek+j+ib_count) = points_ib(4,j)      
    end  
    radek = radek + ib_count;
    misto = findstr(names_ib(i).name, '_');
    jmeno = str2num(names_ib(i).name(1:misto(1)-1));
    jmena_ib(i) = jmeno;
end
vektor_l = vektor_l';


% Nacteni souradnic rohu
dir_rohy  = '_data\_input\rohy\';
names_rohy = dir([dir_rohy,'*_rohy.txt']);
rohy_files = size(names_rohy,1);

for i=1:rohy_files
    points_rohy = (load([dir_rohy,names_rohy(i).name]))';
    rohy_count = size(points_rohy,2);
    for j=1:rohy_count
        Bpomocna(i,1) = points_rohy(1,1);
        Bpomocna(i,2) = points_rohy(2,1);
        Bpomocna(i,3) = points_rohy(1,2);
        Bpomocna(i,4) = points_rohy(2,2);   
        Bpomocna(i,5) = points_rohy(1,3);   
        Bpomocna(i,6) = points_rohy(2,3);   
        Bpomocna(i,7) = points_rohy(1,4);   
        Bpomocna(i,8) = points_rohy(2,4);          
    end
end    



%Podminky pro rohy(1-horni dva, 2-prave dva, 3-dolni dva, 4-leve dva,
%5-bez LD, 6-bez LH, 7-bez PH, 8-bez PD, 9-vse, 10-LH+PD, 11-PH+LD)

dir_klad  = '_data\_input\';
names_klad = dir([dir_klad,'rozlozeni.txt']);
klad = (load([dir_klad,names_klad(1).name]))';
klad = klad';

podminka_typ=zeros(size(klad,1)+1,size(klad,2)+1);
for i=1:(size(klad,1)+1)
    for j=1:(size(klad,2)+1)
        if (i==1) && (j==1)
            podminka_typ(i,j)=0;
        elseif (i==1) && (j==size(klad,2)+1)
            podminka_typ(i,j)=0;
        elseif (i==size(klad,1)+1) && (j==size(klad,2)+1)
            podminka_typ(i,j)=0;
        elseif (i==size(klad,1)+1) && (j==1)
            podminka_typ(i,j)=0;
                
        elseif (i==1) && (j>1) && (j<=size(klad,2))
            if (klad(1,j-1)==0) && (klad(1,j)==0)
                podminka_typ(i,j)=0;
            elseif (klad(1,j-1)==0) && (klad(1,j)~=0)    
                podminka_typ(i,j)=0;
            elseif (klad(1,j-1)~=0) && (klad(1,j)==0)        
                podminka_typ(i,j)=0;
            elseif (klad(1,j-1)~=0) && (klad(1,j)~=0)        
                podminka_typ(i,j)=3;
            end
        elseif (i>1) && (i<=size(klad,1)) && (j==1)
            if (klad(i-1,1)==0) && (klad(i,1)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,1)==0) && (klad(i,1)~=0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,1)~=0) && (klad(i,1)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,1)~=0) && (klad(i,1)~=0)
                podminka_typ(i,j)=2;
            end
        elseif (i==size(klad,1)+1) && (j>1) && (j<=size(klad,2))
            if (klad(i-1,j-1)==0) && (klad(i-1,j)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)~=0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)~=0)
                podminka_typ(i,j)=1;
            end
        elseif (i>1) && (i<=size(klad,1)) && (j==size(klad,2)+1)
            if (klad(i-1,j-1)==0) && (klad(i,j-1)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)==0) && (klad(i,j-1)~=0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i,j-1)==0)
                podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i,j-1)~=0)
                podminka_typ(i,j)=4;
            end
            
        elseif (i>1) && (i<=size(klad,1)) && (j>1) && (j<=size(klad,2))
            if (klad(i-1,j-1)==0) && (klad(i-1,j)==0) && (klad(i,j)==0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)==0) && (klad(i,j)==0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)~=0) && (klad(i,j)==0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)==0) && (klad(i,j)~=0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)==0) && (klad(i,j)==0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=0;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)~=0) && (klad(i,j)==0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=1;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)~=0) && (klad(i,j)~=0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=2;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)==0) && (klad(i,j)~=0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=3;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)==0) && (klad(i,j)==0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=4;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)==0) && (klad(i,j)~=0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=10;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)~=0) && (klad(i,j)==0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=11;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)~=0) && (klad(i,j)~=0) && (klad(i,j-1)==0)
                  podminka_typ(i,j)=5;
            elseif (klad(i-1,j-1)==0) && (klad(i-1,j)~=0) && (klad(i,j)~=0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=6;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)==0) && (klad(i,j)~=0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=7;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)~=0) && (klad(i,j)==0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=8;
            elseif (klad(i-1,j-1)~=0) && (klad(i-1,j)~=0) && (klad(i,j)~=0) && (klad(i,j-1)~=0)
                  podminka_typ(i,j)=9;
            end
        end
    end
end


%Naplneni Matice B dle podminek pro rohy
podminka_pocet=1;
for i=1:(size(klad,1)+1)
    for j=1:(size(klad,2)+1)
        if podminka_typ(i,j)==1
            lh=klad(i-1,j-1);
            ph=klad(i-1,j);
            lh_index=find(jmena_ib==lh);
            ph_index=find(jmena_ib==ph);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (lh_index-1)*12+1)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= -1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= -1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
            
        elseif podminka_typ(i,j)==2
            ph=klad(i-1,j);
            pd=klad(i,j);
            ph_index=find(jmena_ib==ph);
            pd_index=find(jmena_ib==pd);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (pd_index-1)*12+1)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= -1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= -1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;        
                
        elseif podminka_typ(i,j)==3
            pd=klad(i,j);
            ld=klad(i,j-1);
            pd_index=find(jmena_ib==pd);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (pd_index-1)*12+1)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= 1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= 1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;  
            
        elseif podminka_typ(i,j)==4
            lh=klad(i-1,j-1);
            ld=klad(i,j-1);
            lh_index=find(jmena_ib==lh);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (lh_index-1)*12+1)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= 1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= 1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;  
            
        elseif podminka_typ(i,j)==10
            lh=klad(i-1,j-1);
            pd=klad(i,j);
            lh_index=find(jmena_ib==lh);
            pd_index=find(jmena_ib==pd);
            B(podminka_pocet, (lh_index-1)*12+1)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= 1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet, (pd_index-1)*12+1)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= -1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= 1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= -1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2; 
            
        elseif podminka_typ(i,j)==11
            ph=klad(i-1,j);
            ld=klad(i,j);
            ph_index=find(jmena_ib==ph);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;     
            
        elseif podminka_typ(i,j)==5
            lh=klad(i-1,j-1);
            ph=klad(i-1,j);
            lh_index=find(jmena_ib==lh);
            ph_index=find(jmena_ib==ph);
             B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (lh_index-1)*12+1)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= -1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= -1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;            
            podminka_pocet=podminka_pocet+2;
            
            ph=klad(i-1,j);
            pd=klad(i,j);
            ph_index=find(jmena_ib==ph);
            pd_index=find(jmena_ib==pd);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (pd_index-1)*12+1)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= -1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= -1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;     
            
        elseif podminka_typ(i,j)==6
            ph=klad(i-1,j);
            pd=klad(i,j);
            ph_index=find(jmena_ib==ph);
            pd_index=find(jmena_ib==pd);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (pd_index-1)*12+1)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= -1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= -1*Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= -1*Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= -1*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= -1*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= -1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;   
            
            pd=klad(i,j);
            ld=klad(i,j-1);
            pd_index=find(jmena_ib==pd);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (pd_index-1)*12+1)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= 1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= 1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;  
            
        elseif podminka_typ(i,j)==7
            pd=klad(i,j);
            ld=klad(i,j-1);
            pd_index=find(jmena_ib==pd);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (pd_index-1)*12+1)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= 1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= 1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2; 
            
            lh=klad(i-1,j-1);
            ld=klad(i,j-1);
            lh_index=find(jmena_ib==lh);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (lh_index-1)*12+1)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= 1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= 1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
            
        elseif podminka_typ(i,j)==8
            lh=klad(i-1,j-1);
            ph=klad(i-1,j);
            lh_index=find(jmena_ib==lh);
            ph_index=find(jmena_ib==ph);
             B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (lh_index-1)*12+1)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= -1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= -1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;              
            podminka_pocet=podminka_pocet+2;
            
            lh=klad(i-1,j-1);
            ld=klad(i,j-1);
            lh_index=find(jmena_ib==lh);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (lh_index-1)*12+1)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= 1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= 1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
            
        elseif podminka_typ(i,j)==9
            lh=klad(i-1,j-1);
            ph=klad(i-1,j);
            lh_index=find(jmena_ib==lh);
            ph_index=find(jmena_ib==ph);
            B(podminka_pocet, (ph_index-1)*12+1)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+2)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+3)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+4)= Bpomocna(ph_index,7);
            B(podminka_pocet, (ph_index-1)*12+5)= Bpomocna(ph_index,8);
            B(podminka_pocet, (ph_index-1)*12+6)= 1;
            B(podminka_pocet, (ph_index-1)*12+7)= 0;
            B(podminka_pocet, (ph_index-1)*12+8)= 0;
            B(podminka_pocet, (ph_index-1)*12+9)= 0;
            B(podminka_pocet, (ph_index-1)*12+10)= 0;
            B(podminka_pocet, (ph_index-1)*12+11)= 0;
            B(podminka_pocet, (ph_index-1)*12+12)= 0;
            B(podminka_pocet, (lh_index-1)*12+1)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= -1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+7)= Bpomocna(ph_index,7)*Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+8)= Bpomocna(ph_index,8)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+9)= Bpomocna(ph_index,7)*Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+10)= Bpomocna(ph_index,7);
            B(podminka_pocet+1, (ph_index-1)*12+11)= Bpomocna(ph_index,8);
            B(podminka_pocet+1, (ph_index-1)*12+12)= 1;
            B(podminka_pocet+1, (ph_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ph_index-1)*12+6)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= -1*Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= -1*Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= -1*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= -1*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= -1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;              
            podminka_pocet=podminka_pocet+2;
            
            pd=klad(i,j);
            ld=klad(i,j-1);
            pd_index=find(jmena_ib==pd);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (pd_index-1)*12+1)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+2)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+3)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+4)= Bpomocna(pd_index,1);
            B(podminka_pocet, (pd_index-1)*12+5)= Bpomocna(pd_index,2);
            B(podminka_pocet, (pd_index-1)*12+6)= 1;
            B(podminka_pocet, (pd_index-1)*12+7)= 0;
            B(podminka_pocet, (pd_index-1)*12+8)= 0;
            B(podminka_pocet, (pd_index-1)*12+9)= 0;
            B(podminka_pocet, (pd_index-1)*12+10)= 0;
            B(podminka_pocet, (pd_index-1)*12+11)= 0;
            B(podminka_pocet, (pd_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+7)= Bpomocna(pd_index,1)*Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+8)= Bpomocna(pd_index,2)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+9)= Bpomocna(pd_index,1)*Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+10)= Bpomocna(pd_index,1);
            B(podminka_pocet+1, (pd_index-1)*12+11)= Bpomocna(pd_index,2);
            B(podminka_pocet+1, (pd_index-1)*12+12)= 1;
            B(podminka_pocet+1, (pd_index-1)*12+1)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+2)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+3)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+4)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+5)= 0;
            B(podminka_pocet+1, (pd_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2; 
            
            lh=klad(i-1,j-1);
            ld=klad(i,j-1);
            lh_index=find(jmena_ib==lh);
            ld_index=find(jmena_ib==ld);
            B(podminka_pocet, (lh_index-1)*12+1)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+2)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+3)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+4)= Bpomocna(lh_index,5);
            B(podminka_pocet, (lh_index-1)*12+5)= Bpomocna(lh_index,6);
            B(podminka_pocet, (lh_index-1)*12+6)= 1;
            B(podminka_pocet, (lh_index-1)*12+7)= 0;
            B(podminka_pocet, (lh_index-1)*12+8)= 0;
            B(podminka_pocet, (lh_index-1)*12+9)= 0;
            B(podminka_pocet, (lh_index-1)*12+10)= 0;
            B(podminka_pocet, (lh_index-1)*12+11)= 0;
            B(podminka_pocet, (lh_index-1)*12+12)= 0;
            B(podminka_pocet, (ld_index-1)*12+1)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+2)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+3)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+4)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet, (ld_index-1)*12+5)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet, (ld_index-1)*12+6)= -1;
            B(podminka_pocet, (ld_index-1)*12+7)= 0;
            B(podminka_pocet, (ld_index-1)*12+8)= 0;
            B(podminka_pocet, (ld_index-1)*12+9)= 0;
            B(podminka_pocet, (ld_index-1)*12+10)= 0;
            B(podminka_pocet, (ld_index-1)*12+11)= 0;
            B(podminka_pocet, (ld_index-1)*12+12)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+7)= Bpomocna(lh_index,5)*Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+8)= Bpomocna(lh_index,6)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+9)= Bpomocna(lh_index,5)*Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+10)= Bpomocna(lh_index,5);
            B(podminka_pocet+1, (lh_index-1)*12+11)= Bpomocna(lh_index,6);
            B(podminka_pocet+1, (lh_index-1)*12+12)= 1;
            B(podminka_pocet+1, (lh_index-1)*12+1)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+2)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+3)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+4)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+5)= 0;
            B(podminka_pocet+1, (lh_index-1)*12+6)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+7)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+8)= -1*Bpomocna(ld_index,4)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+9)= -1*Bpomocna(ld_index,3)*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+10)= -1*Bpomocna(ld_index,3);
            B(podminka_pocet+1, (ld_index-1)*12+11)= -1*Bpomocna(ld_index,4);
            B(podminka_pocet+1, (ld_index-1)*12+12)= -1;
            B(podminka_pocet+1, (ld_index-1)*12+1)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+2)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+3)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+4)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+5)= 0;
            B(podminka_pocet+1, (ld_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
            
        end
    end
end

%Podminky pro stredy hran
for i=1:(size(klad,1))
    for j=1:(size(klad,2)-1)
        if klad(i,j)~=0 && klad(i,j+1)~=0                 
            l=klad(i,j);
            p=klad(i,j+1);
            l_index=find(jmena_ib==l);
            p_index=find(jmena_ib==p);
            B(podminka_pocet, (p_index-1)*12+1)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2)*((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2);
            B(podminka_pocet, (p_index-1)*12+2)= ((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2)*((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet, (p_index-1)*12+3)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2)*((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet, (p_index-1)*12+4)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2);
            B(podminka_pocet, (p_index-1)*12+5)= ((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet, (p_index-1)*12+6)= 1;
            B(podminka_pocet, (p_index-1)*12+7)= 0;
            B(podminka_pocet, (p_index-1)*12+8)= 0;
            B(podminka_pocet, (p_index-1)*12+9)= 0;
            B(podminka_pocet, (p_index-1)*12+10)= 0;
            B(podminka_pocet, (p_index-1)*12+11)= 0;
            B(podminka_pocet, (p_index-1)*12+12)= 0;
            B(podminka_pocet, (l_index-1)*12+1)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2)*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2);
            B(podminka_pocet, (l_index-1)*12+2)= -1*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2)*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet, (l_index-1)*12+3)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2)*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet, (l_index-1)*12+4)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2);
            B(podminka_pocet, (l_index-1)*12+5)= -1*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet, (l_index-1)*12+6)= -1;
            B(podminka_pocet, (l_index-1)*12+7)= 0;
            B(podminka_pocet, (l_index-1)*12+8)= 0;
            B(podminka_pocet, (l_index-1)*12+9)= 0;
            B(podminka_pocet, (l_index-1)*12+10)= 0;
            B(podminka_pocet, (l_index-1)*12+11)= 0;
            B(podminka_pocet, (l_index-1)*12+12)= 0;
            B(podminka_pocet+1, (p_index-1)*12+7)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2)*((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2);
            B(podminka_pocet+1, (p_index-1)*12+8)= ((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2)*((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet+1, (p_index-1)*12+9)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2)*((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet+1, (p_index-1)*12+10)= ((Bpomocna(p_index,7)+Bpomocna(p_index,1))/2);
            B(podminka_pocet+1, (p_index-1)*12+11)= ((Bpomocna(p_index,8)+Bpomocna(p_index,2))/2);
            B(podminka_pocet+1, (p_index-1)*12+12)= 1;
            B(podminka_pocet+1, (p_index-1)*12+1)= 0;
            B(podminka_pocet+1, (p_index-1)*12+2)= 0;
            B(podminka_pocet+1, (p_index-1)*12+3)= 0;
            B(podminka_pocet+1, (p_index-1)*12+4)= 0;
            B(podminka_pocet+1, (p_index-1)*12+5)= 0;
            B(podminka_pocet+1, (p_index-1)*12+6)= 0;
            B(podminka_pocet+1, (l_index-1)*12+7)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2)*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2);
            B(podminka_pocet+1, (l_index-1)*12+8)= -1*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2)*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet+1, (l_index-1)*12+9)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2)*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet+1, (l_index-1)*12+10)= -1*((Bpomocna(l_index,3)+Bpomocna(l_index,5))/2);
            B(podminka_pocet+1, (l_index-1)*12+11)= -1*((Bpomocna(l_index,4)+Bpomocna(l_index,6))/2);
            B(podminka_pocet+1, (l_index-1)*12+12)= -1;
            B(podminka_pocet+1, (l_index-1)*12+1)= 0;
            B(podminka_pocet+1, (l_index-1)*12+2)= 0;
            B(podminka_pocet+1, (l_index-1)*12+3)= 0;
            B(podminka_pocet+1, (l_index-1)*12+4)= 0;
            B(podminka_pocet+1, (l_index-1)*12+5)= 0;
            B(podminka_pocet+1, (l_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
        end
    end
end          


for j=1:(size(klad,2))
    for i=1:(size(klad,1)-1)
        if klad(i,j)~=0 && klad(i+1,j)~=0   
            h=klad(i,j);
            d=klad(i+1,j);
            h_index=find(jmena_ib==h);
            d_index=find(jmena_ib==d);
            B(podminka_pocet, (h_index-1)*12+1)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2)*((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2);
            B(podminka_pocet, (h_index-1)*12+2)= ((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2)*((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet, (h_index-1)*12+3)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2)*((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet, (h_index-1)*12+4)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2);
            B(podminka_pocet, (h_index-1)*12+5)= ((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet, (h_index-1)*12+6)= 1;
            B(podminka_pocet, (h_index-1)*12+7)= 0;
            B(podminka_pocet, (h_index-1)*12+8)= 0;
            B(podminka_pocet, (h_index-1)*12+9)= 0;
            B(podminka_pocet, (h_index-1)*12+10)= 0;
            B(podminka_pocet, (h_index-1)*12+11)= 0;
            B(podminka_pocet, (h_index-1)*12+12)= 0;
            B(podminka_pocet, (d_index-1)*12+1)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2)*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2);
            B(podminka_pocet, (d_index-1)*12+2)= -1*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2)*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet, (d_index-1)*12+3)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2)*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet, (d_index-1)*12+4)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2);
            B(podminka_pocet, (d_index-1)*12+5)= -1*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet, (d_index-1)*12+6)= -1;
            B(podminka_pocet, (d_index-1)*12+7)= 0;
            B(podminka_pocet, (d_index-1)*12+8)= 0;
            B(podminka_pocet, (d_index-1)*12+9)= 0;
            B(podminka_pocet, (d_index-1)*12+10)= 0;
            B(podminka_pocet, (d_index-1)*12+11)= 0;
            B(podminka_pocet, (d_index-1)*12+12)= 0;
            B(podminka_pocet+1, (h_index-1)*12+7)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2)*((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2);
            B(podminka_pocet+1, (h_index-1)*12+8)= ((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2)*((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet+1, (h_index-1)*12+9)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2)*((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet+1, (h_index-1)*12+10)= ((Bpomocna(h_index,7)+Bpomocna(h_index,5))/2);
            B(podminka_pocet+1, (h_index-1)*12+11)= ((Bpomocna(h_index,8)+Bpomocna(h_index,6))/2);
            B(podminka_pocet+1, (h_index-1)*12+12)= 1;
            B(podminka_pocet+1, (h_index-1)*12+1)= 0;
            B(podminka_pocet+1, (h_index-1)*12+2)= 0;
            B(podminka_pocet+1, (h_index-1)*12+3)= 0;
            B(podminka_pocet+1, (h_index-1)*12+4)= 0;
            B(podminka_pocet+1, (h_index-1)*12+5)= 0;
            B(podminka_pocet+1, (h_index-1)*12+6)= 0;
            B(podminka_pocet+1, (d_index-1)*12+7)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2)*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2);
            B(podminka_pocet+1, (d_index-1)*12+8)= -1*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2)*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet+1, (d_index-1)*12+9)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2)*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet+1, (d_index-1)*12+10)= -1*((Bpomocna(d_index,3)+Bpomocna(d_index,1))/2);
            B(podminka_pocet+1, (d_index-1)*12+11)= -1*((Bpomocna(d_index,4)+Bpomocna(d_index,2))/2);
            B(podminka_pocet+1, (d_index-1)*12+12)= -1;
            B(podminka_pocet+1, (d_index-1)*12+1)= 0;
            B(podminka_pocet+1, (d_index-1)*12+2)= 0;
            B(podminka_pocet+1, (d_index-1)*12+3)= 0;
            B(podminka_pocet+1, (d_index-1)*12+4)= 0;
            B(podminka_pocet+1, (d_index-1)*12+5)= 0;
            B(podminka_pocet+1, (d_index-1)*12+6)= 0;
            podminka_pocet=podminka_pocet+2;
        end
    end
end          



%Vypocty

matice=[A'*A, B'; B, zeros(size(B,1),size(B,1))];
vektor=[A'*vektor_l; zeros(size(B,1),1)];

nezname= inv(matice)*vektor;
x=nezname(1:size(A,2));
v=A*x-vektor_l;
m_nula= sqrt(v'*v/(size(A,1)+size(B,1)-size(A,2)));         







% export IB
    dir_ibtrans  = '_data\_output\ib\';
    for i=1:size(names_ib,1)
        fid = fopen([dir_ibtrans,names_ib(i).name],'w');
        points_ib = (load([dir_ib,names_ib(i).name]));
        localx=x( (12*(i-1)+1):(12*i))
        for j=1:size(points_ib,1)
            vysledek_x=localx(1)*points_ib(j,1)*points_ib(j,1)+localx(2)*points_ib(j,2)*points_ib(j,2)+localx(3)*points_ib(j,1)*points_ib(j,2)+localx(4)*points_ib(j,1)+localx(5)*points_ib(j,2)+localx(6)
            vysledek_y=localx(7)*points_ib(j,1)*points_ib(j,1)+localx(8)*points_ib(j,2)*points_ib(j,2)+localx(9)*points_ib(j,1)*points_ib(j,2)+localx(10)*points_ib(j,1)+localx(11)*points_ib(j,2)+localx(12)
            fprintf(fid,'%.6f\t',points_ib(j,1));
            fprintf(fid,'%.6f\t',points_ib(j,2));
            fprintf(fid,'%.6f\t',vysledek_x);
            fprintf(fid,'%.6f\n',vysledek_y);              
        end
        fclose(fid);
    end

    
 % export ramu    
 dir_rohytrans  = '_data\_output\mezi\';
    for i=1:size(names_rohy,1)
        fid = fopen([dir_rohytrans,names_rohy(i).name],'w');
        points_rohy = (load([dir_rohy,names_rohy(i).name]));
        localx=x( (12*(i-1)+1):(12*i));
        for j=0:20
            pred_x1=points_rohy(1,1)+j/20*(points_rohy(2,1)-points_rohy(1,1));
            pred_y1=points_rohy(1,2)+j/20*(points_rohy(2,2)-points_rohy(1,2));
            vysledek_x1=localx(1)*pred_x1*pred_x1+localx(2)*pred_y1*pred_y1+localx(3)*pred_x1*pred_y1+localx(4)*pred_x1+localx(5)*pred_y1+localx(6);
            vysledek_y1=localx(7)*pred_x1*pred_x1+localx(8)*pred_y1*pred_y1+localx(9)*pred_x1*pred_y1+localx(10)*pred_x1+localx(11)*pred_y1+localx(12);
            fprintf(fid,'%.6f\t',vysledek_x1);
            fprintf(fid,'%.6f\n',vysledek_y1);     
        end    
        for j=0:20
            pred_x2=points_rohy(2,1)+j/20*(points_rohy(3,1)-points_rohy(2,1));
            pred_y2=points_rohy(2,2)+j/20*(points_rohy(3,2)-points_rohy(2,2)); 
            vysledek_x2=localx(1)*pred_x2*pred_x2+localx(2)*pred_y2*pred_y2+localx(3)*pred_x2*pred_y2+localx(4)*pred_x2+localx(5)*pred_y2+localx(6);
            vysledek_y2=localx(7)*pred_x2*pred_x2+localx(8)*pred_y2*pred_y2+localx(9)*pred_x2*pred_y2+localx(10)*pred_x2+localx(11)*pred_y2+localx(12);
            fprintf(fid,'%.6f\t',vysledek_x2);
            fprintf(fid,'%.6f\n',vysledek_y2);         
        end    
        for j=0:20
            pred_x3=points_rohy(3,1)+j/20*(points_rohy(4,1)-points_rohy(3,1));
            pred_y3=points_rohy(3,2)+j/20*(points_rohy(4,2)-points_rohy(3,2));
            vysledek_x3=localx(1)*pred_x3*pred_x3+localx(2)*pred_y3*pred_y3+localx(3)*pred_x3*pred_y3+localx(4)*pred_x3+localx(5)*pred_y3+localx(6);
            vysledek_y3=localx(7)*pred_x3*pred_x3+localx(8)*pred_y3*pred_y3+localx(9)*pred_x3*pred_y3+localx(10)*pred_x3+localx(11)*pred_y3+localx(12);
            fprintf(fid,'%.6f\t',vysledek_x3);
            fprintf(fid,'%.6f\n',vysledek_y3);          
        end    
        for j=0:20
            pred_x4=points_rohy(4,1)+j/20*(points_rohy(1,1)-points_rohy(4,1));
            pred_y4=points_rohy(4,2)+j/20*(points_rohy(1,2)-points_rohy(4,2));  
            vysledek_x4=localx(1)*pred_x4*pred_x4+localx(2)*pred_y4*pred_y4+localx(3)*pred_x4*pred_y4+localx(4)*pred_x4+localx(5)*pred_y4+localx(6);
            vysledek_y4=localx(7)*pred_x4*pred_x4+localx(8)*pred_y4*pred_y4+localx(9)*pred_x4*pred_y4+localx(10)*pred_x4+localx(11)*pred_y4+localx(12);
            fprintf(fid,'%.6f\t',vysledek_x4);
            fprintf(fid,'%.6f\n',vysledek_y4);     
        end
        fclose(fid);
    end





