function c = newtonbruunn(x,y)
if length(x)~=length(y)
    error('Vigrarnir eru ekki af sömu stærð')
end
n = length(x);
fylki = zeros(n,n+1); %Búum til nxn+1 fylki
fylki(:,1) = x';
fylki(:,2) = y'; %Byltum vigrunum x og y svo þetta lýti út eins og mismunakvótataflan sem við gerðum í skiladæmum
lengdy = n-1;
i = 3;
j = 1;
while i<=n+1;
    while  j <= lengdy;
        fylki(j,i) = (fylki(j+1,i-1)-fylki(j,i-1))/(fylki(j+j)-fylki(j));
        j = j+1;
    end
    %While-lykkjum okkur í gegnum fylkið, stak fyrir stak,reiknum og röðum í
    %zeros fylkið sem við bjuggum til í byrjun
    lengdy = lengdy-1;
    i = i+1;
end
c = fylki(1,(2:end)); %prentum út efstu línuna í fylkinu nema stak (1,1)