function [fib,lb,Ab] = kivioj(fi,l,A)
a=6378137;
e2=0.00669438002290;
s=25000;
A=66*pi/180;
fi=0.90805;
l=0.32100;
n=round(s/1000);
ds=s/n
for i=1:n
    M(i)=a*(1-e2)/sqrt((1-e2*(sin(fi(i))^2))^3);
    N(i)=a/sqrt(1-e2*(sin(fi(i)))^2);
    dfi(i)=cos(A(i))*ds/M(i);
    dA(i)=sin(A(i))*tan(fi(i))*ds/N(i);
    fim(i)=fi(i)+(1/2)*dfi(i);
    Am(i)=A(i)+(1/2)*dA(i);
    Mm(i)=a*(1-e2)/sqrt((1-e2*(sin(fim(i))^2))^3);
    
    Nm(i)=a/sqrt(1-e2*(sin(fim(i)))^2);
    dfim(i)=ds*cos(Am(i))/Mm(i);
    dlm(i)=ds*sin(Am(i))/Nm(i)*cos(fim(i));
    dAm(i)=ds*sin(Am(i))*tan(fim(i))/Nm(i);
    fi(i+1)=fi(i)+dfim(i);
    l(i+1)=l(i)+dlm(i);
    A(i+1)=A(i)+dAm(i);
end
fib=fi(n)
lb=l(n)
Ab=A(n)
end