ghreggo

function c = newtonbruun(x,y)
%
%   newtonbruun(x,y)
%   Fallið tekur inn tvo jafnstóra vigra x og y og
%   skilar stuðlum c í Newton formi brúunarmargliðunnar
%   sem fer í gegnum punktana (x,y).
%
% Toluleg greining, vor 2016: Verkefni I, 2 hluti
% Hofundar: Ólafur Ívar Baldvinsson, Hafsteinn Einarsson
% og Agnes Linnet
% Dags:

n=length(x);

m=1; %m vigur sem inniheldur margfeldni stuðlanna

for i=2:n
    if x(i)==x(i-1)
        m(end)=m(end)+1;
    else
        m=[m 1];
    end
end

%Búum til mismunakvótatöflu
%Í fyrsta dálkinn þurfum við fallgildi y

Mkt=zeros(n);
yf=zeros(1,n);
l=length(m);

yf(1:m(1))=y(1);
pick(1:m(1))=1;
for i=2:l
    ind=sum(m(1:i-1))+1;
    yf(ind:ind+m(i)-1)=y(ind);
    pick(ind:n)=ind;
end

Mkt(:,1)=yf;

%fyllum nú út í restina af Mkt
k=n-1;

for j=2:n
    for i=1:k
        if x(i)==x(i+j-1)
            Mkt(i,j)=y(pick(i)+j-1)/factorial(j-1);
        else
            Mkt(i,j)=(Mkt(i+j-(n-k),j-1)-Mkt(i,j-1))/(x(i+j-1)-x(i));
        end
    end
    k=k-1;
end

%Tökum nú fyrstu línuna í Mtk og skilum sem stuðlafylkinu c

c=zeros(1,n);

for i=1:n
    c(i)=Mkt(1,i);
end

end