function [body_cary] = nulova_cara

% Ze souboru nc.txt,
% ve kterem jsou pod sebou body mezi nimiz povede 
% nulova cara (vzdy tedy po dvojici),
% vytvori soubor, kde jsou spocteny body cary
% cislovane od parametru fce
%
%      body_cary=nulova_cara(cislovat_od)

format bank
load nc.txt
pocet_bodu=size(nc);
pocet_bodu=pocet_bodu(1,1);
for i=1:2:pocet_bodu
    y1=nc(i,2);
    x1=nc(i,3);
    z1=nc(i,4);
    y2=nc(i+1,2);
    x2=nc(i+1,3);
    z2=nc(i+1,4);
    
    s12=((y1-y2)^2+(x1-x2)^2)^(1/2);
    v1m=-z1;
    v12=z2-z1;
    y12=y2-y1;
    x12=x2-x1;
    
    s1m=s12*v1m/v12;
    ym=y1+s1m*y12/s12;
    xm=x1+s1m*x12/s12;
    
    body_cary((i+1)/2,1)=(i+1)/2;
    body_cary((i+1)/2,2)=ym;
    body_cary((i+1)/2,3)=xm;
end

save body_cary.txt body_cary -ascii;