duljina luka

function [G] = duluka(fi)

%%
  a = 6378137.0; % earth semimajor axis in meters
  finv = 298.257223563; % reciprocal flattening
  f=1/finv;
  b=a*(1-f);
  e=sqrt((a^2-b^2)/a^2);
  e2=e^2;
%%
A=1+3/4*e2+45/64*e2*e2;
B=3/4*e2+15/16*e2*e2;
C=15/64*e2*e2;

for i=1:1:max(size(fi))
  if i>=2

    G=a*(1-e2)*(A(fi(i)-fi(i-1))-B/2*(sin*2*fi(i)-sin*2*fi(i-1))+C/4*(sin*4*fi(i)-sin*4*fi(i-1)));
  end
end