%
%  Les gögn
%
load glm1.txt;
%
%
%
nn=rows(glm1);
einn=ones(nn,1);
x=[einn,glm1(:,2)];
y=glm1(:,1);
%
%  skilgreini fall og byrjunargildi
%
b=inv(x'*x)*x'*y
%
% b fengið með venjulegri regression
%

for ii=1:1000
xib=x*b;
ynytt=y./xib;
xnytt=[x(:,1)./xib, x(:,2)./xib];
bnytt=inv(xnytt'*xnytt)*xnytt'*ynytt
if (abs(b-bnytt)<0.000001), break, end
b=bnytt
end
%

