function [ b_fi, b_la, A_BA ] = kivioj(fi, la, A, s, a, e2)

n = round(s/1000); 
ds = s/n; 
for i = 1:n   
    M = M_p(fi,a,e2);
    N = N_p(fi,a,e2);
    dfi_p = (ds*cos(A))/M;
    dA_p = (ds*sin(A)*tan(fi))/N;
    fi_m = fi + dfi_p/2;
    A_m = A + dA_p/2;
    M_m = M_p(fi_m,a,e2);
    N_m = N_p(fi_m,a,e2);
    dfi = (ds*cos(A_m))/M_m;
    dla = (ds*sin(A_m))/(N_m*cos(fi_m));
    dA = (ds*sin(A_m)*tan(fi_m))/N_m;
    
    fi = fi + dfi;
    la = la + dla;
    A = A + dA;
end;

b_fi = fi;
b_la = la;

if A < pi 
    A_BA = A + pi;
else
    A_BA = A - pi;
end;

end